func5p2.m

来自「数值分析最常用的四十种算法」· M 代码 · 共 11 行

%Integration,using Simpson methom
x=linspace(0.1,3.0,30)
y=linspace(0.1,3.0,30)
for i=1:30
    z=sub5p2(x(i))
    disp(z)
    y(i)=(sub5p2(x(i))+4*sub5p2((x(i)+4)/2)+sub5p2(4))/((2*pi)^(1/2)*6)*(4-x(i))
end
%算得值总体上(特别离原点较远位置)与真实值还是比较接近的，一下举例说明：x=0.1,y=0.3850,true
%value:0.4602;x=1,y=0.1561;true value:0.1587;x=3,y=0.0013,true
%value:0.00135