%Integration,using adaptive Simpson methom(quad8)
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
%算得值總體上(特別離原點較遠(yuǎn)位置)與真實值還是比較接近的，一下舉例說明：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