The MATLAB code is written as follows:

syms delm phim kpm nup y1

phi=(pi/180)*39;

del=(pi/180)*26;

beta=0.5;

phim=phi*beta*(1-(y1));

delm=((phim)/(phi))*del;

nup=((asin(sin(delm))/(sin(phim)))+delm)/2;

kpm=(cos(delm)/(1-sin(phim)))*(cos(delm)+sqrt((sin(phim))^2-...

(sin(delm))^2))*exp(2*nup*tan(phim));

fu=kpm*cos(delm)*y1*y1;

int(fu,0.1,0.2)

More Akshay Pratap Singh's questions See All
Similar questions and discussions