extern double a, b, h;
extern int M, N, m, j, k;
#undef PI
extern double PI;
extern VEC   *u_vec;