我想在Matlab中用可变精度算法评估一个定积分。可以使用符号数学工具箱以这种方式完成:
syms x
f = (x.^2000).*((1-x).^4000)
vpa(int(f,0,1))
这给出了变量精度算术积分的答案。
但我喜欢在没有符号数学工具箱的情况下评估积分。我可以使用命令'integral'来计算积分,但由于积分是以固定精度计算的,所以它返回零,即下面代码的输出为零。
f = @(x) (x.^2000).*((1-x).^4000)
integral(f,0,1)
如何在不使用符号数学工具箱的情况下将vpa与数值积分相结合?