咖玛函数Γ(x) = 1/3 ，求x的解

永远

杯具了，Wolfram Mathematica 9尽然求不出，还是我编程错了，求指导。如果本软件求不出果断向elim老师学他的软件

永远

走过，路过的高手帮忙看看吧……

elim

Mathematica 目前还有很多事情干不了, 干不完全. 它在不断改进中。它改进的速度没有它体积膨胀的速度快：）
由

可见Mathematica 可以算反Gamma函数值，只是反函数是多值的情况，不知道如何控制返回的是哪个值。

1. (14:57) gp > \p 100
   
2.    realprecision = 105 significant digits (100 digits displayed)
   
3. (14:58) gp > InG(x,a,b,h)=my([s,t]=[a,b]);while((abs(gamma(s)-x)>h)&&(w=((s+t)/2.))&&(v=(gamma(s)-x)*(gamma(w)-x)),if(v> 0,s=w, v< 0, t=w, s=t=w));return((s+t)/2.);
   
4. (14:58) gp > InG(1./3,-3.3,-3.4,10^(-101))
   
5. %2 = -3.390724059691994004341696638396287543328933107815136438770280261925836049189099333821049034411524956
   
6. (14:58) gp > gamma(%)
   
7. %3 = 0.3333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333
   
8. (14:58) gp >

复制代码

elim

我的程序也能得到 Mathematica 的解。精度也可以任意提高.

永远

elim 发表于 2018-7-25 05:06
我的程序也能得到 Mathematica 的解。精度也可以任意提高.

果断的选择学习老师楼下这个软件！！…

elim

数学软件各有优缺点．

永远

elim 发表于 2018-7-25 07:52
数学软件各有优缺点．

可有相关课件，不然我都不知道那是那……

Ysu2008

PARI/GP 强项在哪个方面？

jzkyllcjl

我已经 称赞了elim的 算法，但是 我仍然坚持 这个无尽小数是一个数列性质的变数，它的极限才满足Г（α）=1/3。对于这个 无尽小数，只能写出 全能近似等式 Г（-3.3907240……）~1/3。
在elim提出这个解法之后，笔者查看了菲赫金哥尔茨《微积分学教程》二卷三分册716-717页渐进级数一节，根据这一节对无穷渐进无穷级数的等价表达符号以及有穷级数的近似表达式，笔者对无尽小数的上述改革意见是应当提出的。

Ysu2008

jzkyllcjl 发表于 2018-7-25 11:11
我已经 称赞了elim的 算法，但是 我仍然坚持 这个无尽小数是一个数列性质的变数，它的极限才满足Г（α）=1 ...

你的理论不能说服任何人，哪怕一个支持者也没有。
这么好几十年下来，最后还是光杆司令。
你就没想过什么环节出了问题？