<Example Program>

Program log_gamma_example
    Real :: x(8) = [ 2.0, 20.0, 200.0, 2000.0, 2e4, 2e5, 2e6, 2e7 ]
    Print 1, x
    Print 1, Log_Gamma(x)
1   Format(999ES9.2)
End Program

Execution Results

 2.00E+00 2.00E+01 2.00E+02 2.00E+03 2.00E+04 2.00E+05 2.00E+06 2.00E+07
 0.00E+00 3.93E+01 8.58E+02 1.32E+04 1.78E+05 2.24E+06 2.70E+07 3.16E+08