#include <stdlib.h>
#include <stdio.h>
#include "matrix_io.h"

int main(int argc, char *argv[]) {

  if (argc != 4) {
    fprintf(stderr,"Usage: %s A b x\n", argv[0]);
    return EXIT_FAILURE;
  }

char *A = argv[1];
char *b = argv[2];
char *x = argv[3];

  /* Insert your code here */

// Read input
vector_t *b_ = read_vector(b);
matrix_t *A_ = read_matrix(A);

/****************Koden virker hertil**************/
// Call call_dgels
double info, work = call_dgels(A_,b_);

// write output
printf("%f",work);


unsigned long n = b_->n;

struct vector *work_;
work_->n = n;
work_->v = &work_[0];

write_vector(x, work_);

  return EXIT_SUCCESS;
}