[ruby-list:1366] Re: Rational and Complex

原です。

> なるほど... newton法だと大部簡単なアルゴリズムになるんですね. bignumの
> 割算を実装するには開平法のように実装するしか方法がないので, 平方根も開
> 平法がよいに違いないと信じていたのですが...

bignum の割算て難しいですね。昔 C++ で作ろうとして、挫折、
結局 KNUTH か何かをカンニングしました。わしは割算も出来ない
のかと情けなかった。:-)


> ところで, 冪乗根は, 
> 
> def invpow(n, d)
>   x = d
>   loop do
>     a = ((n - 1) * x**(n-1) + d/x)/n
>     return x if a >= x
>     x = a
>   end
> end
> 
> でいいんですかね??


式は出沢さんのが正しいんですが^^;、n > 2 の時も収束判定がこれでい
いかどうかわかりません。
＃また宿題が増えてしまった。
---===-=-=-=-=-=-=-=-=======--=-=-=-==-=-===-=-=-=-=-=-=--=-==-=--
                            Shin-ichro Hara(Nagaoka Univ.of Tech.)