simpsonint.m

来自「基于matlab的数值计算试验程序。equ_direct」· M 代码 · 共 7 行

function I=int(n)
%simpson
x=0:1/8:1;
for i=1:9,
y(i)=x(i)/(4+x(i)^2);
end
I=1/6*(2/8)*(y(1)+y(9)+2*(y(3)+y(5)+y(7))+4*(y(2)+y(4)+y(6)+y(8)))