octave:9> format compact octave:10> octave:10> octave:10> octave:10> x = randn(6,1) x = -1.23105 1.13420 0.36789 0.98213 0.20994 0.30009 octave:11> [v, bet] = house1(x); octave:12> octave:12> P = eye(6)-bet*v*v'; octave:13> P*x ans = 2.00897 -0.00000 0.00000 0.00000 -0.00000 0.00000 octave:14> norm(x) ans = 2.0090 octave:15> octave:15> [v, bet] = house(x); octave:16> P = eye(6)-bet*v*v'; octave:17> P*x ans = 2.00897 -0.00000 0.00000 0.00000 -0.00000 0.00000 octave:18> x = randn(6,1) x = -1.015781 0.801436 -0.877283 -1.977747 0.093801 0.784859 octave:19> [v, bet] = house(x); octave:20> P = eye(6)-bet*v*v'; octave:21> P*x ans = 2.64197 -0.00000 -0.00000 0.00000 -0.00000 0.00000 octave:22> ); parse error: syntax error >>> ); ^ octave:22> octave:22> x = randn(6,1) x = 1.12069 0.26540 0.84943 -0.92029 -2.56679 1.16383 octave:23> [v, bet] = house(x); octave:24> P = eye(6)-bet*v*v'; octave:25> P*x ans = -3.29207 0.00000 -0.00000 -0.00000 0.00000 0.00000 octave:26> octave:26> octave:26> X = randn(8,3) X = -0.6457719 -0.4002955 -0.5434182 -0.4892368 0.2244585 -1.2836340 0.3403562 -1.1364485 1.0754215 -0.5828899 -0.0028010 0.5452226 -0.3667113 -0.8529507 0.1246327 0.1363842 0.6273785 -0.8192728 1.6491485 0.1126639 0.5463845 -0.2559651 0.5470596 -0.9272934 octave:27> [B, bet] = hoQR(X); octave:28> v1 = [1; B(2:8,1)]; octave:29> P1 = eye(8) - bet(1)*v1*v1'; octave:30> X1 = P1*X X1 = 2.0125e+00 1.0318e-01 9.9782e-01 5.8981e-17 3.1712e-01 -9.9998e-01 7.7195e-17 -1.2009e+00 8.7809e-01 1.2143e-17 1.0760e-01 8.8317e-01 7.0256e-17 -7.8350e-01 3.3725e-01 4.1200e-18 6.0155e-01 -8.9835e-01 -3.4001e-16 -1.9968e-01 -4.0976e-01 5.5511e-17 5.9554e-01 -7.7889e-01 octave:31> v2 = [0; 1; B(3:8,2)]; octave:32> P1 = eye(8) - bet(2)*v2*v2'; octave:33> X2 = P2*X1 error: 'P2' undefined near line 1 column 6 octave:33> X2 = P1*X1 X2 = 2.0125e+00 1.0318e-01 9.9782e-01 -1.4216e-17 1.7101e+00 -1.4404e+00 1.4094e-17 -1.6653e-16 4.9843e-01 1.7797e-17 1.7347e-17 9.1719e-01 2.9088e-17 -2.2204e-16 8.9546e-02 3.5728e-17 3.1919e-16 -7.0817e-01 -3.5050e-16 1.0408e-17 -4.7289e-01 8.6804e-17 1.1102e-16 -5.9061e-01 octave:34> octave:34> octave:34> octave:34> [Q R] = qr(X); octave:35> Q Q = Columns 1 through 5: -0.32087622 0.21471225 -0.36135576 -0.13669717 -0.43611274 -0.24309579 -0.14591762 -0.56383419 0.50706305 0.02180536 0.16911884 0.67473711 -0.04424877 0.38942702 -0.41578776 -0.28963092 -0.01583607 0.58157163 0.65733140 0.11378029 -0.18221439 0.48776632 -0.26880309 -0.02056001 0.77834093 0.06776767 -0.36276860 -0.24725551 0.22359264 0.01756800 0.81944179 -0.01644146 -0.16800864 0.25126030 0.12954323 -0.12718592 -0.32756433 -0.22296235 0.16405747 -0.02449689 Columns 6 through 8: 0.18148288 0.69176752 -0.00057291 -0.42073667 -0.13799196 -0.38477780 0.12756314 -0.37257959 0.18587611 0.07434480 0.35648254 0.00033012 0.16545716 0.06700665 0.13610126 0.84584279 -0.15186983 -0.11770486 -0.10551160 0.45645896 0.02484487 -0.11952542 -0.02472735 0.88566383 octave:36> R R = 2.01253 0.10318 0.99782 0.00000 -1.71014 1.44038 0.00000 0.00000 1.47365 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 octave:37> diary off