#Métodos Numéricos - 2023. IMERL, Fing
#Código de ejemplo del ejercicio 6 del práctico 2

A=10.*rand(6)-5;
B=A;
B(:,4)=ones(6,1);

u=A(:,4)-ones(6,1);
v=zeros(6,1);
v(4,1)=1;

Ainv=inv(A);
c=v.'*Ainv*u;
alfa=1/(1-c);

Binv=Ainv+alfa.*((Ainv*u*v.'*Ainv));

%Verifico
disp('Verificación')
sum(sum(B*Binv-eye(6,6)))
sum(sum(Binv*B-eye(6,6)))