program-id.y.data division.working-storage section. 1 i pic 9999. 1 j pic 9. 1 k pic 9999. 1 n pic 9(7). procedure division.a.add 1 to j compute k=k/j compute n=n+function mod(k,1+j)*10^j if k>0 go a.display n(7- function log10(n):)if i<1000 add 1 to i move i to k move 0 to n j go a.
Note that non-ascii characters in the above source code will be escaped (such as \x9f).