I'm trying to find out wether a matrix is orthonormal. I begin by checking if the vectors are normal by doing

for j=1:2;
    if norm(S:,j) ~= 1;
        return; % Not normal vector
    end
end

But when norm returns 1.0000 comparing that to 1 is true and the function returns, which is not what i want. Any ideas?

Thx