implemented householder transformation in both C++ and octave/Matlab