comsimpon.m

来自「基于极大嫡谱估计准则对变形数据进行预测 基于MATLAB平台编制了相应程序」· M 代码 · 共 12 行

function y=comsimpon(f,a,b,n)
s=evalf(f,a)-evalf(f,b);
h=(b-a)/(2*n);
x=a;
for i=1:2:(2*n-1)
    r=x+h;
    s=s+4*evalf(f,x);
    x=x+h;
    s=s+2*evalf(f,x);
end
s=s*h/3.0;
y=s;