u=[-10,2,8]
a=[3,12,-1]
disp("Vector u:",u)
disp("Vector a:",a)
disp("u*a=",u*a')
disp("norm(a)^2:",norm(a)^2);
disp("Projection of u along a:",((u*a')/(norm(a)^2))*a);