Harj. 3 loppuviikko, teht. 4
============================

> with(linalg):
Warning: new definition for   norm
Warning: new definition for   trace

--------------------------------------------------------------------------------
> v1:=matrix(2,1,[1,1]);

                                         [ 1 ]
                                   v1 := [   ]
                                         [ 1 ]
--------------------------------------------------------------------------------
> v2:=matrix(2,1,[2,-1]);

                                        [  2 ]
                                  v2 := [    ]
                                        [ -1 ]
--------------------------------------------------------------------------------
> Fv1:=matrix(2,1,[9,-3]);

                                         [  9 ]
                                  Fv1 := [    ]
                                         [ -3 ]
--------------------------------------------------------------------------------
> Fv2:=matrix(2,1,[3,6]);

                                         [ 3 ]
                                  Fv2 := [   ]
                                         [ 6 ]
--------------------------------------------------------------------------------
> ##Annettu:  A v1 = Fv1, A v2 = Fv2;
--------------------------------------------------------------------------------
> ## Talloin:  A [v1 | v2]=[A v1 | A v2] = [Fv1 | Fv2]
> # 
> V:=augment(v1,v2);

                                      [ 1   2 ]
                                 V := [       ]
                                      [ 1  -1 ]
> FV:=augment(Fv1,Fv2);

                                       [  9  3 ]
                                 FV := [       ]
                                       [ -3  6 ]
--------------------------------------------------------------------------------
> Siis:
--------------------------------------------------------------------------------
> A:=evalm(FV&*inverse(V));

                                      [ 4   5 ]
                                 A := [       ]
                                      [ 1  -4 ]
--------------------------------------------------------------------------------
> # Tarkistus:
> evalm(A&*v1);

                                     [  9 ]
                                     [    ]
                                     [ -3 ]
--------------------------------------------------------------------------------
> evalm(A&*v2);

                                      [ 3 ]
                                      [   ]
                                      [ 6 ]
--------------------------------------------------------------------------------
> # No oikeinhan se!
>