作业帮 > 综合 > 作业

如何用matlab计算x^3/(exp(x)-1)的积分

来源:学生作业帮 编辑:大师作文网作业帮 分类:综合作业 时间:2024/09/24 07:18:45
如何用matlab计算x^3/(exp(x)-1)的积分
如何用matlab计算x^3/(exp(x)-1)的积分
syms x
>> int(x^3/(exp(x)-1))
ans =
x^3*log(1 - exp(x)) - 6*x*polylog(3,exp(x)) + 6*polylog(4,exp(x)) + 3*x^2*polylog(2,exp(x)) - x^4/4
再问: 如果积分的定义域为[0,1.1433],怎么让最后的值变为一个具体的数字呢?谢谢
再答: 由于我的matlab 是2010b版本 没有maple引擎 不支持polylog函数 所以没法完全转化为数值 syms x y=int(x^3/(exp(x)-1),x,0,1.1433) y = (1494447319737*log(1 - exp(11433/10000)))/1000000000000 + (392140467*polylog(2, exp(11433/10000)))/100000000 - (34299*polylog(3, exp(11433/10000)))/5000 + 6*polylog(4, exp(11433/10000)) - limit(x^3*log(1 - exp(x)) - 6*x*polylog(3, exp(x)) + 6*polylog(4, exp(x)) + 3*x^2*polylog(2, exp(x)) - x^4/4, x = 0, Right) - 17086016206553121/40000000000000000 >> vpa(y) ans = 6.8105760247635834601326874832302 - 1.0*limit(x^3*log(1 - exp(x)) - 6*x*polylog(3, exp(x)) + 6*polylog(4, exp(x)) + 3*x^2*polylog(2, exp(x)) - x^4/4, x = 0, Right) 数值积分的方法可以 > quad(@(x)x.^3./(exp(x)-1),0,1.1433) ans = 0.3166