#include <math.h>
#include <stdio.h>

#define AXIS 50
#define NUMBER 10000

double gx[NUMBER][AXIS];

int main(int argc, char *argv[]){
int i,k,number,axis,itmp,id;
double gxm[AXIS],dgx[AXIS];
FILE *fin,*fx;

fin=fopen(argv[1],"r");
fx=fopen(argv[2],"w");

(void)printf("number of configs ?\n");
(void)scanf("%d",&number);
if(number>NUMBER){printf("oops - redimension arrays!\n");exit(1);}
(void)printf("lattice length ?\n");
(void)scanf("%d",&axis);
if(axis>AXIS){printf("oops - redimension arrays!\n");exit(1);}

for(k=0;k<number;k++)
for(i=0;i<axis;i++)
(void)fscanf(fin,"%d%lg",&id,&gx[k][i]);

for(i=0;i<axis;i++){
gxm[i]=0.0;
for(k=0;k<number;k++){
gxm[i]+=gx[k][i];
}
gxm[i]/=number;
}

for(i=0;i<axis;i++){
dgx[i]=0.0;
for(k=0;k<number;k++){
dgx[i]+=(gx[k][i]-gxm[i])*(gx[k][i]-gxm[i]);
}
dgx[i]=sqrt(dgx[i])/number;
}

for(i=0;i<axis;i++){
(void)fprintf(fx,"%d %lg %lg\n",i,gxm[i],dgx[i]);
}

return(0);
}