Double Precision Function Fortran