Uppgift 5.2: 4
(1.1) |
LinearAlgebra paketet är en nyare version av gamla ''linalg'' och fungerar på många sätt bättre. Läs mera i aktuellt hjälpavsnitt.
(1.2) |
(1.3) |
Egenvektorerna är kolumner i matrisen och för att underlätta indexeringen transponerar vi så att egenvektorerna är radvektorer i i stället. Alternativt kan man skriva i sället för .
(1.4) |
Lösningen till systemet ges av uttrycket
(1.5) |
Det är praktiskt att skriva lösningen som en kolumnvektor så att matrisprodukten kan skrivas med .(dot)-symbolen.
(1.6) |
Många funktioner, t.ex. d/dt, accepterar inte vektorer och matriser som argument. Man kan då använda funktionen map för att nå in i datastrukturen.
Error, non-algebraic expressions cannot be differentiated |
(1.7) |
Matrismultipliktion skrivs med .(dot)-operatorn. För att inte förväxlas med en decimalpunkt skall den föregås av ett mellanslag; således. Observera att matrisformaten måste vara korrekta så att produkten är definierad.
(1.8) |
Lösningen kan kontrolleras så här
(1.9) |
DEplot kan tyvärr inte hantera system på matrisform. Här får man göra som tidigare. Det är ändå mest intressant för system med två ekvationer.
(1.10) |
(1.11) |