The MATLAB code is written as follows:
syms delm phim kpm nup y1phi=(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)
Best Answer