Sponsored links: Algebra eBooks
 

Related

get_lu_factors-matrix-rat

m: matrix([e1,2*e2,e3...

rat(get_lu_factors(lu...

Calculate

get_lu_factors-matrix

m: matrix([e1,e2],[1,...

get_lu_factors(lu_fac...

Calculate

get_lu_factors-matrix-rat

m1: matrix([e1,2*e2,e...

m2: matrix([2*e1,e2,e...

rat(get_lu_factors(lu...

Calculate

get_lu_factors-matrix-rat

m1: matrix([e1,2*e2,e...

m2: matrix([2*e1,e2,e...

rat(get_lu_factors(lu...

Calculate

get_lu_factors-matrix-rat

m: matrix([e1,2*e2,e3...

rat(get_lu_factors(lu...

Calculate

get_lu_factors-matrix

m: matrix([e1,e2],[1,...

get_lu_factors(lu_fac...

Calculate

get_lu_factors-matrix-rat

m1: matrix([e1,2*e2,e...

m2: matrix([2*e1,e2,e...

rat(get_lu_factors(lu...

Calculate

get_lu_factors-matrix

m: matrix([e1,e2],[1,...

get_lu_factors(lu_fac...

Calculate

get_lu_factors-matrix-rat

m: matrix([e1,2*e2,e3...

rat(get_lu_factors(lu...

Calculate

get_lu_factors-matrix

m: matrix([e1,e2],[1,...

get_lu_factors(lu_fac...

Calculate

get_lu_factors

Run Example
(%i1)A:matrix([2,2,1],[1,1,2],[2,-1,1]);
                                 [ 2   2   1 ]
                                 [           ]
(%o1)                            [ 1   1   2 ]
                                 [           ]
                                 [ 2  - 1  1 ]
(%i2) lu_factor(A);
                     [ 2   2   1 ]
                     [           ]
                     [ 1       3 ]
(%o2)               [[ -   0   - ], [1, 3, 2], generalring]
                     [ 2       2 ]
                     [           ]
                     [ 1  - 3  0 ]
(%i3) get_lu_factors(lu_factor(A));
                                 [ 1  0  0 ]  [ 2   2   1 ]
                    [ 1  0  0 ]  [         ]  [           ]
                    [         ]  [ 1  1  0 ]  [ 0  - 3  0 ]
(%o3)              [[ 0  0  1 ], [         ], [           ]]
                    [         ]  [ 1       ]  [         3 ]
                    [ 0  1  0 ]  [ -  0  1 ]  [ 0   0   - ]
                                 [ 2       ]  [         2 ]
(%i4) 
Run Example
A:matrix([2,2,1],[2,-1,1],[1,1,2]);
                                 [ 2   2   1 ]
                                 [           ]
(%o1)                            [ 2  - 1  1 ]
                                 [           ]
                                 [ 1   1   2 ]
(%i2) lu_factor(A);
                     [ 2   2   1 ]
                     [           ]
                     [ 1  - 3  0 ]
(%o2)               [[           ], [1, 2, 3], generalring]
                     [ 1       3 ]
                     [ -   0   - ]
                     [ 2       2 ]
(%i3) get_lu_factors(lu_factor(A));
                                 [ 1  0  0 ]  [ 2   2   1 ]
                    [ 1  0  0 ]  [         ]  [           ]
                    [         ]  [ 1  1  0 ]  [ 0  - 3  0 ]
(%o3)              [[ 0  1  0 ], [         ], [           ]]
                    [         ]  [ 1       ]  [         3 ]
                    [ 0  0  1 ]  [ -  0  1 ]  [ 0   0   - ]
                                 [ 2       ]  [         2 ]
(%i4) 
Run Example
block(a: matrix ( [1,3,4],[2,3,7],[6,8,7]),A:lu_factor(a),L:get_lu_factors(A)[2] , U:get_lu_factors(A)[3],print("deco is" , L ,U));
        [ 1  0   0 ] [ 1   3    4   ]
        [          ] [              ]
        [ 2  1   0 ] [ 0  - 3  - 1  ]
deco is [          ] [              ] 
        [    10    ] [           41 ]
        [ 6  --  1 ] [ 0   0   - -- ]
        [    3     ] [           3  ]
                               [ 1   3    4   ]
                               [              ]
                               [ 0  - 3  - 1  ]
(%o1)                          [              ]
                               [           41 ]
                               [ 0   0   - -- ]
                               [           3  ]
(%i2) 

Related Help

Help for Get_lu_factors