clear all; close all; clc;
% REPASO DE TÉCNICAS NUMÉRICAS

% Supongamos que tenemos que resolver un sistema de ecuaciones lineales Ax+Bu=0 . Por ej.el modelo de una columna de 
% absorción de 5 etapas da lugar a un sistema de ese tipo. Sea

A=[-.325 .125 0 0 0;.2 -.325 .125 0 0;
0 .2 -.325 .125 0;0 0 .2 -.325 .125;0 0 0 .2 -.325];

B=[.2 0;0 0;0 0;0 0; 0 .25];

u=[0 .1]';

% para resolver el sistema puede utilizarse el comando "inv" para invertir la matriz A:

x = -inv(A)*B*u;
disp(x);

pause
% pero también puede utilizarse el operador "\"
x1 = -A\(B*u);
disp(x1);