void poisson_2d_exact ( int n, double x[], double y[], double u[], double ux[], 
  double uy[], double uxx[], double uxy[], double uyy[] );