PROGRAM expressions
  IMPLICIT NONE

  DOUBLE PRECISION :: x, y, z, d

  READ *, x, y, z

  d = x**2 + y**2 + z**2
  PRINT *, "x^2+y^2+z^2 = ", d

  d = x**2/y+z
  PRINT *, "x^2/y+z = ", d

  d = 2.45d0*(x+1.5d0) + 3.1d0*(y+0.4d0) + 5.2d0 - z/2d0
  PRINT  *, "2.45(x+1.5) + 3.1(y+0.4) + 5.2 - z/2 = ", d

  d = (12.8d0*x + 5d0*y) / (11.3d0*y + 4d0*z)
  PRINT *, "(12.8x + 5y) / (11.3 y + 4z) =", d

  d = x**(2.0d0/3.0d0) + y**(2.0d0/3.0d0) + z**(2.0d0/3.0d0)
  PRINT *, "x^(2/3)+y^(2/3)+z^(2/3) = ", d

END PROGRAM expressions
