io_status read_matrix(double*, int, int, const char*);
void print_matrix(const double*, int, int, int);
void init_matrix(double*, int, int, int);
double f(int, int, int, int, int);
void init_vector(double*, double*, int);