!-----  Ex1_4: 関数の値の計算   -----
      REAL :: f, x, a
      f(x) = x**5 + 6.0*x**4  + 2.0*x**3 + x**2 + 4.0
      READ*, a
      PRINT*, ' f(', a, ')=', f(a)
      END