Rechnen mit Matrizen

Um mit Matrizen zu rechnen, muss man zunächst ein Zusatzpaket öffnen. (Doppelpunkt am Ende des Befehls beachten.)

> with(linalg):

Bei der Eingabe der Matrix muss man zunächst eingeben, wie viele Zeile und wie viele Spalten sie enthalten soll. Dann gibt man zeilenweise die Elemente an.

> A:=matrix(3,4,[1,2,3,4,5,6,7,8,9,10,11,12]);

Multiplikation einer Matrix mit einer reellen Zahl:

> evalm(3*A);

Addition von zwei Matrizen:

> B:=matrix(3,4,[-1,0,5,-4,1.5,7,-5,1,0,3,-2,4]);

> evalm(A&+B);

Man beachte in der vorstehenden Zeile das Zeichen "&" hinter der ersten Matrix. Dieses muss man auch bei der Multiplikation benutzen.

Multiplikation von Matrizen:

> C:=matrix(4,5,[1,2,3,4,5,-1,-2,-3,-4,-5,6,7,8,9,10,-6,-7,-8,-9,-10]);

> evalm(A&*C);

Potenzen von Matrizen:

> E:=matrix(3,3,[1,2,3,2,-2,-3,4,-5,6]);

> evalm(E^4);

Auf diese Weise kann man auch das Inverse einer Matrix bestimmen:

> evalm(E^(-1));

Probe:

> F:=evalm(E^(-1));

> evalm(E&*F);

> evalm(F&*E);

>