euler.m

来自「微分的数值计算综合了算法在导弹中的应用」· M 代码 · 共 22 行

%euler.m
function [x,y]=euler();
input('请按下列步骤输入初始值!');
x0=input('x0:');
y0=input('y0:');
h=input('步长:');
e=input('最大范围:');

i=1:h:e;
n=max(size(i));
x(1)=x0;
y(1)=y0;
for i=1:n-1
x(i+1)=x(i)+h*(x(i)-4*y(i)*sqrt(abs(x(i)*y(i))));
y(i+1)=y(i)+h*(-y(i)+4*x(i)*sqrt(abs(x(i)*y(i))));
end
hold on
plot(x,y)
title('欧拉方法')