计算16个数的标准差（在标准差是整数的情况下）

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

double ret=0;

my_pow(double z)
{
    return (z - ret) * (z - ret);
}

int main()
{
    double a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p = 0;
    scanf("%lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf", &a, &b, &c, &d, &e, &f, &g, &h, &i, &j, &k, &l, &m, &n, &o, &p);
    ret = (a + b + c + d + e + f + g + h + i + j + k + l + m + n + o + p) / 16;
    double zong1 = 0;
    zong1 = (my_pow(a) + my_pow(b) + my_pow(c) + my_pow(d) + my_pow(e) + my_pow(f) + my_pow(g) + my_pow(h) + my_pow(i) + my_pow(j) + my_pow(k) + my_pow(l) + my_pow(m) + my_pow(n) + my_pow(o) + my_pow(p));
    double zong2=0;
    zong2=sqrt(zong1*0.0625);
    printf("%lf",zong2);
    return 0;
}

因为double的精度不够，方法是这么个方法