void poisson_2d ( int nx, int ny, void rhs ( int nx, int ny, double **f ), 
  double u_exact ( double x, double y ) );
void jacobi ( int nx, int ny, double dx, double dy, double **f, 
  double **u, double **unew );
double rms_norm ( int m, int n, double **a );