Расчёт дисперсии

    for (int i=0; i<10000; i++) {
                fscanf(f, "%lf  %lf %lf", &a, &w0[i], &c);
                //printf("%g\n", w0[i]);
                //s1 = s1 + w0[i];
                a1 = a1 + w0[i]*w0[i];
                a2 = a2 + w0[i];
            }   
            //printf("s1 = %lf\n", s1/20000);
            //printf("d1 = %lf\n", a1/20000);
            wf[0] = w0[0];
            p[0] = 0;
            for (int i=1; i<10000; i++) {
                    p[i] = (1 - (p[(i-1)]/((p[(i-1)]+r)))) * p[(i-1)] + q;
                    wf[i] = (1 - (p[(i-1)]/((p[(i-1)]+r)))) * wf[(i-1)] + p[(i-1)]*w0[i]/((p[(i-1)]+r));
                    //printf("%g\n", wf[i]);
                    //s = s + wf[i];    
                    b1 = b1 + wf[i]*wf[i];
                    b2 = b2 + wf[i];
                }
            //printf("s = %lf\n", s/20000);
            //printf("d1 = %lf\n", d/20000);
            d1 = a1/10000 - (((a2*a2)/10000)/10000);
            printf("d1 = %lf\n", d1);
            d = b1/10000 - (((b2*b2)/10000)/10000);
            printf("d = %lf\n", d);
            fprintf(f1, "%lf %lf %lf %lf\n", q, r, s/10000);
            fprintf(f2, "%lf %lf %lf %lf\n", q, r, d1/d);

s1 = s1 + ...

s1+=...

#include <stdio.h>
 
int main () {
 
    double w0[20000];
    double wf[20000];
    double p[20000];
    double a, b, c;
    double q;
    double r;
    double s = 0;
    double s1 = 0;
    double d = 0;
    double d1 = 0;
    double h=1.0;
    double a1=0;
    double a2=0;
    double b1=0;
    double b2=0;
    //printf("Enter q\n");
    //scanf("%lf", &q);
    //printf("Enter r\n");
    //scanf("%lf", &r);
    printf("starting... \n");
    FILE *f;
    FILE *f1;
    FILE *f2;
    //f = fopen ("kalman.txt", "r");    
    if ((f1 = fopen("s.txt", "w")) == NULL)
        printf("file1 was not opened\n");
    if ((f2 = fopen("d.txt", "w")) == NULL)
        printf("file2 was not opened\n");
    for (double q=0.01; q<10; q+=0.01) {
        for (double r=0.0001; r<300; r+=10.0) {
            //f = fopen("data.txt", "r");
             if( (f  = fopen( "data.txt", "r" )) == NULL ) 
                 printf( "file was not opened\n" );
             //else
                 //printf("file was opened successfully");
 
            for (int i=0; i<10000; i++) {
                fscanf(f, "%lf  %lf %lf", &a, &w0[i], &c);
                //printf("%g\n", w0[i]);
                s1 = s1 + w0[i];
                a1 = a1 + w0[i]*w0[i];
                a2 = a2 + w0[i];
            }   
            //printf("s1 = %lf\n", s1/20000);
            //printf("d1 = %lf\n", a1/20000);
            wf[0] = w0[0];
            p[0] = 0;
            for (int i=1; i<10000; i++) {
                    p[i] = (1 - (p[(i-1)]/((p[(i-1)]+r)))) * p[(i-1)] + q;
                    wf[i] = (1 - (p[(i-1)]/((p[(i-1)]+r)))) * wf[(i-1)] + p[(i-1)]*w0[i]/((p[(i-1)]+r));
                    //printf("%g\n", wf[i]);
                    s = s + wf[i];  
                    b1 = b1 + wf[i]*wf[i];
                    b2 = b2 + wf[i];
                }
            //printf("s = %lf\n", s/20000);
            //printf("d1 = %lf\n", d/20000);
            d1 = a1/10000 - (((a2*a2)/10000)/10000);
            printf("d1 = %lf\n", d1);
            d = b1/10000 - (((b2*b2)/10000)/10000);
            printf("d = %lf\n", d);
            fprintf(f1, "%lf %lf %lf %lf\n", q, r, s/10000);
            fprintf(f2, "%lf %lf %lf %lf\n", q, r, d1/d);
             //if ( fclose( f ) )
            fclose(f);
            d=0;
            d1=0;
            a1=0;
            a2=0;
            b1=0;
            b2=0;
            s=0;
            s1=0;
     // {
      //   printf( "The file 'crt_fopen.c' was not closed\n" );
     // }
 
        }
    }  
}
//}

#include <stdio.h>
#include <math.h>
#include <stdio.h>
 
#define  MAX_ARR_SIZE            (20000)
#define  OUTER_LOOP_CNT          (1000)
#define  INNERR_LOOP_CNT         (30)
#define  INPUT_LOOP_CNT          (10000)
#define  CALC_LOOP_CNT           (10000)
 
int main(int argc,char** argv)
{
   FILE*    pIn = fopen("data.txt","rt");
 
   if (!pIn)
   {
      printf("file In wasn't opened.\n");
      return 0;
   }
 
   double   pW1[MAX_ARR_SIZE];
   double   pW2[MAX_ARR_SIZE];
   double   pPP[MAX_ARR_SIZE];
 
   memset(pW1,0,sizeof(double) * MAX_ARR_SIZE);
   memset(pW2,0,sizeof(double) * MAX_ARR_SIZE);
   memset(pPP,0,sizeof(double) * MAX_ARR_SIZE);
 
   double   fSum  = 0.0;
   double   fSum2 = 0.0;
 
   for (int kk = 0; kk < INPUT_LOOP_CNT; ++kk)
   {
      char     pszTemp[MAX_PATH + 1];
 
      fgets(pszTemp,MAX_PATH,pIn);
 
      sscanf(pszTemp,"%Lf",&pW1[kk]);
 
      fSum  += pW1[kk];
      fSum2 += pW1[kk] * pW1[kk];
   }
 
   fclose(pIn);
   pIn = NULL;
 
   FILE*    pOut1 = fopen("s.txt","wt");
 
   if (!pOut1)
   {
      printf("File Out1 wasn't opened.\n");
      return 0;
   }
 
   FILE*    pOut2 = fopen("d.txt","wt");
 
   if (!pOut2)
   {
      printf("File Out2 wasn't opened.\n");
      return 0;
   }
 
   double   fA2 = fSum;
 
   double   qq = 0.01;
 
   for (int ii = 0; ii < OUTER_LOOP_CNT; ++ii, qq += 0.01)
   {
      double   rr = 0.0001;
 
      for (int jj = 0; jj < INNERR_LOOP_CNT; ++jj, rr += 10.0)
      {
         pW2[0] = pW1[0];
         pPP[0] = 0;
 
         double   fSum3 = fSum;
 
         double   fSum21 = 0.0;
         double   fSum22 = 0.0;
 
         for (int ll = 1; ll < CALC_LOOP_CNT; ++ll)
         {
            pPP[ll] = qq + pPP[(ll - 1)] * (1 - (pPP[(ll - 1)] / ((pPP[(ll - 1)] + rr))));
 
            pW2[ll] = (1 - (pPP[(ll - 1)] / ((pPP[(ll - 1)] + rr)))) * pW2[(ll - 1)] + pPP[(ll - 1)] * pW1[ll] / ((pPP[(ll - 1)] + rr));
 
            fSum3 += pW2[ll];
 
            fSum21 += pW2[ll] * pW2[ll];
            fSum22 += pW2[ll];
         }
 
         double   fD1 = 1.E-4 * fSum   - 1.E-8 * fA2    * fA2;
         double   fD2 = 1.E-4 * fSum21 - 1.E-8 * fSum22 * fSum22;
 
         fprintf(pOut1,"%lf %lf %lf %lf\n",qq,rr,1.E-4 * fSum3);
         fprintf(pOut2,"%lf %lf %lf %lf\n",qq,rr,fD1 / fD2);
      }
   }
 
   fclose(pOut1);
   fclose(pOut2);
 
   return 0;
}