simpson.m

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

function simpson(a,b,M)
%Simpson求积公式
%Input   f is the function
%        a,b是积分下限与上限
%        M是f的四阶导在(a,b)间的最大值 
%Output   
%         S 求积值
%         R 余项值
syms x;
f=input('f=');
S=(b-a)*(subs(f,x,a)+4*subs(f,x,(a+b)/2)+subs(f,x,b))/6;
R=-M*(b-a)^5/2880;
fprintf('S=%12.5f\n',S)
fprintf('R=%12.5f\n',R)