int solve_matrix_equation(int nrows, int ncols, 
		          double *matrix_a, double *vector_b, 
			  double *vector_x, double *vector_s);