on rounded$ %% numeric, rather than symbolic, computation precision 25$ %% floating-point decimal precision golden_ratio := (1.0 + SQRT(5.0))/2.0$ lo := hi := n := 1$ WRITE n, " ", hi$ WHILE hi < 2^63 do << n := n + 1; ratio := hi / lo; WRITE n, " ", hi, " ", ratio, " ", (ratio - golden_ratio); hi := lo + hi; lo := hi - lo; >>;