% Tämä koodi piirtää vektorin projektion toisen vektorin virittämän 
% aliavaruuden suhteen.

% Dollarimerkkien perään voi kirjoittaa kommentteja, jotka eivät vaikuta
% koodiin.

% Määritellään projisoitava vektori.
v = [3 ; 1];
% Määritellään vektori jonka aliavaruuden suhteen projektio otetaan.
w = [1;2];
% Lasketaan projektio.
vw = dot(v,w);
ww = dot(w,w);
projv = (vw/ww)*w;

% Seuraava koodinpätkä piirtää vektorit v ja projv. Sinun ei tarvitse 
% ymmärtää sen kaikkia yksityiskohtia.
figure(1)
clf
% Piirretään vektorit v ja projv.
plot([0 v(1)],[0 v(2)],'r','LineWidth',3)
hold on
plot([0 projv(1)],[0 projv(2)],'b','LineWidth',3)
% Kirjoitetaan kuvaan vektorien nimet.
legend('v','proj_w(v)','Location','northwest')
% Piirretään koordinaattiakselit.
plot([-8 8],[0,0],'--k')
hold on
plot([0 0],[-6,6],'--k')
hold on
% Piirretään projektiosuora
plot([-6*w(1) 6*w(1)],[-6*w(2),6*w(2)],':b')
hold on
% Määritellään akseleiden pituudet.
axis([-6 6 -6 6])
axis equal