double precision function fcuad(x)
  implicit none
  double precision, intent(in) :: x
  fcuad = x*x
end function fcuad