xy =[1 ; 1]
phi = 3.14159/4
R = [cos(phi) -sin(phi) ; sin(phi) cos(phi)]
xyc = R*xy