#include <stdio.h>
#include <stdlib.h>
#include <gsl/gsl_blas.h>

#define MAXLEN 80
#define TOL 10.0E-6
#define MAXITE 100

int cgaf(gsl_matrix_float *A,gsl_vector_float *x,gsl_vector_float *b,float Tol,int MaxIte);
void vecsub(gsl_vector_float *a,gsl_vector_float *b,gsl_vector_float *ans);
void vecsca(gsl_vector_float *a,double x,gsl_vector_float *ans);
void vecread(gsl_vector_float *a,gsl_vector_float *b,gsl_vector_float *ans);

int kain(gsl_matrix_float **A,gsl_vector_float **b);
int kaiout(gsl_vector_float *x);