clear;clc;close all
k=1;
a=-10; %积分下限
b=10; %积分上限
global x
for x=-10:0.5:10;
while (1)
y1=quad(@myfun1,a,-2)+quad('myfun1',2,b);
y2=quad(@myfun1,2*a,-2)+quad('myfun1',2,2*b);
if abs((y1-y2)/(y1+eps))<1e-5
y(k)=y2;
k=k+1;
break
else
a=2*a;
b=2*b;
end
end
end
x=-10:0.5:10;
plot(x,y)
xlabel('x')
ylabel('y')
grid on
