C言語での平均値、分散、標準偏差の計算例

#include<stdio.h>
#include "math.h" 
 double sum = 0;//  
 double array[4] = {1.2,2.1,3.1,4.1};
 int length = 0;//    
 double average = 0;//    
 double var = 0; //   
 double standard = 0; //     
 int main()
 {
  for (int i = 0; i <= 3;i++)
  {
   sum += array[i];//  
  }
 
  length = sizeof(array) / sizeof(array[0]);//     
  average = sum / length;//    
 
  for (int j = 0; j <= 3;j++)
  {
   var += pow(array[j]-average,2)/length;//   
  } 
  standard = pow(var,0.5);//    
 
  printf("sum = %g
",sum);//         printf("%f
",sum);               6 
  printf("length of array[4] = %d
",length);
  printf("average = %g
",average);
  printf("var = %g
",var);
  printf("standard = %g
",standard);
  while (1);
  return 0;
 }

#include<math.h>
#include<stdio.h>
double fun(double x[10])
{
 double p = 0.0,f=0.0,g=0.0;
 for (int i = 0; i < 10; i++)
 {
  p += x[i];
 }
 p = p / 10;
 printf("%f
", p);
 for (int j = 0; j < 10; j++)
 {
  g = x[j] - p;
  f += g*g;
  printf("%f
", f);
 }
 f = f / 10;
 f = sqrt(f);
 return f;
}
void main()
{
 FILE *wf;
 double s,x[10]={95.0,89.0,76.0,65.0,88.0,72.0,85.0,81.0,90.0,56.0};
 int i;
 printf("
The original data is:
");
 for(i=0;i<10;i++)
  printf("%6.1f ",x[i]);
 printf("

 ");
 s=fun(x);
 printf("s=%f

 ",s);
/******************************/
 wf=fopen("out.dat","w");
 fprintf(wf,"%f",s);
 fclose(wf);
/*****************************/
}