m=[1 0;0 4];
k=[3 -1;-1 1];
newk=m^(1/2)*k*m^(-1/2)
