Ingenieria

! % ! #& " '! , ! ! !( / ! ! ! " ' ! !( # 1 $ ./ ! #! $ ! ! ! 0 (

!

!"

#!

!

!

$

#

$ # $ $ , !" ! ! 1(! - / $ 1 ! $ ( $ / #! $

! ) ! *+ ( ' , .% /" ( ! / ' (! ! * $ $ - ' !(! $ , " $ $ $ !1( $ #! %

$

2% - ! 3 % ( $ $ $ $ ! *!0 & ./ " 1( 4(5 " * ( $ ! ! '($ $ ! .' $ (5 & ) ! % ( $ & ( ! " " ' # !! ! 4( $ ! & & ! % 6% - 2 ! , - # # ! & $ !(' !1( & ! ( ! ! 4( ( ! ! $$( $ / 4( ! - 6 ,( $ ! !( - $ ! ! $ ! 5 $ $ " ! * 5 $ $ + ( ', .% 7% - 7 1 " & # 3 ! $ ! ! 1 $ $ ( ' 4( $ # 0 ( ! ! !( 8 #! / 9 $$ ! $ $ ! % ! $ $ ( &( $ ,( $ with()% 1( $ $ ! 4( &( $ ,( 0, *$ ! ,(!1 $ ! % 0, $ 1 !1 # ! !% . 4( $ # 1 $ 0 $ # ! 4( $ ! /9 '" ( $

!

!

!

( - : > restart; with(linalg):
Warning, the protected names norm and trace have been redefined and unprotected

$4(

" !1 8

/ ! 4( ! ( $ $ $ ! !( 4( % >

! 4(

" # $ % $ ( 9 ,( " 1( ( (

"

# ! $

$ ! $ 4( 4(0 " ) ! !1 # ! !'$ & ! '($ $ ! !1 !

%&
$ ! '(
( $ , ( ( ! ! $ matrix() ' ! $ , ! $! *! ! $ ! ,! $ ! .' > A:= matrix([[1,-alpha,2/3],[-1,0,1],[beta/3,2,-1]]); 2 −α 6 7 2

= :; matrix(3,3,[1,-alpha,2/3,-1,0,1,beta/3,2,-1]); 2 A[3,1]; $ )*+$ ( −α = 6 / 6 7 2

.%

"
" $ ( 4(' ,! /( # !( !

'
(

"
$ , ! % $ (

! *

*"
' ./ " ! (! β 7 6 addcol(A1,1,2,-1); 2 = I > addrow(A1,2,3,-9); 2 = I > solve({x1-x2+x3=1,4*x1-5*x2-5*x3=4,2*x1+x2-x3=2,x1+2*x2-2*x3=1},{x1,x2,x3}); = =/ = 2 N # 1 /! $ ( ! ! #5 ( $ ( ! ( ! !% ! ! $ linsolve(A,B)/ 4( ! 4( !" $ ! ( ! = : > A:= matrix([[1,2,2],[1,1,3],[-2,0,1]]); M 2 :; 2 B:=matrix([[3],[2],[1]]); :; > linsolve(A,B);B:=col(Am,4); > linsolve(A,B);

9

! ! )!

2 B:=col(C,4); A:=delcols(C,4..4); :; @ 2/ =/ =/ = A 2 = :; = = > backsub(A,B); ! 1 " # P (#% > / ! $ @ 2/ =/ = A $ forwardsub(A,b) ( !" ! ( $ /' ( ! FO;B 31( ! 4( !

4
,

5

!
! !" " " 1 ( /! , (! #! $ ( - $% $ ( ( &( $ , ' 5 2 = 2/ 6 = 6/%%%/ ) ! $ ! ! 1(! .: '$%( = / $ )=% $) ' ./ 4( 0/ ( " )=% 1(! ( " 4( , ( $ = / % ! $ ! $ ! : 5 $

' "...