﻿功能:有限区间到正无穷区间积分

格式:
[y,error,np]=IntegralInfinite2(f,a,n,error0,xn)
[y,error,np]=IntegralInfinite2(f,a,n,error0)
[y,error,np]=IntegralInfinite2(f,a,n)
[y,error,np]=IntegralInfinite2(f,a)
f:符号变量存储的积分表达式
a:积分下限,可以为数值,可以为矩阵变量
n:复化积分数等分最大深度值(等分数=64*2^n),默认为8
error0:表示控制的相对误差,默认1E-12
xn:积分变量名称,默认为x

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

原理:采用复化后的2点高斯积分求解

注意:本函数与IntegralInfinite()函数不同的是,其变换方式不一样

例子:

//对sin(x)^3*cos(x)/x在[0,inf]上进行积分,这个积分使用函数{IntegralInfinite<矩阵运算\IntegralInfinite>}会计算不出
f = "sin(x)^3*cos(x)/x";
[y,er,np]=IntegralInfinite2(f,0)//回车得到如下结果,真实积分结果=π/16=0.19634954084936
y =
[ 0.19075279277901 ]
er =
[ 0.00496619343024 ]
np =
[ 128.000000000000 ]

 