Fortran program for area of circle