program-id.y.data division.working-storage section. 1 m pic 99. 1 n pic 9(4). 1 s pic x(5). procedure division.a.accept s move s(1:2)to m move s(3:)to n compute n=n*(1+n)*(1+m*3- n)/6 display n(4- function log10(n):)if n=819 display 53 else if n=968 display"50 50"else go a.
Note that non-ascii characters in the above source code will be escaped (such as \x9f).