X = [2,3,4 ; 3,3,4 ; 4,4,4 ; 5,4,4 ; 5,4,5 ; 4,4,5 ; 3,3,5 ; 2,3,5 ; 2,3,4];
Xprim = [2 0 0 0; 0 2 0 0; 0 0 2 0; 0 0 0 1];
Xprim2 = [1/2 0 0 0; 0 1/2 0 0; 0 0 1/2 0; 0 0 0 1];

plot3(X(:,1), X(:,2), X(:,3));
hold on

xj = [X ones(9,1)];

res = xj * Xprim;
plot3(res(:,1), res(:,2), res(:,3));

res2 = xj * Xprim2;
plot3(res2(:,1), res2(:,2), res2(:,3));