disp('Matrix-vector product')
clear;
zero_mat = zeros(3)
pause
identity_mat = eye(4)
pause
A = rand(5)
pause
disp('2 A =')
disp(2*A)
pause
x = [1 2 3 4 5]'
pause
c = 3
pause
disp('(c A) x = ')
disp((c*A)*x)
pause
disp('A (c x) = ')
disp(A*(c*x))
pause
disp('An elementary matrix')
pause
E21 = eye(5)
pause
E21(2,1) = 1
pause
E21_x = E21*x
pause
disp('A block matrix')
pause
A11 = [2 -1;-1 2]
pause
A12 = zeros(2,3)
pause
A21 = A12'
pause
A22 = [2 -1 0;-1 2 -1;0 -1 2]
pause
disp('AA = [A11 A12;A21 A22]')
disp([A11 A12;A21 A22])