|
|
- 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.);
复制代码
上面这段代码是我写的。定义了函数 InG(x,a,b,h). 其中 x 取自方程 Γ(s) = x, (a,b) 使得 (Γ(a)-x)(Γ(b)-x) < 0 且 a,b 之间 Γ(s) = x 只有一个解 s. 最后 h 是精度控制参数. 使用这个函数前要将其定义告诉软件。
下载安装 Pari/GP 后,软件自带一个简单的教程。这个法国的自由软件主要是为搞数论的人设计的。其语言语法有点怪,但可处理不少 Mathematica 搞不了的问题。有问题可以提出讨论。数学研发站大约有比较懂这个软件的朋友。 |
|