

format long e
fprintf ('-------------------------------------------------\n') ;
x = [-3 4 5]'
[v, beta, s] = house (x) 
x - v*(beta*(v'*x))
fprintf ('-------------------------------------------------\n') ;


x = [3 4 5]'
[v, beta, s] = house (x) 
x - v*(beta*(v'*x))
fprintf ('-------------------------------------------------\n') ;

x = [1 eps]'
[v, beta, s] = house (x) 
x - v*(beta*(v'*x))
fprintf ('-------------------------------------------------\n') ;

x = pi
[v, beta, s] = house (x) 
x - v*(beta*(v'*x))
fprintf ('-------------------------------------------------\n') ;

x = -pi
[v, beta, s] = house (x) 
x - v*(beta*(v'*x))
fprintf ('-------------------------------------------------\n') ;

x = [1 0 0]'
[v, beta, s] = house (x) 
x - v*(beta*(v'*x))
fprintf ('-------------------------------------------------\n') ;
