adv_trap.m

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

function Adv_Trap
%复化梯形公式 
format long;
syms x;
f=input('被积函数式f(x)=');
a=input('积分下限a=');
b=input('积分下限b=');
n=input('节点数n=');
h=(b-a)/(n-1);
sum=0;
F=zeros(n+1,1);
F(1)=subs(f,x,a);
F(n+1)=subs(f,x,b);
for i=2:n
    F(i)=subs(f,x,(a+(i-1)*h));
    sum=sum+F(i);
end
F
T=(F(1)+2*sum+F(n+1))*h/2;