FORTRAN (1954)



       A=0
       DO 100 I=1,1000,1
       X=REAL(I)/1000
       Y=SQRT(1-X*X)
       A=A+Y/1000
 100 CONTINUE
       WRITE (6,*) 'pi = ',4*A
       END