romberg4.m

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

function q = romberg4(a,b,n,randa)
h=(b-a)/n;
q=0;
x=a;
for i=1:n
    x=x+(n-1)*h;
    x0=x;
    x1=x+h/8;
    x2=x+2*h/8;
    x3=x+3*h/8;
    x4=x+4*h/8;
    x5=x+5*h/8;
    x6=x+6*h/8;
    x7=x+7*h/8;
    x8=x+8*h/8;
    zz=h/5670*(441*fx4(x0,randa)+2048*fx4(x1,randa)+736*fx4(x2,randa)+2048*fx4(x3,randa)+2036*fx4(x4,randa)+2048*fx4(x5,randa)+736*fx4(x6,randa)+2048*fx4(x7,randa)+441*fx4(x8,randa));
    q=q+zz;
end