ep2_f0.m

此代码展示了Euler计算方法与差分方法

%  Ep2_f0: << 实验二 >>  算式(2.13)所对应的函数
%  Designed by FGH

function y= ep2_f0(r, e)
   a= 30.6;
   b= 29.6;
   R= 30;
   % formula (2.13)
   y= a*b*sqrt( r.^2 + R^2*(1-r.^2)*((cos(e)/a)^2+(sin(e)/b)^2) );
   % formula (2.3)
   %y= a*b*r*sqrt( 1 + (R*r).^2*((cos(e)/a)^2+(sin(e)/b)^2)/(1-r.^2) );