LU factorization for matrices in quasiseparable form via orthogonal transformations