数值分析（代码）-梯形公式、辛普森公式、复合梯形公式、复合辛普森公式

习题5.4
运行代码：
main.h:
clc;clear;format long;
h=1/6;
%h=1/8;%例5.4
N=(1-0)/h;
x=0:h:1;
s_t=f(x(1));
for n=1:N-1
    s_t=s_t+2*f(x(n)+h);
end
s_t=s_t+f(x(N+1));
s_t=s_t*h/2;
s_t
s=0;
for n=1:4
    s=s+2*f((2*n-1)/8)+f(n/4);
end
s=2*s+f(0)-f(1);
s=s/24

f.m:
function f=f(x)
if x==0
    f=1;
else 
   % f=sin(x)/x;
     f=x/log(x+1);
end

运行结果： 

这一章节的其他代码见资源