<PURE Example Program>

Program pure_example
  Real, Dimension(2,3) :: a
  Forall (i=1:2,j=1:3) a(i,j) = f(i,j)
  Print 1, a
1 Format(6(2X,F0.5))
Contains
  Pure Real Function f(i,j)
    Integer, Intent(In) :: i, j
    f = i + j*0.1
  End Function
End Program

■ Execution Results

  1.10000  2.10000  1.20000  2.20000  1.30000  2.30000