function h(z)
  integer, parameter :: DP=kind(0.0D0)
  real(kind=DP):: z,h
  real(kind=DP),parameter :: twopi=6.28318530717959_DP
  h=sin(twopi*z)
end function h

function F(x,y)
  integer, parameter :: DP=kind(0.0D0)
  real(kind=DP):: F,x,y
  F=2*exp(x+y)+3*x*x+6*y*y+sin((x+y)/2)+cos((x+y)/2)
end function F

function up(x,y)
  integer, parameter :: DP=kind(0.0D0)
  real(kind=DP):: up,x,y
  up=exp(x+y)+x*x*x+2*y*y*y+sin((x+y)/2)-cos((x+y)/2)
end function up