﻿功能:无限区间到无限区间积分

格式:
[y,error,np]=IntegralInfiniteD(f,n,error0,xn)
[y,error,np]=IntegralInfiniteD(f,n,error0)
[y,error,np]=IntegralInfiniteD(f,n)

f:符号变量存储的积分表达式
n:复化积分数等分最大深度值(等分数=64*2^n),默认为8
error0:表示控制的相对误差,默认1E-12
xn:积分变量名称,默认为x

y:返回的积分值
error:返回的积分相对误差
np:积分时的复化等分数

原理:采用两点高斯积分求解

例子:

//对x^2*exp(-x^2)函数x在-inf到inf的定积分
f = "x^2*exp(-x^2)";
[y,er,np]=IntegralInfiniteD(f)//回车得到如下结果
y =
[ 0.88622692545276 ]
er =
[ 1.5043521983E-13 ]
np =
[ 1024.00000000000 ]