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

FUNCTION inva (x)
    k = 0
    WHILE k * (k ^ 2 + 1) / 2 + 2 < x
        k = k + 1
    WEND
    inva = 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)
PRINT
FOR e = 1 TO 20
    PRINT USING "10^##: #########"; e; inva(10 ^ e)
NEXT e
END