计算无穷限积分(广义积分)的函数,用Matlab实现,发给大家参考。
1.函数quadFromInf用于计算从负无穷大到一个确定值a的积分。
2.函数quadToInf用于计算从一个确定值a到正无穷大的积分。
--------------------------------------
%name:quadFromInf.m
function Isum = quadFromInf(fun,a,dx0,tol,method)
if nargin < 2 ,a=0 ;end
if nargin < 3 ,dx0=0.5 ;end
if nargin < 4 ,tol = 5e-4 ;end
if nargin < 5 ,method = 1
;end
j=0;dx = dx0;Isum = 0;x2 = a; maxint = 35;
%fprintf('\n j dx x2 I_j Isum\n');
while j
x1 = x2 -
dx;
switch
method
case 1, I = q