|
|
本帖最后由 elim 于 2018-8-1 13:49 编辑
- g(d,b,e,s)=my(v=1.);my(h=log(b));while((w=v-(b^v+s*v^e-d)/(h*b^v+s*e*v^(e-1.)))&&(abs(w-v) > 1./10^101),v=w);return(v);
- (14:43) gp > g(3,Pi,Pi,1)
- %29 = 0.80107847391239416228251685689215059553685369472464879054225588916736944452088584841...
- (14:44) gp > g(250,Pi,Pi,1)
- %30 = 4.36337526396209710142356027640407705627054882291557618687506721328919018960974808208...
- (14:44) gp > g(0,sqrt(Pi),sqrt(Pi),-1)
- %31 = 1.77245385090551602729816748334114518279754945612238712821380778985291128459103218137...
复制代码 |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有帐号?注册
x
|