■ペル数(その3)
同様にして,√3の最良近似では
(2+√3)^n=an+bn√3
(2−√3)^n=an−bn√3
より
an+1=4an−an-1,bn+1=4bn−bn-1
α,βを2次方程式x^2−4x+1=0の根2±√3として,初期値をa1=1,a2=2,a3=7,b1=0,b2=1,b3=4とすると
an/bn→ √3
となります.
近似分数列{an/bn}で非常によく近似できる実数αについて
|α−an/bn|<1/bn^2
が成立するならば,αは無理数です(ディリクレの定理).
===================================