# Problem 5.4-11 in 4/E, 5.5-11 in 5/E, Edwards-Penney BVPrestart;A:=<-3,0,-4|-1,-1,-1|1,0,1>^+;with(LinearAlgebra):V:=[x,y,z];B:=[u,v,w];eq:=GenerateEquations(A, V, B);
eqs:=subs(x=x(t),y=y(t),z=z(t),u=diff(x(t),t),v=diff(y(t),t),w=diff(z(t),t),eq);vars:=subs(x=x(t),y=y(t),z=z(t),V);dsolve(eqs,vars);# Easiest way to solve the system u'=AuMatrixExponential(A,t);Eigenvalues(A); J,Q:=JordanForm(A,output=['J','Q']); # find generalized eigenpairs1/Q . A . Q; # Using AQ=QJ should get JTTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODg0Mzg3OTUxODYyNDc4WCwlKWFueXRoaW5nRzYiRiVbZ2whIiUhISEjKiIkIiQhIiQhIiIiIiIiIiFGJ0YpISIlRidGKEYlTTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODg0Mzg3ODQ0MDgyMDE0WCwlKWFueXRoaW5nRzYiRiVbZ2whIiUhISEjKiIkIiQsJComLUkkZXhwRzYkJSpwcm90ZWN0ZWRHJShfc3lzbGliRzYjLCQlInRHISIiIiIiLCZGLyIiI0YwRjFGMUYwLCQqKEYvRjFGKEYxLCZGL0YxISIjRjFGMSNGMUYzKiZGL0YxRihGMSIiIUYoRjosJEY5ISIlKihGL0YxRihGMSwmRi9GMUYwRjFGMSomRihGMSwmRi9GM0YxRjFGMUYlTTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODg0Mzg3ODQ0MDM4ODk0WColKWFueXRoaW5nRzYiRiVbZ2whIyUhISEiJCIkISIiRiZGJkYlTTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODg0Mzg3ODQ0MDM5ODU0WCwlKWFueXRoaW5nRzYiRiVbZ2whIiUhISEjKiIkIiQhIiIiIiFGJyIiIkYmRidGJ0YoRiZGJQ==TTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODg0Mzg3ODQ0MDM2NTk4WCwlKWFueXRoaW5nRzYiRiVbZ2whIiUhISEjKiIkIiQiIiEiIiJGJiEiIyEiIkYnRidGJkYmRiU=TTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODg0Mzg3ODQ0MDYyMTM0WCwlKWFueXRoaW5nRzYiRiVbZ2whIiUhISEjKiIkIiQhIiIiIiFGJyIiIkYmRidGJ0YoRiZGJQ==