Exercise 2.2: More Statistics

Given: four three-element arrays, a, b, c, and d, all containing real numbers;
       the formulas for the mean and standard deviation
            mean = (a+b+c+d)/4;
            standard deviation = (((a-mean)**2 + (b-mean)**2 +
                                   (c-mean)**2 + (d-mean)**2) / 3)**0.5

let mean_a := (a[1] + a[2] + a[3])/3.0;
    mean_b := (b[1] + b[2] + b[3])/3.0;
    mean_c := (c[1] + c[2] + c[3])/3.0;
    mean_d := (d[1] + d[2] + d[3])/3.0;

    a_diff_sq_sum := (a[1] - mean_a) * (a[1] - mean_a) +
                     (a[2] - mean_a) * (a[2] - mean_a) +
                     (a[3] - mean_a) * (a[3] - mean_a);
    b_diff_sq_sum := (b[1] - mean_b) * (b[1] - mean_b) +
                     (b[2] - mean_b) * (b[2] - mean_b) +
                     (b[3] - mean_b) * (b[3] - mean_b);
    c_diff_sq_sum := (c[1] - mean_c) * (c[1] - mean_c) +
                     (c[2] - mean_c) * (c[2] - mean_c) +
                     (c[3] - mean_c) * (c[3] - mean_c);
    d_diff_sq_sum := (d[1] - mean_d) * (d[1] - mean_d) +
                     (d[2] - mean_d) * (d[2] - mean_d) +
                     (d[3] - mean_d) * (d[3] - mean_d);

    sigma_a := sqrt( a_diff_sq_sum / 2.0 );
    sigma_b := sqrt( b_diff_sq_sum / 2.0 );
    sigma_c := sqrt( c_diff_sq_sum / 2.0 );
    sigma_d := sqrt( d_diff_sq_sum / 2.0 );

    a_record := [ mean: mean_a; st_dev: sigma_a ];
    b_record := [ mean: mean_b; st_dev: sigma_b ];
    c_record := [ mean: mean_c; st_dev: sigma_c ];
    d_record := [ mean: mean_d; st_dev: sigma_d ];

in  array [ 1: a_record, b_record, c_record, d_record ]
end let