trapezia.m

数值分析常用几种方法

function S=trapezia(f,a,b,f0,h)
%f为函数f(x,y);
%a,b为取值范围的左右端点；
%f0为函数初值y(0)，h为步长,n为弟迭代步数；
%feval函数实现输入变量按照指定函数计算的函数值.
n=(b-a)/h;
T=zeros(1,n+1);
Y=zeros(1,n+1);
T=a:h:b;
Y(1)=f0;
for j=1:n
    z0=Y(j)+h*feval('f',T(j),Y(j));
    Y(j+1)=Y(j)+h*(feval('f',T(j),Y(j))+feval('f',T(j+1),z0))/2;
end
t=[0:.1:5];
y=1./exp(10*t);
figure(1)
plot(T,Y,'*',T,Y,'r-',t,y,'k-');
hold on
T=[T' Y']