最近,我在julia有一个与反拉普拉斯有关的问题。
这是我在茱莉亚写的代码。
using SpecialFunctions
using InverseLaplace
value = Talbot(u -> gamma_inc(1,u,2)[2], 80)
value(2)
此代码会产生一些错误消息,例如
MethodError: no method matching gamma_inc(::Float64, ::BigFloat, ::Int64, ::String)
实际上,我已经可以基于mpmath在python中计算出这种功能了。
import mpmath as mp
laplace = lambda u: mp.gammainc(1,u)
result = mp.invertlaplace(laplace, t = 2, method= 'talbot',degree = 80)
print(result)
结果是
5.6730679553874e-43
我想在茱莉亚上解决这个问题,有人可以给我任何建议吗?
p.s。对不起,我的英语不好。希望你能理解。.