clear
tic
f=@(x,y) y/x-y^2/x;
h=0.01;
y0=0.5;
x0=1;

while x0<=1.5
      y1=y0+h*f(x0,y0);
      x1=x0+h;
      y2=y0+h*(f(x0,y0)+f(x1,y1))/2;
      %梯形方法迭代公式
      y0=y1;
      y1=y2;
      x0=x1;
end

disp('梯形方法解：');y0
disp('解析解：');
g=@(x) x/(1+x);y=g(1.5)
disp('计算用时：');
toc
在这里只需把f(x,y)和g(x)修改一下即可，此版本是基于R2008a写的，不适用于低于6.5版本一下。