euler.m

来自「基于Matlab的计算方法中常用的算法程序」· M 代码 · 共 18 行

function s=Euler(fun,x0,xn,y0,n)
  % 用Euler法计算常微分方程初值问题
  % x0为初值条件y(x0)=y0中的x0，y0为初值条件中的y0
  % xn为x取值区间的最后一个节点的横坐标值
  % n为把区间分成的等份数
  if nargin<5
      error
      return
  end
  h=(xn-x0)/n;
  for i=1:n
      yp=y0+h*feval(fun,x0,y0);
      x0=x0+h;
      yc=y0+h*feval(fun,x0,yp);
      y0=(yp+yc)/2;
  end
  s=y0;
  return