function [XYZ2] = T12(r,T_21)

r
s=size(XYZ1)
n=s(1)

XYZ2=[]
 for i=4:n
    a=XYZ(i)
    b=T_21*a';
    XYZ2=[XYZ2;b]
 end
end

