計算機科学のブログ

大きい数を含む計算とWolfram AlphaとMathematica、Wolfram Engine

Wolfram AlphaではなくWolfram Engineで、ローカルではどうなるのか気になったから試してみた。

コード(Wolfram Language)

Integrate[(x^(5 10^12) - 1) / (Sqrt[x] Log[x]), {x, 0, 1}]
PolynomialGCD::lrgexp: Exponent is out of bounds for function PolynomialGCD.



PolynomialGCD::lrgexp: Exponent is out of bounds for function PolynomialGCD.



PolynomialGCD::lrgexp: Exponent is out of bounds for function PolynomialGCD.



Further output of `1` will be suppressed during this calculation.: Further output of PolynomialGCD::lrgexp will be suppressed during this calculation.
Output
f[n_] := Integrate[(x^n-1) / (Sqrt[x] Log[x]), {x, 0, 1}]
f[n]
Output
f[5 10^12]
PolynomialGCD::lrgexp: Exponent is out of bounds for function PolynomialGCD.



PolynomialGCD::lrgexp: Exponent is out of bounds for function PolynomialGCD.



PolynomialGCD::lrgexp: Exponent is out of bounds for function PolynomialGCD.



Further output of `1` will be suppressed during this calculation.: Further output of PolynomialGCD::lrgexp will be suppressed during this calculation.
Output
Simplify[f[n], Element[n, PositiveIntegers]] // TraditionalForm
Output

ということで、多少時間はかかるものの、計算できた。るめなるさん(@Re_menal2)log(10000000000001)、log(10兆1)と一致。また、pepperさん(@pepper_aobuta)の文字を使う場合について、その条件も表示してくれた。さらに、nを正の整数と指定すれば(Element[n, PositiveIntegers])、Log(2n + 1)という結果も表示してくれた。