%
% Experimentando con matrices.
%

fprintf('\nMatrices. Ejemplo 1\n');

A = [ 1 0 -1; 0 1 0; 2 0 -1 ]

d = det(A)

B = A'

C = A * B
D = A .* B

fprintf('\nMatrices. Ejemplo 2\n');

v = [ 10; 20; 30 ];

u = A * v

% Producto en forma explícita

for i = 1:3
  u(i) = 0;
  for j = 1:3
    u(i) = u(i) + A(i, j) * v(j);
  end
end

u

fprintf('\nMatrices. Ejemplo 3\n');

% Sistemas de ecuaciones lineales (pequeños).
%
%    A X = b
%
% Si det(A) no es cero, entonces X = A^(-1) * b
% es decir: X = A\b  (división "a la izquierda")
%
% Ej:
%      x + y = 2
%      x - y = -2
%
%      Solución x = 0, y = 2.

A = [ 1 1; 1 -1];
b = [2; -2];

X = A \ b


