编码如下，平均值结果正确，然而标准差是错误的。

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

// declare functions
float mean(float array[], int n);
float stdev(float array[], int n);

// define function
float mean(float array[], int n)
{
 float arraysum = 0.;
 int i;
 for (i=0; i<n; i++)
{
 arraysum = arraysum + array[i];
}
 return arraysum/n;
}

//define function
float stdev(float array[], int n)
{
 float arraysum = 0.;
 float avgl = mean(array,n);
 int i;
 for (i=0; i<n; i++)
   {
      arraysum = (arraysum + array[i]-avg1)*(arraysum + array[i]-avg1);
   }
 return sqrt(arraysum/n);
}
  
int main()
{
 float input[10];
 float avg;
 float sd;

 printf("enter ten floats\n");
 scanf("%f %f %f %f %f %f %f %f %f %f",&input[0],&input[1],&input[2],&input[3],&input[4],&input[5],&input[6],&input[7],&input[8],&input[9]);
 printf("%f %f %f %f %f %f %f %f %f %f\n",&input[0],input[1],input[2],input[3],input[4],input[5],input[6],input[7],input[8],input[9]);  
 // call function 1
 avg = mean(input,10);
 printf("array mean = %f\n",avg);
 // call function 2
 sd = stdev(input,10);
 printf("array stdev = %f\n",sd);
 return 0;
}

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

float mean( const float array[], size_t n )
{
    double sum = 0.;
    for( size_t i=0; i!=n; ++i )
        sum += array[i];
    return sum/n;
}

float stdev( const float array[], size_t n )
{
    double sum = 0.;
    double avg = mean(array,n);
    for( size_t i=0; i!=n; ++i )
        sum += (array[i]-avg)*(array[i]-avg);
    return sqrt(sum/n);
}

int main( void )
{
    float input[10];
    printf("enter ten doubles\n");
    scanf("%f %f %f %f %f %f %f %f %f %f",&input[0],&input[1],&input[2],&input[3],&input[4],&input[5],&input[6],&input[7],&input[8],&input[9]);
    printf("%f %f %f %f %f %f %f %f %f %f\n",input[0],input[1],input[2],input[3],input[4],input[5],input[6],input[7],input[8],input[9]); 

    float avg = mean(input,10);
    printf( "array mean = %f\n", avg );

    float sd = stdev(input,10);
    printf( "array stdev = %f\n", sd );

    return 0;
}

怎么我把
arraysum = (arraysum + array[i]-avg1)*(arraysum + array[i]-avg1);
改成
arraysum += (array[i]-avg1)*(+ array[i]-avg1);
就可以输出标准差的值，而原本的是无穷呢，他们不是一样的吗？
 

这都是啥啊