#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define n 1000
#define MU 5
#define ML 5

main(){
  int i,j;
  static double A[n][n],B[n];

FILE *fpr,*fpb;
fpr=fopen("a.dat","w");
fpb=fopen("b.dat","w");

  for(i=0;i<n;i++){
    for(j=0;j<n;j++){
      A[i][j]=0.0;
    }
    B[i]=0.0;
  }

  for(i=0;i<n;i++){
    for( j=0 ; j<n;j++){
      A[i][j]=rand()%10;
    }
  }

  for(i=0;i<n;i++){
    for(j=0;j<n;j++){
      B[i]=B[i]+A[i][j];
    }
  }
   
  for(i=0;i<n;i++){
    fprintf(fpb,"%f ",B[i]);
  }

  for(j=0;j<n;j++){
      for(i=0;i<n;i++){
      fprintf(fpr,"%.1f ",A[i][j]); 
    }
    fprintf(fpr,"\n");
  }
  fclose(fpr);
  fclose(fpb);
}


int min(int former, int latter){
  return (former<latter)? former : latter;
}  
int max(int former, int latter){
  return (former>latter)? former : latter;
}