N -> Número de muestras
media(x) = sum( x[ i ] ) / N
varianza2 = sum( pow(x[ i ] - media(x), 2) ) / N = sum( pow( x[ i ], 2 ) ) / N - pow( media(x), 2 )
varianza = sqrt( varianza2 )
Código
double calcularDesviacion (int serie[]) { double sumatorio = 0; double sumatorio2 = 0; double promedio = 0; double desviacion = 0; int t; for ( t=0; t < MAXTIRADAS; t++ ) { sumatorio += serie[ t ]; sumatorio2 += pow( serie[ t ], 2 ); } promedio = sumatorio / MAXTIRADAS; double desviacion2 = sumatorio2 / MAXTIRADAS - pow( promedio, 2 ); desviacion = sqrt( desviacion2 ); return desviacion; }
EDITO: He corregido la ecuación. Espero que no tenga más fallos