I want to translate my matlab code fortran?

theta=0:10:60;  % Note: if you reduce the step size you will get a smoother plot
rho=0.465;
%d=25; wrong
%k=0.99*(1+(1i*0.66));
k=real(0.97+0.6402i);
poftheta=sqrt((k.^2)-(cos(theta).^2));
y=(rho*sin(theta))./(poftheta);
roftheta=((y-1)./(y+1));
             % ^---- missed dot
r0=-20*log10(mod(roftheta));
            % ^----- here
grid on
plot(theta,r0)

[madhan ravi]

Duplicate of https://www.mathworks.com/matlabcentral/answers/437069-hi-i-am-trying-to-plot-the-following-equation-its-showing-error-in-mod-function-kindly-clarify-me

[Image Analyst]

What do you mean by "fortran" in your subject line and tag. Do you want to translate this code INTO the FORTRAN programming language? If so, why, and it looks simple enough so why can't you do it?

Answer 1

[KSSV]

Read the documentation of mod. It takes two inputs, you have to decide what should be the input. For demo, below I have taken mod(n,2).

theta=0:10:60;  % Note: if you reduce the step size you will get a smoother plot
rho=0.465;
%d=25; wrong
%k=0.99*(1+(1i*0.66));
k=real(0.97+0.6402i);
poftheta=sqrt((k.^2)-(cos(theta).^2));
y=(rho*sin(theta))./(poftheta);
roftheta=((y-1)./(y+1));
             % ^---- missed dot
r0=-20*log10(mod(roftheta,2));
            % ^----- here
grid on
plot(theta,r0)

[Anu G]

Thank you sir