program main implicit none integer,parameter :: a=10000 integer :: c,b,g,d,e integer,dimension(8401) :: f do b=0,8400 f(b)=a/5 end do e=0 do c=8400,14,-14 do b=c-1,1,-1 g=2*b-1 d=d*b+f(b)*a f(b)=mod(d,g) d=d/g end do if(c==8400) then write(6,'(f5.3)',advance='no') real(e+d/a)/1000 else write(6,'(i4.4)',advance='no') e+d/a end if e=mod(d,a) end do end program main