20/02/01 11:30:55 6JPjRL8V.net
Maxima:
Brentせんせーの真似
expm1(x):=if abs(x) < 1e-8 then x*x/2+x else block([y:expm1(x/2)], y*y+2*y);
一方、英語版Wikipediaに拠ると
URLﾘﾝｸ(en.wikipedia.org)
expm1(x):=2*tanh(x/2)/(1-tanh(x/2));
tanh()があって、expm1()が無い場合はこうすればいいのか…。
20/02/01 11:30:55 6JPjRL8V.net
Maxima:
Brentせんせーの真似
expm1(x):=if abs(x) < 1e-8 then x*x/2+x else block([y:expm1(x/2)], y*y+2*y);
一方、英語版Wikipediaに拠ると
URLﾘﾝｸ(en.wikipedia.org)
expm1(x):=2*tanh(x/2)/(1-tanh(x/2));
tanh()があって、expm1()が無い場合はこうすればいいのか…。