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

A=[1, 1, 1;1, 2, 4;1, 3, 9;1, 4, 16];

#Obs: para contestar lo que se pregunta no es necesario
#calcular las transformaciones.

#Calcular H1
x1=A(:,1);
e1=zeros(4,1);
e1(1)=1;
u1=x1-sign(x1(1))*norm(x1)*e1;
u1=u1/norm(u1);

H1=eye(4)-2*u1*u1';
A2=H1*A

#Calcular H2
x2=A2(2:end,2);
e2=zeros(3,1);
e2(1)=1;
u2=x2-sign(x2(1))*norm(x2)*e2;
u2=u2/norm(u2);

H2_=eye(3)-2*u2*u2';
H2=eye(4);
H2(2:end,2:end)=H2_;

A3=H2*A2