program-id.y.data division.working-storage section.
1 i pic 999.
1 k pic 999.
procedure division.a.compute k=function mod(1+i,3)compute k=i+k+i*0^k
display k(3- function log10(k):)add 1 to i
if i<401 go a.