program jastrow
  implicit none
  print *, 'The total Jastrow factor'
  print *, dexp(factor_ee + factor_en + factor_een)

end program