Pell's equation: Difference between revisions

Content added Content deleted
m (added whitespace.)
m (→‎{{header|REXX}}: used a better form of over/under statement placement.)
Line 987: Line 987:
z= 1; r= x + x
z= 1; r= x + x
do until ( (e2 + x*f2)**2 - n*f2*f2) == 1
do until ( (e2 + x*f2)**2 - n*f2*f2) == 1
y= r*z - y
y= r*z - y; z= floor( (n - y*y) / z)
z= floor( (n - y*y) / z)
r= floor( (x + y ) / z)
r= floor( (x + y ) / z)
parse value e2 r*e2 + e1 with e1 e2
parse value e2 r*e2 + e1 with e1 e2
parse value f2 r*f2 + f1 with f1 f2
parse value f2 r*f2 + f1 with f1 f2
end /*until*/
end /*until*/
return e2 + x * f2 f2</lang>
return e2 + x * f2 f2</lang>