ketesi.m

来自「数值分析方法,主要里面涵盖了计算方法方面matlab程序！」· M 代码 · 共 20 行

%function KeTeSi(a,b) %a为积分下限,b为积分上限。
clear
clc
eps=0.5E-6;
a=0;
b=1;
n=1;
s1=0;
while(1)
    x=a:(b-a)/(4*n):b;
    m=4*n+1;
    h=(b-a)/n;
    s=(h/90)*(7*f(a)+2*sum(f(x(3:2:m-2)))+4*sum(f(x(2:2:m-1)))+7*f(b));
    fprintf('s(%d)=%f\n',n,s);
    if abs(s1-s)<eps
        break;
    end
    s1=s;
    n=n*2;
end