■無理数(その9)

【3】√3のディオファントス近似とペル方程式

an^2−3bn^2=−1

は解をもちませんが,

an^2−3bn^2=1

が成り立つ最小解は(a,b)=(2,1)であることから,√3に収束する分数列を構成すると

  (2+√3)^n=an+bn√3

  (2−√3)^n=an−bn√3

より

  an+1+√3bn+1=(2+√3)(an+√3bn)

          =(2an+3bn)+√3(an+2bn)

  an+1=2an+3bn

  bn+1=an+2bn

  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

となります.

===================================