simpsonint4.m

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

 function I4=int(m)
%simpson4
x=0:1/4:1;
for i=1:5,
y(i)=x(i)/(4+x(i)^2);
end
I4=1/6*(1/2)*(y(1)+y(5)+2*y(3)+4*(y(2)+y(4)))