void gauss_seidel ( int n, double r[], double u[], double &dif_l1 );
void poisson_1d ( int n, double a, double b, double ua, double ub,
  double force ( double x ), double exact ( double x ), int &it_num, 
  double u[] );