P=[1,2,3];
Q=[2,4,-1];
u=Q-P;
u1=[u(1),0,0];
u2=[0,u(2),0];
u3=[0,0,u(3)];
nu1=norm(u1);
nu2=norm(u2);
nu3=norm(u3);


PQ=P+u1;
PQQ=PQ+u2;
PQQQ=Q;
murtoviiva=[P',PQ',PQQ',Q'];
plot3(murtoviiva(1,:),murtoviiva(2,:),murtoviiva(3,:))
hold on
plot3(P(1),P(2),P(3),'*')
plot3(Q(1),Q(2),Q(3),'*')