N = 50 ;
B = rand ; % a constant
a =linspace(0,4,N) ;
phi = linspace(0,2,N)*pi/180 ;
[A, Phi] = meshgrid(a,phi) ;
R = A.*sin(Phi/2)./(B*cos(Phi/2).^2) ;
surf(A,Phi,R)
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!