% Strang 2.4/26 outer products
% June 2007, M. Kawski
%
A=[1 0;2 4; 2 1]
B=[3 3 0;1 2 1]
AB=A*B
pause
S=zeros(3,3);
for i=1:2
    A(:,i)
    B(i,:)
    N=A(:,i)*B(i,:)
    pause
    S=S+N;
end
sum_of_all_4=S
compare_AB=AB