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

#define MAXFILELEN 80
#define TOLERANCE 10.0e-6
#define MAXITERATIONS 80

// read a matrix and vector, and create a vector


gsl_vector_float CreateVector(int VectorLength);

gsl_vector_float ReadVector();

gsl_matrix_float ReadMatrix();


gsl_vector_float ConjugateGradient(gsl_matrix_float A,gsl_vector_float x,gsl_vector_float b,float Tolerance,int MaxIte);  
void VectorSubtraction(gsl_vector_float *a,gsl_vector_float *b,gsl_vector_float *ans);  
void VectorAddition(gsl_vector_float *a,gsl_vector_float *b,gsl_vector_float *ans);  
void VectorScalar(gsl_vector_float *a,float x,gsl_vector_float *ans);