e = 2.718281828459
func stirling x .
   return sqrt (2 * pi / x) * pow (x / e) x
.
print " X    Stirling"
for i to 20
   d = i / 10
   numfmt 2 4
   write d & "    "
   numfmt 3 4
   print stirling d
.