改进欧拉格式的matlab实现

改进欧拉格式的matlab实现

Funtion E4=Euler_4(fun,x0,y0,xN,N)

%改进的Euler公式，其中，

%fun为一阶微分方程的函数

%x0,y0为初始条件

%xN为取值范围的一个断点

%h为区间步长

%N为区间的个数

%x为Xn构成的向量

%y为Yn构成的向量

x=zeros(1,N+1);y=zeros(1,N+1);

x(1)=x0;y(1)=y0;

h=(xN-x0)/N;

for n=1:N

x(n+1)=x(n)+h;

y(n+1)=y(n)+h/2*(feval(fun,x(n),y(n))+ feval(fun,x(n+1),z0); end

T=[x’,y’]

建立一个M-文件，其内容为

Funtion z=f(x,y)

z=x+y;

运行窗口中输入>>Euler_4(‘f’,0,1,0.1,5)