梯形法的matlab程序

梯形法的matlab程序functiontixing(f,x0,y0,a,b,n)f=inline(f);h=(b-a)/n;x=a:h:b;y(1)=y0;fori=1:ny(i+1)=y(i)+h*f(x(i),y(i));fork=1:2*ny(i+1)=y(i)+h*(f(x(i),y(i))+f(x(i+1),y(i+1)))/2;endendy