function testhouseholder()
    A = [1;2;3;4;5]
    B = [6;7;8;9;10]
    u = householder(A,B)
    disp(u)
endfunction