P=[0,0,2];
A=[1,0,0];
B=[0,0,1];
yli=1;
suora1=[(A-yli*(B-A))',(A+(1+yli)*(B-A))'];
pisteet=[P',A',B'];
plot3(pisteet(1,:),pisteet(2,:),pisteet(3,:),'*')
hold on, axis equal
plot3(suora1(1,:),suora1(2,:),suora1(3,:))

u=B-A;
uhat=u/norm(u);
v=A-P;
v1=(v*uhat')*uhat;
v2=v-v1;
etaisyys=norm(v2);

%vektorit
vek{1}=[A',B'];nimi{1}='u';
vek{2}=[P',A'];nimi{2}='v';
vek{3}=[P',(P+v1)'];nimi{3}='v1';
vek{4}=[P',(P+v2)'];nimi{4}='v2';

for k=1:4
plot3(vek{k}(1,:),vek{k}(2,:),vek{k}(3,:))
text(mean(vek{k}(1,:)),mean(vek{k}(2,:)),mean(vek{k}(3,:)),nimi{k})
end

%pisteet
pos{1}=P;nimi{1}='P';
pos{2}=A;nimi{2}='A';
pos{3}=B;nimi{3}='B';

for k=1:3
plot3(pos{k}(1),pos{k}(2),pos{k}(3),'*')
text(pos{k}(1),pos{k}(2),pos{k}(3),nimi{k})
end