program tmp

  ! Antoine の式を利用. H2O の係数は化学工学便覧参照.
  real(8), parameter  :: A = 23.1964d0
  real(8), parameter  :: B = 3816.44d0
  real(8), parameter  :: C = -46.13d0
  real(8) :: temp
  integer :: i

  do i = 1, 10
    temp = 273.15d0 + 10.0 * (i - 1)

    write(*,*)  temp, dexp( A - B / (C + temp)  )

!    write(*,*)  temp, ( B / (C + temp)** 2.0d0 ) *  8.31 * (temp**2.0d0)
  end do

end program tmp