1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
|
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') ;
|
