5B1132 Amelia 1 för P ht 2003

Att lösa ekvationssystem med maple

Här är några tips för den som vill lösa linjära ekvationssystem med datorprogrammet maple. Starta maple t ex genom att skriva xmaple i terminalfönstret på din sun-station. Börja sedan maple-sessionen med att skriva (obs håll koll på kolon, semikolon och komman):

>with(plots);

och tryck på retur. Skriv sedan

>with(LinearAlgebra);

och tryck på retur. Sedan kan ni prova nedanstående kommandon och se vad som händer. Missa inte semikolon, men prompten behöver ni förstås inte skriva. Avsluta raden med att trycka på retur. Därefter kan ni modifiera inmatningarna så att de passar till Grupparbete 1.

Man kan alltid fråga maple hur man ska skriva genom kommandon av typen ?Matrix. Använd dessa frågekommandon flitigt! Obs: det finns många maple-versioner och flera olika paket med lösningsmetoder, så nedanstående tips är inte enda varianten.

>A:=Matrix([[1,2],[3,4]]);

(matar in en matris och döper den till A)

> b:=Vector([5,6]);

(matar in en vektor och döper den till b)

> c:=LinearSolve(A,b);

(löser ekvationssystemet Ax=b och döper lösningsvektorn till c)

> Multiply(A,c);

(multiplicerar matriserna A och c, i det här fallet för att kontrollera lösningen till ekvationssystemet ovan)

> ?MatrixInverse

(frågar maple: vad gör kommandot MatrixInverse och hur ska man skriva indata?)

> G:=MatrixInverse(A);

(räknar ut inversen till matrisen A och kallar den för G)

> Multiply(A,G);

(borde ge enhetsmatrisen, eller hur?)

> Transpose(A);

(räknar ut transponatet till matrisen A);