FUNCTION a(n)
    LET n = n + 2
    LET a = n*(n^2 + 1)/2
END FUNCTION

FUNCTION inv_a(x)
    LET k = 0
    DO WHILE k*(k^2+1)/2+2 < x
       LET k = k + 1
    LOOP
    LET inv_a = k
END FUNCTION

PRINT "The first 20 magic constants are: ";
FOR n = 1 TO 20
    PRINT a(n);" ";
NEXT n
PRINT
PRINT "The 1,000th magic constant is "; a(1000)

FOR e = 1 TO 20
    PRINT USING "10^##": e;
    PRINT USING": #########": inv_a(10^e)
NEXT e
END