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)