<Example Program>

Program erf_example
    Real :: x(6) = [ 0.5,1.0,2.0,3.0,4.0,5.0 ]
    Print 1, x
    Print 1, Erf(x)
1   Format(999F9.6)
End Program

Execution Results

 0.500000 1.000000 2.000000 3.000000 4.000000 5.000000
 0.520500 0.842701 0.995322 0.999978 1.000000 1.000000