> #===============================================
> # Testaufgabe 4
> #===============================================
> restart; with(linalg):
Warning, new definition for norm
Warning, new definition for trace
> m:=3; n:=2; A:=matrix(m,n, [3,7,0,12,4,1]);
>
> read ` givens.src`;
Spalte 1, Zeile 3: tau= .75000, s= .80000, c= .60000
R: 1. und 3. Zeile
5.00000 5.00000
0.00000 -5.00000
Q': 1. und 3. Zeile
.60000 0.00000 .80000
-.80000 0.00000 .60000
Spalte 2, Zeile 3: tau= -.41667, c= .92308, s= -.38462
R: 2. und 3. Zeile
13.00000
0.00000
Q': 2. und 3. Zeile
.30769 .92308 -.23077
-.73846 .38462 .55385
> multiply(Q,R); # Test
>