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

int main(int argc, char **argv)
{
    FILE *in;
    int n = 0, s;
    double x, X, sum = 0, sumsq = 0;
    if (argc != 2) {
	printf("Usage: %s <filename>\n", argv[0]);
	return 1;
    }
    in = fopen(argv[1], "r");
    if (in == NULL) {
	printf("Unable to open %s\n", argv[1]);
	return 1;
    }
    s = fscanf(in, "%lf", &x);
    while (s == 1) {
	n++;
	sum += x;
	sumsq += (x * x);
	// printf("%f - %d\n", x, s);
	s = fscanf(in, "%lf", &x);
    }
    printf("mu_mle = %lf, sigma_mle = %lf\n", X = sum/n, sqrt(sumsq/n - (X)*(X)));
}