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