setSlice(1);
width = getWidth();
height = getHeight();

if (nSlices>2){
setBatchMode(true);
	for (z=0; z<nSlices; z++){
		/////////////////////////////////////
		acc1 = 0.0;
		sigNum = 0.0;

		for (y=0; y<height; y++) {
			for (x=0; x<width; x++) {
				i = getPixel(x,y);
				if (i!=0){
					sigNum = sigNum + 1;
					acc1 = acc1 + i;
				}
			}
		}
		//print(acc1);
		//print(sigNum);

		mean = acc1/sigNum;
		//print(mean);
		//""+acc1;
		/////////////////////////////////////
		acc2 = 0.0;

		for (y=0; y<height; y++) {
			for (x=0; x<width; x++) {
				i = getPixel(x,y);
				sa = i - mean;
				if (i!=0){acc2 = acc2 + pow(sa, 2);}
			}
		}
		variance = acc2/sigNum;
		stdev = sqrt(variance);
		stdev3 = pow(stdev, 3);
		//print(stdev);
		/////////////////////////////////////
		acc3 = 0.0;

		for (y=0; y<height; y++) {
			for (x=0; x<width; x++) {
				i = getPixel(x,y);
				sa = i - mean;
				if (i!=0){acc3 = acc3 + pow(sa, 3);}
			}
		}
		skewness1 = acc3/stdev3;
		skewness2 = skewness1/sigNum;

		//print(skewness2);
		/////////////////////////////////////
			
		row = nResults();
		setResult("sigNum", row, sigNum);
		setResult("skewness", row, skewness2);
		updateResults();
		run("Next Slice [>]");
	}
setBatchMode(false);

}else {
/////////////////////////////////////
acc1 = 0.0;
sigNum = 0.0;
for (y=0; y<height; y++) {
	for (x=0; x<width; x++) {
		i = getPixel(x,y);
		if (i!=0){
			sigNum = sigNum + 1;
			acc1 = acc1 + i;
		}
	}
}
//print(acc1);
//print(sigNum);
mean = acc1/sigNum;
//print(mean);
//""+acc1;
/////////////////////////////////////
acc2 = 0.0;
for (y=0; y<height; y++) {
	for (x=0; x<width; x++) {
		i = getPixel(x,y);
		sa = i - mean;
		if (i!=0){acc2 = acc2 + pow(sa, 2);}
		}
	}
variance = acc2/sigNum;
stdev = sqrt(variance);
stdev3 = pow(stdev, 3);
//print(stdev);
/////////////////////////////////////
acc3 = 0.0;
for (y=0; y<height; y++) {
	for (x=0; x<width; x++) {
		i = getPixel(x,y);
		sa = i - mean;
		if (i!=0){acc3 = acc3 + pow(sa, 3);}
	}
}
skewness1 = acc3/stdev3;
skewness2 = skewness1/sigNum;
//print(skewness2);
/////////////////////////////////////
			
row = nResults();
setResult("sigNum", row, sigNum);
setResult("skewness", row, skewness2);
updateResults();

}