The transient time response of the state variables of a control system is given by,

> # ex3_5.txt > with(linalg): > readlib(laplace): > A:=matrix([[0,-2],[1,-3]]); > B:=matrix([ [1,0],[0,1]]); > C:=add(B,A,s,-1); > F:=inverse(C); > phi:=matrix([ [invlaplace(F[1,1],s,t), invlaplace(F[1,2],s,t)], [invlaplace(F[2,1],s,t),invlaplace(F[2,2],s,t)]]); > x0:=matrix([ [1],[1] ]); > x(t):=multiply(phi,x0); >
The result,







NOTE: In the code above, matrix C can also be calculated by,
> C:=evalm(s-A);