<!-- Original:  Phil Pavlik (ppavlik@nmu.edu) -->
<!-- Contributor URL:  http://chemistry.nmu.edu/ppavlik  -->

<!-- This script and many more are available free online at -->
<!-- The JavaScript Source!! http://javascript.internet.com -->

<!-- Begin
function calculate() {
var mean, sdev, N, sign, sum, str;
var y = Array(), x = Array(), freq = Array();
mean = document.statForm.m.value * 1;
sdev = document.statForm.s.value * 1;
N = document.statForm.n.value * 1;
sdevsq = sdev * sdev;
for(var i = 0; i <= N; i++) {
yr = 1;
yxr = 0;
while(yr > yxr) {
xr = 500 * sdev * Math.random();
yr = 100 * Math.random();
yxr = 100 * Math.exp(-(xr) * (xr) / (2 * sdevsq));
}
s = 1;
if(Math.random() < .5)
s = -1;
x[i] = mean + s * xr;
}
str = "";
for(var i = 1; i <= N; i++) {
str = str + x[i] + "\n";
}
sum = 0;
for(var i = 1; i <= N; i++) {
sum = sum + x[i];
}
mean = sum / N;
sum = 0;
for(var i = 1; i <= N; i++) {
sum = sum + (x[i] - mean) * (x[i] - mean)
}
sdev = Math.sqrt((sum) / (N - 1));
document.statForm.vals.value = str;	
document.statForm.mm.value = mean;
document.statForm.ss.value = sdev;
}
//  End -->