script T3 CN

clc
clear

m = input('Nr. linii: ');
n = input('Nr. coloane: ');

A = zeros(m, n);
rel = 1;
while(rel)
    while(rel)
        for i = 1: m
            fprintf('Linia %g: ', i);
            A(i,:) = input(' ');
        end
        rel = input('Introduceti 1 pentru reluare citire sau 0 pentru continuare: ');
    end
    R = A;
    Q = eye(m);
    for k = 1: n
        x = zeros(m, 1);
        x(k : m, 1) = R(k : m, k);
        nx = norm(x);
        u = x;
        u(k) = x(k) + nx;
        nu = norm(u);
        if nu ~= 0
            v = u / nu;
            w = 2 * R' * v;
            R = R - v * w';
            Q = Q - 2 * Q * v * v';
        end
    end
    Q
    R
    rel = input('Introduceti 1 pentru reluare program sau 0 pentru inchidere: ');
end