fhtrap.m

来自「数值分析中的插值计算」· M 代码 · 共 23 行

function fhtrap(a,b,n)
%复化梯形公式
%Input   f is the function
%        a,b是积分下限与上限
%        n是节点个数
%Output   
%         T 求积值
syms x;
f=input('f=');
n=n-1;
h=(b-a)/n;
sum=0;
F=zeros(n,1);
F(n+1)=subs(f,x,a);
F(n)=subs(f,x,b);
for k=1:n-1
    F(k)=subs(f,x,(a+k*h));
    sum=sum+F(k);
end
T=h*(F(n+1)+F(n)+2*sum)/2;
F
T