what is the error in this program

for k=5:5:50

z=zeros(1,length(k));

Pr1=(10*0.75*0.75*(exp(-0.984*k/cos(pi/36)))*0.01*cos(pi/36)/((2*pi*k.^2)*(1-cos(pi/3))));

%Pr_dbm1(k)=10*log10(Pr1)

z(k)=10*log10(Pr1);

plot(k,z)

xlabel('distance')

ylabel('Power')

title('Rcvd Power vs Distance')

%Pr1

end

Answer 1

[Stephen23]

k = 5:5:50;
Pr1 = (10*0.75*0.75*(exp(-0.984*k./cos(pi/36))) .* 0.01 .* cos(pi/36) ./ ((2*pi*k.^2)*(1-cos(pi/3))));
z=10.*log10(Pr1);
plot(k,z,'-*')

[sanjiv kumar]

thank you for helping me............

Answer 2

[Alex Mcaulley]

There is no error in your code, the problem is that it doesn't do whatt you want. I guess you want something like this (without loops):

k=5:5:50;
Pr1=(10*0.75*0.75*(exp(-0.984*k/cos(pi/36)))*0.01*cos(pi/36)./((2*pi*k.^2)*(1-cos(pi/3))));
%Pr_dbm1(k)=10*log10(Pr1)
z=10*log10(Pr1);
plot(k,z)
xlabel('distance')
ylabel('Power')
title('Rcvd Power vs Distance')