function K = transformation_mat(th)

K = 2/3*[ cos(th)   cos(th-2*pi/3)   cos(th+2*pi/3)
          sin(th)   sin(th-2*pi/3)   sin(th+2*pi/3)
            1/2          1/2              1/2       ];