function SPDmaticeTransformujiKruhNaElipsu()

A = rand(2,2); A = 0.5*(A+A');
[V,Lambda] = eig(A);

t = 0:0.01:2*pi; x = [cos(t);sin(t)]; Ax = A*x;
plot(x(1,:),x(2,:),'r',Ax(1,:),Ax(2,:),'b')
axis equal
hold on; plot([0,V(1,1)],[0,V(2,1)],'k')
hold on; plot([0,V(1,2)],[0,V(2,2)],'k')