public static double[][] powerSymmMatrix(
	double[][] inMatrix,
	double power )
{
	EigenValueDecompositionSymm eigenDeco =
		new EigenValueDecompositionSymm(inMatrix);
	int m = Matrix.getNumOfRows(inMatrix);
	double[][] eigenVectors = eigenDeco.getEigenVectors();
	double[] eigenValues = eigenDeco.getEigenValues();
	for ( int i = 0; i < m; ++i )
		eigenValues[i] = Math.pow(eigenValues[i], power);
	return(
		Matrix.mult(
			Matrix.mult(eigenVectors, Matrix.diag(eigenValues)),
			Matrix.transpose(eigenVectors)));
}