Stohastic / Resistansechain / main

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define _USE_MATH_DEFINES
#define N 20//количество точек в которых находился элетрон
#define C 10
float matozhid(float *x)
{
    float M = 0;
    int k;

for(k=0;k<C;k++)
 M += x[k]*1.;//тут может быть массив вероятностей

return M/C;
}

float dispers(float *x,float m)//dyspersia
{
    float Dy = 0.;
    int k;

for(k=0;k<C;k++)
 Dy += (x[k]-m)*(x[k]-m);

return sqrt(Dy/C);
}

int main()
{
int i,j,k=1;
float R[N],Rab[C],Rm=0;
FILE *fp = fopen("resist.txt","w");


for(i=0;i<C;i++)
Rab[i] = 0;

for(i=0;i<N;i++)
    {
    R[i] = 0.7+(1.*rand())/(RAND_MAX+0.)*0.6;
    printf("%f\n",R[i]);
    }

printf("\n");
fprintf(fp,"%d\n",C);

for(i=0;i<C;i++)
    {
    for(j=0;j<=i*2;j+=2)
        Rab[i] += R[j];
    Rab[i] += R[i+k];
    k++;
    printf("%f\n",Rab[i]);
    fprintf(fp,"%f\n",Rab[i]);
    Rm += Rab[i];
    }

printf("sredneye: %f \n",matozhid(Rab));
printf("dispersia %f \n",dispers(Rab,matozhid(Rab)));


fclose(fp);
//fclose(fp1);
return 0;
}