func5p1.asv

数值分析最常用的四十种算法

%Integration,using N-L methom
I1=4*(atan(1)-atan(0))
%Trapz Rule
I2=1/2*(sub5p1(1)+sub5p1(0))
%Simpson Rule
I3=1/6*(sub5p1(1)+4*sub5p1(1/2)+sub5p1(0))
%Composite Trapezoidal Rule
x=linspace(0,1,100)
I4=0
for i=1:99
I4=I4+sub5p1(x(i))+sub5p1(x(i+1))
end
I4=I4/2/100
%3.1416,3, 3.1333,3.1102