P7 - Gram Schmidt Process.sce

clc;clear;
A= [0 0 1;0 1 1;1 1 1];
disp(A, "The given matrix A is:");
[m,n] = size(A);
for k=1:n
    V(:,k) = A(:,k);
    for j=1:k-1
        R(j,k) = V(:,j)'*A(:,k);
        V(:,k) = V(:,k)-R(j,k)*V(:,j);
    end
    R(k,k) = norm(V(:,k));
    V(:,k) = V(:,k)/R(k,k);
end
disp(V,"Q: ");