what is wrong with this code?

 
functiony=pdf_rician(k,c,sigma)
 
%% pdf of process of various values of system parameters 
 
for x=0:0.0625:10
temp=0; 
pdf_rici=0;
 
for p=0:10
    for n=0:(p+1)
        %    closed term expansion 
        b=nchoosek(p+0.5,n);  
        
        % gamma term calculation 
        f=factorial(p);
        g=(gamma(p+1))*(gamma(c));
        g_ma= 2*sqrt(2)*sqrt(pi)*((x.^(2*p+1))*(k.^c))/((sigma.^(p+0.5))*g*f);   
        
        % bessel function term calculation  
        z=2*sqrt(k*(sigma+x.^2)/sigma);
        K = besselk(p-n+c,z); 
        
        % Exponential term calculation
        exponential=exp(-x.^2/sigma); 
        % Last term 
        q=(p-n+c)/2;
        last_term=((k*sigma/(sigma+(x.^2))).^q);
        
        % Final pdf function 
        
        temp1= b*g_ma*K*exponential*last_term;
        temp=temp+temp1;
    end
    pdf_rici=pdf_rici+temp;
end
 
y(1,x/(.0625)+1)=pdf_rici; 
 
end 

[Walter Roberson]

What leads you to say there is a problem with the code? Is there an error message?

[Atrolita afra]

yes i could not get the figure and there was an error message syaing,"" Matrix dimensions must agree.

Answer 1

[Jamal Nasir]

the variable temp1 is not single value. it is a vector so need to check the equation

[Atrolita afra]

here is the equation,that i am trying to plot