[MATLAB项目实例源码]求解无穷区间定积分问题 源程序代码.zip

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 msgbox('MATLAB编程答疑，请加QQ: 1530497909','MATLAB答疑','help') web http://url.cn/NSFcAs -browser