function factor(){

var sn = "tfac";
var w=window.open("tfactors.gif",sn,'scrollbars=no,width=250,height=275');
var d = w.document;
d.close();
w.focus();
}

function printit(){

var i = 1;
var wd = 0;
var ht = 0;
var sumW = 0;
var sumH = 0;
var sumHW = 0;
var sumWW = 0;
var sumHH = 0;
var sumHWHW = 0;
var sumP = 0;
var sumPP = 0
var avgH = 0;
var avgW = 0;
var avgHW = 0;
var avgP = 0;
var sigmaH = 0;
var sigmaW = 0;
var sigmaHW = 0;
var sigmaP = 0;
var ablank = "______";
var dummy1 = "";
var dummy2 = "";
var dummy3 = "";
var dummy4 = "";
var zz = 1;
var zz2 = 1;

var XX = parent.top.BB.myX;
var YY = parent.top.BB.myY;


var howmany = parent.top.BB.whichtry;
var header = '<td valign = bottom align = center><b>';
var tailer = '</b></td>';
var isay = '<table border = 0 style = "font-size:12 pt"><tr>' + header + "#" + tailer + header + 'Width(px):' + tailer +header + 'Height(px):' + tailer + header + 'Perim(px):'+ tailer + header + 'Area(px<sup>2</sup>):'+ tailer + '</tr>';

var myName = prompt("Please enter your name","Joe Miner");
var today = new Date();
today = today.toLocaleString(); 
sn = "A";
w=window.open("",sn,"resizable,menubar,scrollbars")
d = w.document

if (document.layers){
d.write('<BR><center><font size = +1><b>STATISTICAL ANALYSIS OF EXPERIMENTAL DATA</b></font>');
}else{
d.write('<BR><center><b>STATISTICAL ANALYSIS OF EXPERIMENTAL DATA</b>');
}
d.write('<BR><BR><b>' + today + '&nbsp;&nbsp;_______________________________&nbsp;&nbsp;' + myName + '</b><BR><BR>');



	for (i==1;i <= howmany;i++){
wd = parent.top.BB.cuts[i].ww;
ht = parent.top.BB.cuts[i].hh; 
sumH += ht;
sumW += wd;
sumHH += ht*ht;
sumWW += wd*wd;
	sumHW += ht*wd;
	sumHWHW += ht*ht*wd*wd;
sumP += 2*wd + 2*ht;
sumPP += 4*(wd + ht)*(wd + ht);  
	isay += '<tr>' + header + i + tailer + header + wd + tailer;
	isay += header +  ht + tailer;
		isay += header + (2*ht + 2*wd) + tailer; 
		isay += header + ht*wd + tailer + '</tr>';
	}  
i = 1;
isay += '<tr><td>&nbsp;</td><td></td><td></td><td></td></tr>';
isay += '<tr>' + header + "Sum:" + tailer + header  + sumW + tailer;
isay += header + sumH + tailer;
isay += header + sumP + tailer;
isay += header + sumHW + tailer + '<tr>';		
avgH = sumH/howmany;
avgW = sumW/howmany;
avgHW = sumHW/howmany;
avgP = sumP/howmany;

sigmaH = Math.sqrt((sumHH - howmany*avgH*avgH)/(howmany - 1));
sigmaW = Math.sqrt((sumWW - howmany*avgW*avgW)/(howmany - 1));
sigmaHW = Math.sqrt((sumHWHW - howmany*avgHW*avgHW)/(howmany - 1));
sigmaP = Math.sqrt((sumPP - howmany*avgP*avgP)/(howmany - 1));  

isay += '<tr>' + header + "Mean:" + tailer + header + formatit(avgW,1) + tailer;
isay += header + formatit(avgH,1) + tailer;

			isay += header + formatit(avgP,1) + tailer;
			isay += header + formatit(avgHW,0) + tailer + '<tr>';
isay += '<tr>' + header + "Target:" + tailer + header + XX + tailer + header + YY + tailer + header + Math.round(2*XX + 2*YY) + tailer + header + Math.round(XX*YY) + tailer + '</tr>';
 		
d.write(isay + '</table></center><BR>');
//=================
 isay = '<center><table style = "font-size:12 pt"><tr>' + header + '<sub>&nbsp;</sub>#<sup>&nbsp;</sup>' + tailer + header + 'Width(px)<sub>&nbsp;</sub><sup>&nbsp;</sup>' + tailer + header + '(W - W<sub>mean</sub>)&nbsp;<sup>&nbsp;</sup>' + tailer + header + '(W - W<sub>mean</sub>)<sup>2</sup>'+ tailer;
isay += header + '&nbsp;<sub>&nbsp;</sub>||<sup>&nbsp;</sup>&nbsp;' + tailer;
isay +=  header + 'Height(px)<sub>&nbsp;</sub><sup>&nbsp;</sup>' + tailer + header + '(H - H<sub>mean</sub>)&nbsp;<sup>&nbsp;</sup>' + tailer + header + '(H - H<sub>mean</sub>)<sup>2</sup>'+ tailer + '</tr>'; 

i = 1;
sumH = 0;
sumW = 0;
sumHH = 0;
sumWW = 0;
	for (i==1;i <= howmany;i++){

wd = parent.top.BB.cuts[i].ww;
ht = parent.top.BB.cuts[i].hh;
isay += '<tr>' + header + i + tailer + header + wd + tailer;
dummy1 = ablank;
dummy2 = ablank;
if (mode == "PRACTICE"){
		zz = wd - avgW;
		 dummy1 = formatit(zz,1);
		zz = zz*zz;
		 dummy2 = formatit(zz,1);
sumW +=  parseFloat(dummy1);
sumWW += parseFloat(dummy2); 
}
isay += header + dummy1 + tailer + header + dummy2 + tailer + header + '&nbsp;<sub>&nbsp;</sub>||<sup>&nbsp;</sup>&nbsp;' + tailer;

dummy1 = ablank;
dummy2 = ablank;
if (mode == "PRACTICE"){
			zz = ht - avgH;
		 dummy1 = formatit(zz,1);
			zz = zz*zz;
		 dummy2 = formatit(zz,1);
sumH +=  parseFloat(dummy1);
sumHH += parseFloat(dummy2);
}
isay += header + ht + tailer + header + dummy1 + tailer + header + dummy2 + tailer + '<tr>';
}

isay += '<tr><td COLSPAN = 8><hr width = "100%"></td></tr>';
isay += '<tr>' + header + "Sum:" + tailer + header + '&nbsp;';
dummy1 = ablank;
dummy2 = ablank;
dummy3 = ablank;
dummy4 = ablank;
if (mode == "PRACTICE"){
	dummy1 = formatit(sumW,1);
	dummy2 = formatit(sumWW,1);
dummy3 = formatit(sumH,1);
	dummy4 = formatit(sumHH,1);
}
isay += header + dummy1 + tailer + header + dummy2 + tailer + header + '&nbsp;<sub>&nbsp;</sub>||<sup>&nbsp;</sup>&nbsp;' + tailer + header + "&nbsp;" + tailer + header + dummy3 + tailer + header + dummy4 + tailer + '</tr>';
isay += '<tr><td COLSPAN = 8>&nbsp;</td></tr>';
dummy1 = ablank;
dummy2 = ablank;
dummy3 = ablank;
dummy4 = ablank;
if (mode == "PRACTICE"){
dummy1 = formatit(sumWW/(howmany - 1),1);
dummy2 = formatit(sumHH/(howmany - 1),1);
zz = Math.sqrt(sumWW/(howmany -1));
dummy3 = formatit(zz,1);
zz = Math.sqrt(sumHH/(howmany -1));
dummy4 = formatit(zz,1);
} 
isay += '<tr>' + header + '<sub>&nbsp;</sub><font face = "Symbol">s</font><sup>2</sup> = ' + tailer + '<td COLSPAN = 2 valign = bottom align = center><b>' + "Sum(W - W<sub>mean</sub>)<sup>2</sup>/(N - 1)" + tailer + header + '<sub>&nbsp;</sub><sup>&nbsp;</sup>' + dummy1 + tailer;



isay += header + '&nbsp;<sub>&nbsp;</sub>||<sup>&nbsp;</sup>&nbsp;' + tailer + '<td COLSPAN = 2 valign = bottom align = center><b>' + "Sum(H - H<sub>mean</sub>)<sup>2</sup>/(N - 1)" + tailer + header + '<sub>&nbsp;</sub><sup>&nbsp;</sup>' + dummy2 + tailer + '</tr>';

isay += '<tr>' + header + '<font face = "Symbol">s</font> = ' + tailer + '<td COLSPAN = 2 valign = bottom align = center><b>&nbsp;' + tailer + header + dummy3 + tailer;

isay += header + '&nbsp;<sub>&nbsp;</sub>||<sup>&nbsp;</sup>&nbsp;' + tailer + '<td COLSPAN = 2 valign = bottom align = center><b>&nbsp;'  + tailer + header + dummy4 + tailer + '</tr>';   
//======
dummy1 = ablank;
dummy2 = ablank;
if (mode == "PRACTICE"){
zz = Math.sqrt(sumWW/(howmany -1)/howmany);
dummy1 = formatit(zz,1);
zz2 = Math.sqrt(sumHH/(howmany -1)/howmany);
dummy2 = formatit(zz2,1);
} 
isay += '<tr>' + header + '<font face = "Symbol">s</font><sub>mean</sub> = ' + tailer + '<td COLSPAN = 2 valign = bottom align = center><b>' + '<font face = "Symbol">s</font>/sqrt(N)<sub>&nbsp;</sub>' + tailer + header + '<sub>&nbsp;</sub>' + dummy1 + tailer;

isay += header + '&nbsp;<sub>&nbsp;</sub>||&nbsp;&nbsp;' + tailer + '<td COLSPAN = 2 valign = bottom align = center><b>&nbsp;'  + tailer + header + '<sub>&nbsp;</sub>' + dummy2 + tailer + '</tr>';  


dummy1 = 0;
dummy3 = 0;
dummy4 = 0;
if (howmany < 11){
	dummy2 = parent.tfac[(howmany - 1)];
	}else{
		dummy2 = parent.tfactor(howmany - 1);
}
zz = dummy2*zz;
zz2 = dummy2*zz2;

if (mode == "PRACTICE"){
dummy1 = howmany -1;
}else{
	while (dummy1 != howmany - 1) dummy1 = prompt("How many degrees of freedom (df = N - 1) for the t-factor?","");

	parent.factor();

		dummy3 = prompt("What is the t-factor at the 95% Confidence Interval?","");
		if (dummy3 != dummy2) alert("Check closer, the value is " + dummy2);
}
		
isay += '<tr>' + header + 'df:' + tailer + '<td COLSPAN = 2 valign = bottom align = center><b>(degrees of freedom)' + tailer + header + dummy1 + tailer;
isay += header + '&nbsp;<sub>&nbsp;</sub>&nbsp;&nbsp;&nbsp;' + tailer + '<td COLSPAN = 2 valign = bottom align = center><b>t factor (95%):'  + tailer + header  + dummy2 + tailer + '</tr>';  


dummy1 = ablank;
dummy2 = ablank;
if (mode == "PRACTICE"){
	dummy1 = formatit(zz,1);
	dummy2 = formatit(zz2,1);
}
isay += '<tr>' + header + '<font face = "Symbol">d</font><sub>95%</sub> = ' + tailer + '<td COLSPAN = 2 valign = bottom align = center><b>' + 't-factor x <font face = "Symbol">s</font><sub>mean</sub> = ' + tailer + header + '<sub>&nbsp;</sub>' + dummy1 + tailer;

isay += header + '&nbsp;<sub>&nbsp;</sub>||&nbsp;&nbsp;' + tailer + '<td COLSPAN = 2 valign = bottom align = center><b>&nbsp;'  + tailer + header + '<sub>&nbsp;</sub>' + dummy2 + tailer + '</tr>';  
//========

dummy1 = ablank;
dummy2 = ablank;
if (mode == "PRACTICE"){
	dummy1 = 2*Math.sqrt(zz*zz + zz2*zz2);
	dummy2 = avgH*avgW*Math.sqrt(zz*zz/(avgW*avgW) + zz2*zz2/(avgH*avgH));
dummy1 = formatit(dummy1,1);
dummy2 = Math.round(dummy2);
}
isay += '<tr>' + header + '(<font face = "Symbol">d</font><sub>95%</sub>)<sub>peri</sub> = ' + tailer + '<td COLSPAN = 2 valign = bottom align = center><b>' + '&nbsp;' + tailer + header + '<sub>&nbsp;</sub>' + dummy1 + tailer;

isay += header + '&nbsp;<sub>&nbsp;</sub>&nbsp;&nbsp;&nbsp;' + tailer + header + '(<font face = "Symbol">d</font><sub>95%</sub>)<sub>area</sub> = ' + tailer  + header + '&nbsp;'  + tailer + header + '<sub>&nbsp;</sub>' + dummy2 + tailer + '</tr>';  
//========
isay += '<tr><td COLSPAN = 8><hr width = "100%"></td></tr>';

d.write(isay + '</table></center>');
if (mode != "PRACTICE"){
isay =   "<center><b>perimeter = ___________ ± __________ px (95% level of confidence)<BR><BR>";
isay += "area = ___________ ± __________ px<sup>2</sup> (95% level of confidence)</b></center>"
d.write(isay);
}
d.close();

}

function formatit(AA,BB){
var LL = 0;
var XX = 0;
var mult = " ";
if (AA < 0){
		AA = Math.abs(AA);
		mult = "-"
}
if (BB > 0){
	AA = Math.round(AA*Math.pow(10,BB)).toString();
	LL = AA.length;
if (LL > 1){
	XX = AA.substring(0,LL-1) + "." + AA.substring(LL-1);
	}else{
	XX = "0." + AA;
	}
}else{
	XX = Math.round(AA);
}

	return  mult + XX;
}