#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define N 1000

double A[N*N];
double x[N];

main(){
  FILE *fpa,*fpb;
  int i,j,k,c;
  double dum,pivot,aik,akj,max;
  static long int n=N,inc=1,info,piv[N];
  /* input */

  fpa=fopen("a.dat","r");
  fpb=fopen("b.dat","r");

  for(i=0;i<n*n;i++)
    fscanf(fpa,"%lf",&A[i]);

  fclose(fpa);

  for(i=0;i<n;i++){
      fscanf(fpb,"%lf",&x[i]);
    }
 
  fclose(fpb);


/* check :pivot */

  printf("N = %d\n",N);
  dgesv_(&n,&inc,A,&n,piv,x,&n,&info);
  for(i=0; i<N; ++i) printf("%lf\n", x[i]);
  return(0);

}