[MATLAB基础] 求——给定精度求积分的matlab程序

[yfl]

function f=erfen(x)
clear,clc
f=@(x)cos(x);
a=-pi/2;b=pi/2;error=0.0001;
n=100;sum=(f(a)+f(b))*(b-a)/2;
for i=1:n
k=2^i;T=0;
h=(b-a)/k;
X=a+h*[1:k-1];
for j=1:k-1
T=f(X(j))+T;
end
sum(i+1)=h/2*(f(a)+f(b))+T*h;
E=sum(i+1)-sum(i);
if E<error
S=sum(i+1);
break
end
end
disp(S)