//page 156
clear;
close;
clc;
a=[1;1;1];
disp(a,'a=');
P=(a*a')/(a'*a);
disp(P,'Matrix that projects onto a line through a=(1,1,1) is');
//end