% This does the plotting for 2/4/8-PSKinput_angle_vector = [pi/2 pi/3 pi/15 pi/40 ]'temp_size=length(input_angle_vector);angle_vector=input_angle_vector(2:temp_size);angle_vector(1)=(pi/2)-angle_vector(1);size_of_constellation = length(angle_vector);err_vector=zeros(30,size_of_constellation);for snr=-5:1:24    snr;	snr_no=10^(snr/10);	noise_qty = sqrt(1/(snr_no));	for i=1:size_of_constellation+1        i;		err_vector(snr+6,i)=recurse_other(i,angle_vector,noise_qty);	end;end;err_vectorfor i=1:size_of_constellation+1	k=-5:1:24;	semilogy(k,err_vector(k+6,i));	hold on;end;grid;hold off;xlabel('Carrier to noise ratio E_{s}/N_{0} [dB]');ylabel('Average bit error probability');