//@code
objref m, q, e
m = new Matrix(5,5)
m.setdiag(0, 2)
m.setdiag(-1, -1)
m.setdiag(1, -1)
m.printf

q = new Matrix(1,1)
e = m.symmeig(q)
print "eigenvectors"
q.printf

print "eigenvalues"
e.printf

print "qt*m*q"
q.transpose.mulm(m).mulm(q).printf

print "qt*q"
q.transpose.mulm(q).printf