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

#define N_kezdeti 50
#define m  5000
#define j  50

FILE *fp;

float k_N_medium[j+1];

main()
 {
   int k, r, i, ii, ij, N;
	for(i=1; i<j+1; i++)
	 {
	 N=i*N_kezdeti;
	 for(ii=1; ii<m+1	; ii++)
		{ k=0;
		for(ij=1; ij<N+1; ij++)
			{
			k=k+((int)((float)(rand())/(RAND_MAX+1.0)*2)*2-1);
			}
		 k_N_medium[i]=(float)(k_N_medium[i]*(ii-1)+k*k)/(float)(ii);
		}
	 }
  fp=fopen("result.dat","a");
  for(ii=1; ii<j+1; ii++)
	 fprintf(fp,"%d    %f\n", ii*N_kezdeti, k_N_medium[ii]);
  fclose(fp);
	 
}