theta=0;
r=[1 0 0; 0 cos(theta) -sin(theta);0 sin(theta) cos(theta)];
a=r*[2;0;0]
b=r*[0;2;0]
c=r*[0;0;2] 
