% This is the main program.priority_vector1 = [9 4 2 1];priority_vector2 = [9 4 2 1];alpha1 = priority_vector1*(priority_vector1)';alpha2 = priority_vector2*(priority_vector2)';size_of_in_phase = length(priority_vector1);size_of_q_phase = length(priority_vector2);err_vector1=zeros(30,size_of_in_phase);err_vector2=zeros(30,size_of_q_phase);     for snr=-5:1:24	snr_no=10^(snr/10);	noise_qty = sqrt( (alpha1+alpha2) /(snr_no));	for i=1:size_of_in_phase		err_vector1(snr+6,i)=recurse_other(i,priority_vector1,noise_qty);	end;end;for snr=-5:1:24        snr_no=10^(snr/10);        noise_qty = sqrt( (alpha1+alpha2) /(snr_no));        for i=1:size_of_q_phase                err_vector2(snr+6,i)=recurse_other(i,priority_vector2,noise_qty);        end;end;figure;%%% Plot the vectors nowfor i=1:size_of_in_phase             k=-5:1:24;        semilogy(k,err_vector1(2+(k+4),i));        hold on;end;for i=1:size_of_q_phase	k=-5:1:24;	semilogy(k,err_vector2(1+(k+5),i));	hold on;end;grid;hold off;xlabel('Carrier to noise ratio E_{s}/N_{0} [dB]');ylabel('Average bit error probability');