Juliaで波動関数の確率密度を計算してプロットしたい. 次の関数を評価するために、Julia コードの小さなスニペットを書きました。
ジュリア(不完全)コードは次のとおりです。
set_bigfloat_precision(100)
A = 10
C = 5
m = BigFloat(9.10938356e-31)
ℏ = BigFloat(1.054571800e-34)
t = exp(-(sqrt(C * m) / ℏ))
私が評価する最後の行t
は0.000000000000...
. BigFloat の精度も設定してみました。運が悪い!私は何を間違っていますか?助けていただければ幸いです。