double media (double *p, int elementos){

    int i;
    double sumatorio=0.0, media=0.0;

    for (i = 0; i < elementos; i++){
        sumatorio = sumatorio + *p;
        p++;
    }

    media = sumatorio / elementos;

    return(media);
}


double varianza (double *p, int elementos){

    int i;
    double sumatorio = 0.0, var = 0.0, med = 0.0;

    med = media (p, elementos);

    for (i = 0; i < elementos; i++){
        sumatorio = sumatorio + ((*p - med)*(*p - med));
        p++;
    }

    var = sumatorio / (elementos -1);
    return (var);

}