Fortran language, call IMSL5.0 function library, when a separate call function is a function, can correct solution, but, in a cycle, it is not normal,

As shown, two calculation results,