Super Magic Box

Put in a value for D to compute its Super Magic Box:

n
A n = a n 1 C n 1 A n 1
C n = ( D A n 2 ) / C n 1
a n = ( A n + a 0 ) / C n
p n = a n p n 1 + p n 2
q n = a n q n 1 + q n 2
p n 2 D q n 2

The fundamental unit in [ D ] is a + b D , with norm N

A factorization of D is given by g c d ( p + C , D ) = d D = d * ( p + C )