51 arraystats(double *x, int limit, double *mu, double *sigma, double *rho)  in arraystats()  argument
62 	*mu = sum/(double)n;  in arraystats()
63 	*sigma = sqrt((sumsquare - (double)n*(*mu)*(*mu))/(double)(n-1));  in arraystats()
66 		top += ((double)x[i]- *mu)*((double)x[i-1]- *mu);  in arraystats()
67 		sigma2 += ((double)x[i-1] - *mu)*((double)x[i-1] - *mu);  in arraystats()
93 makedist(double *x, int limit, double mu, double sigma)  in makedist()  argument
107 		input = (x[i]-mu)/sigma;  in makedist()
200 	double mu, sigma, rho;  in main()  local
219 	arraystats(x, limit, &mu, &sigma, &rho);  in main()
222 		limit, mu, sigma, rho);  in main()
225 	table = makedist(x, limit, mu, sigma);  in main()