// 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)
