Question

请问如何用matlab完成这道题？

Answer 1

clear
clc
f=@(x)4/(1+cos(x))
a=0;b=1;h=1;esp=10^-5;
fa=f(a);fb=f(b);
S1=h/6*(fa+fb+4*f(a/2+b/2));
h=h/2;
S2=h/6*(fa+4*f(a+h/2)+2*f(a+h)+4*f(a+3*h/2)+fb);
while abs(S1-S2)>esp
S1=S2;
h=h/2;
s=0;x=a;
N=(b-a)/h;
for i=1:N
x=x+h/2;
s=s+f(x)*4;
x=x+h/2;
s=s+f(x)*2;
end
S2=h/6*(s+fa-fb);
end
S=S2

追问

可以给我讲解一下你的思路吗