adv_simpson.m

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

function Adv_Simpson
%复化Simpson公式
format long;
syms x;
f=input('被积函数式f(x)=');
a=input('积分下限a=');
b=input('积分下限b=');
n=input('节点数n=');

h=(b-a)/(n-1);
sum0=0;
sum1=0;
start=subs(f,x,a);
bottom=subs(f,x,b);
F=zeros(n,2);
for i=1:n-1
    F(i,1)=subs(f,x,(a+i*h));
    sum0=sum0+F(i,1);
end
for i=1:n
    F(i,2)=subs(f,x,(a+(k-0.5)*h));
    sum1=sum1+F(i,2);
end
S=h*(start+4*sum1+2*sum0+bottom)/6;
I=S