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 4 2
   write d & "    "
   numfmt 4 3
   print stirling d
.