ep3_p1.m

来自「此代码展示了Euler计算方法与差分方法」· M 代码 · 共 24 行

%  Ep3_p1: << 实验三 >>  Euler法 (降阶式)
%  Designed by FGH

n= 266;
H= 120; Vw= 450; Ve= 90;
lanmuda= Ve/Vw;
h= H/n;

clear x p;
y= 0:h:H;
% (3.17)
x(1)= 0; p(1)= 0;
for k= 1:n
   % (3.15)
   x(k+1)= x(k) + h*p(k);
   % (3.16)
   p(k+1)= p(k) + h*lanmuda*sqrt(1+p(k)^2)/(H-y(k));
end

y
x
p
L= x(n+1)
T= x(n+1)/Ve