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[] );