fortran: first try to numerical calculate pi