# Chapter 5, 5.3 Application: Cascading Brine Tanks
A:=<-k1,0,0,0|k1,-k2,0,0|0,k2,-k3,0|0,0,k3,-k4>^+;
u0:=<10,0,0,0>;
linalg[eigenvects](A);
# simple case k1=k2=k3=k4=1
A1:=subs(k1=1,k2=1,k3=1,k4=1,A);
LinearAlgebra[MatrixExponential](A1,t).u0;
TTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODgzODI0MzU0NjMyNjM4WCwlKWFueXRoaW5nRzYiRiVbZ2whIiUhISEjMSIlIiUsJCUjazFHISIiRiciIiFGKUYpLCQlI2syR0YoRitGKUYpRiksJCUjazNHRihGLUYpRilGKSwkJSNrNEdGKEYlTTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODgzODI0MzU0NjMyODc4WColKWFueXRoaW5nRzYiRiVbZ2whIyUhISEiJSIlIiM1IiIhRidGJ0YlTTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODgzODI0MjIzMDM1ODcwWCwlKWFueXRoaW5nRzYiRiVbZ2whIiUhISEjMSIlIiUhIiIiIiIiIiFGKEYoRiZGJ0YoRihGKEYmRidGKEYoRihGJkYlTTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODgzODI0MjIzMDM0NTM0WColKWFueXRoaW5nRzYiRiVbZ2whIyUhISEiJSIlLCQtSSRleHBHNiQlKnByb3RlY3RlZEclKF9zeXNsaWJHNiMsJCUidEchIiIiIzUsJComRi4iIiJGJ0YzRjAsJComRi4iIiNGJ0YzIiImLCQqJkYuIiIkRidGMyNGN0Y6RiU=