home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD1.mdf / basic / compiler / ubasic / malm / rennet.ub < prev    next >
Encoding:
Text File  |  1990-08-22  |  669 b   |  18 lines
  1.    10   *Rennet(&Cf(),Leng,&A,&B,&C)
  2.    20   ' Periodic continued fraction to quadratic.
  3.    30   ' CF() is an input parameter - passed by reference for efficiency.
  4.    40   ' Modeled on the Pascal version.
  5.    45   ' 7 May 1990
  6.    50   local I,H1=0,H2=1,K1=1,K2=0,Te
  7.    60   if Leng<0 then A=0:B=0:C=0:return endif
  8.    70   if Leng=0 then A=1:B=0:C=-Cf(0):return endif
  9.    80   for I=1 to Leng
  10.    90   Te=H1+Cf(I)*H2:H1=H2:H2=Te
  11.   100   Te=K1+Cf(I)*K2:K1=K2:K2=Te
  12.   110   next I
  13.   120   A=H1:B=(H2-K1)-2*H1*Cf(0)
  14.   140   C=Cf(0)*(Cf(0)*H1-(H2-K1))-K2
  15.   150   Te=gcd(B,A):Te=gcd(C,Te)
  16.   160   A=A\Te:B=B\Te:C=C\Te
  17.   170   return ' End of subroutine Rennet.
  18.