
%%piirretään akselit
lmin=-1;lmax=1;
plot3([lmin,lmax],[0,0],[0,0],'k',"linewidth",1),hold on
plot3([0,0],[lmin,lmax],[0,0],'k',"linewidth",1)
plot3([0,0],[0,0],[lmin,lmax],'k',"linewidth",1)

%%nimetään akselit
text(lmax,0,0,['x'])
text(0,lmax,0,['y'])
text(0,0,lmax,['z'])
view(120,30)

%%tehdään vektorit
u=[0.4,0.2,0]';
v=[0.1,0.5,0.05]';
w=cross(u,v);

%%siirretään vektorit pisteeseen p
p=[1,1,0]';
pu=[p,p+u];
pv=[p,p+v];
pw=[p,p+w];

%%piirretään p
t=plot3(pu(1,:),pu(2,:),pu(3,:),"linewidth",2,'color','red')
text(pu(1,end),pu(2,end),pu(3,end),['u'],'color','red')
hold on
plot3(pv(1,:),pv(2,:),pv(3,:),'r',"linewidth",2,'color','green')
text(pv(1,end),pv(2,end),pv(3,end),['v'],'color','green')
plot3(pw(1,:),pw(2,:),pw(3,:),'r',"linewidth",2,'color','blue')
text(pw(1,end),pw(2,end),pw(3,end),['w'],'color','blue')

%%tarkastelusuuntaa voisi muuttaa näin
%s=0;
for s=1:6
az=120+15*(cos(s));el=40+15*(sin(s));view(az,el);
end
Ristitulo
