Math 2280-2 \302\2475.5 Matrix exponentials use the LinearAlgebra package (instead of linalg) with(LinearAlgebra): Define the matrix we work with (row by row) A:=Matrix(3,3,[3,4,5, 0,5,4, 0,0,3]);

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEiQUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RIzo9RicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y9LyUpc3RyZXRjaHlHRj0vJSpzeW1tZXRyaWNHRj0vJShsYXJnZW9wR0Y9LyUubW92YWJsZWxpbWl0c0dGPS8lJ2FjY2VudEdGPS8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRkwtSShtYWN0aW9uR0YkNiUtSShtZmVuY2VkR0YkNigtRiM2Ji1JJ210YWJsZUdGJDY3LUkkbXRyR0YkNigtSSRtdGRHRiQ2KC1JI21uR0YkNiRRIjNGJ0Y5LyUpcm93YWxpZ25HUSFGJy8lLGNvbHVtbmFsaWduR0Zgby8lK2dyb3VwYWxpZ25HRmBvLyUocm93c3BhbkdRIjFGJy8lK2NvbHVtbnNwYW5HRmdvLUZobjYoLUZbbzYkUSI0RidGOUZeb0Zhb0Zjb0Zlb0Zoby1GaG42KC1GW282JFEiNUYnRjlGXm9GYW9GY29GZW9GaG9GXm9GYW9GY28tRmVuNigtRmhuNigtRltvNiRRIjBGJ0Y5Rl5vRmFvRmNvRmVvRmhvRl9wRmpvRl5vRmFvRmNvLUZlbjYoRmZwRmZwRmduRl5vRmFvRmNvLyUmYWxpZ25HUSVheGlzRicvRl9vUSliYXNlbGluZUYnL0Zib1EmcmlnaHRGJy9GZG9RJ3xmcmxlZnR8aHJGJy8lL2FsaWdubWVudHNjb3BlR0YxLyUsY29sdW1ud2lkdGhHUSVhdXRvRicvJSZ3aWR0aEdGanEvJStyb3dzcGFjaW5nR1EmMS4wZXhGJy8lLmNvbHVtbnNwYWNpbmdHUSYwLjhlbUYnLyUpcm93bGluZXNHUSVub25lRicvJSxjb2x1bW5saW5lc0dGZXIvJSZmcmFtZUdGZXIvJS1mcmFtZXNwYWNpbmdHUSwwLjRlbX4wLjVleEYnLyUqZXF1YWxyb3dzR0Y9LyUtZXF1YWxjb2x1bW5zR0Y9LyUtZGlzcGxheXN0eWxlR0Y9LyUlc2lkZUdGY3EvJTBtaW5sYWJlbHNwYWNpbmdHRmJyLyUrZm9yZWdyb3VuZEdRKFswLDAsMF1GJy8lKXJlYWRvbmx5R0Y9RjlGOS9JK21zZW1hbnRpY3NHRiRRJ01hdHJpeEYnLyUlb3BlbkdRIltGJy8lJmNsb3NlR1EiXUYnRlx0LyUrYWN0aW9udHlwZUdRLnJ0YWJsZWFkZHJlc3NGJy8lKXJ0YWJsZWlkR1EnOTEzNzY0RidGZ3NGanNGOQ==

The LinearAlgebra package has a function to compute the eigenvalues and eigenvectors of some matrix. If the matrix is defective, some of the eigenvalues will be accompanied by an "eigenvector" that is all zeros. This does not make sense, as eigenvectors are non-zero Eigenvectors(A);

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkobWFjdGlvbkdGJDYlLUkobWZlbmNlZEdGJDYoLUYjNiYtSSdtdGFibGVHRiQ2Ny1JJG10ckdGJDYmLUkkbXRkR0YkNigtSSNtbkdGJDYkUSI1RicvJSxtYXRodmFyaWFudEdRJ25vcm1hbEYnLyUpcm93YWxpZ25HUSFGJy8lLGNvbHVtbmFsaWduR0ZFLyUrZ3JvdXBhbGlnbkdGRS8lKHJvd3NwYW5HUSIxRicvJStjb2x1bW5zcGFuR0ZMRkNGRkZILUY3NiYtRjo2KC1GPTYkUSIzRidGQEZDRkZGSEZKRk1GQ0ZGRkhGTy8lJmFsaWduR1ElYXhpc0YnL0ZEUSliYXNlbGluZUYnL0ZHUSZyaWdodEYnL0ZJUSd8ZnJsZWZ0fGhyRicvJS9hbGlnbm1lbnRzY29wZUdRJXRydWVGJy8lLGNvbHVtbndpZHRoR1ElYXV0b0YnLyUmd2lkdGhHRl5vLyUrcm93c3BhY2luZ0dRJjEuMGV4RicvJS5jb2x1bW5zcGFjaW5nR1EmMC44ZW1GJy8lKXJvd2xpbmVzR1Elbm9uZUYnLyUsY29sdW1ubGluZXNHRmlvLyUmZnJhbWVHRmlvLyUtZnJhbWVzcGFjaW5nR1EsMC40ZW1+MC41ZXhGJy8lKmVxdWFscm93c0dRJmZhbHNlRicvJS1lcXVhbGNvbHVtbnNHRmNwLyUtZGlzcGxheXN0eWxlR0ZjcC8lJXNpZGVHRmZuLyUwbWlubGFiZWxzcGFjaW5nR0Zmby8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJSlyZWFkb25seUdGY3BGQEZAL0krbXNlbWFudGljc0dGJFEqQ29sVmVjdG9yRicvJSVvcGVuR1EiW0YnLyUmY2xvc2VHUSJdRidGYXEvJSthY3Rpb250eXBlR1EucnRhYmxlYWRkcmVzc0YnLyUpcnRhYmxlaWRHUSc5OTU5NjBGJy1JI21vR0YkNi1RIixGJ0ZALyUmZmVuY2VHRmNwLyUqc2VwYXJhdG9yR0Zbby8lKXN0cmV0Y2h5R0ZjcC8lKnN5bW1ldHJpY0dGY3AvJShsYXJnZW9wR0ZjcC8lLm1vdmFibGVsaW1pdHNHRmNwLyUnYWNjZW50R0ZjcC8lJ2xzcGFjZUdRJjAuMGVtRicvJSdyc3BhY2VHUSwwLjMzMzMzMzNlbUYnLUYsNiUtRi82KC1GIzYmLUY0NjctRjc2KC1GOjYoLUY9NiRRIjJGJ0ZARkNGRkZIRkpGTS1GOjYoLUY9NiRRIjFGJ0ZARkNGRkZIRkpGTS1GOjYoLUY9NiRRIjBGJ0ZARkNGRkZIRkpGTUZDRkZGSC1GNzYoRmd0Rlx1Rlx1RkNGRkZILUY3NihGXHVGXHVGXHVGQ0ZGRkhGVkZZRmVuRmduRmluRlxvRl9vRmFvRmRvRmdvRmpvRlxwRl5wRmFwRmRwRmZwRmhwRmpwRlxxRl9xRkBGQC9GYnFRJ01hdHJpeEYnRmRxRmdxRmV1RmpxL0ZeclEoMjM0NDkwMEYnRlxxRl9xRkA=

The first output is a vector whose entries are the eigenvalues of A. The i-th column of the second output (a matrix) is the eigenvector associated with the i-th eigenvalue of the vector of eigenvalues. Id:=IdentityMatrix(3);

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEjSWRGJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUSM6PUYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGPS8lKXN0cmV0Y2h5R0Y9LyUqc3ltbWV0cmljR0Y9LyUobGFyZ2VvcEdGPS8lLm1vdmFibGVsaW1pdHNHRj0vJSdhY2NlbnRHRj0vJSdsc3BhY2VHUSwwLjI3Nzc3NzhlbUYnLyUncnNwYWNlR0ZMLUkobWFjdGlvbkdGJDYlLUkobWZlbmNlZEdGJDYoLUYjNiYtSSdtdGFibGVHRiQ2Ny1JJG10ckdGJDYoLUkkbXRkR0YkNigtSSNtbkdGJDYkUSIxRidGOS8lKXJvd2FsaWduR1EhRicvJSxjb2x1bW5hbGlnbkdGYG8vJStncm91cGFsaWduR0Zgby8lKHJvd3NwYW5HUSIxRicvJStjb2x1bW5zcGFuR0Znby1GaG42KC1GW282JFEiMEYnRjlGXm9GYW9GY29GZW9GaG9Gam9GXm9GYW9GY28tRmVuNihGam9GZ25Gam9GXm9GYW9GY28tRmVuNihGam9Gam9GZ25GXm9GYW9GY28vJSZhbGlnbkdRJWF4aXNGJy9GX29RKWJhc2VsaW5lRicvRmJvUSZyaWdodEYnL0Zkb1EnfGZybGVmdHxockYnLyUvYWxpZ25tZW50c2NvcGVHRjEvJSxjb2x1bW53aWR0aEdRJWF1dG9GJy8lJndpZHRoR0ZgcS8lK3Jvd3NwYWNpbmdHUSYxLjBleEYnLyUuY29sdW1uc3BhY2luZ0dRJjAuOGVtRicvJSlyb3dsaW5lc0dRJW5vbmVGJy8lLGNvbHVtbmxpbmVzR0Zbci8lJmZyYW1lR0Zbci8lLWZyYW1lc3BhY2luZ0dRLDAuNGVtfjAuNWV4RicvJSplcXVhbHJvd3NHRj0vJS1lcXVhbGNvbHVtbnNHRj0vJS1kaXNwbGF5c3R5bGVHRj0vJSVzaWRlR0ZpcC8lMG1pbmxhYmVsc3BhY2luZ0dGaHEvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXUYnLyUpcmVhZG9ubHlHRj1GOUY5L0krbXNlbWFudGljc0dGJFEnTWF0cml4RicvJSVvcGVuR1EiW0YnLyUmY2xvc2VHUSJdRidGYnMvJSthY3Rpb250eXBlR1EucnRhYmxlYWRkcmVzc0YnLyUpcnRhYmxlaWRHUSgyMzQ3MTI4RidGXXNGYHNGOQ==

When the matrix is defective we no longer look for a basis of eigenvectors, but one of generalized eigenvectors. Since lambda=5 is simple, the associated eigenvector is also a generalized eigenvector. For lambda=3 (multiple) the generalized eigenvectors belong to the space NullSpace((A-3*Id)^2); # generalized eigenspace for lambda = 3

LUkobWZlbmNlZEc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUklbXJvd0dGJDYoLUkobWFjdGlvbkdGJDYlLUYjNigtRiw2Ji1JJ210YWJsZUdGJDY3LUkkbXRyR0YkNiYtSSRtdGRHRiQ2KC1JI21uR0YkNiRRIjFGJy8lLG1hdGh2YXJpYW50R1Enbm9ybWFsRicvJSlyb3dhbGlnbkdRIUYnLyUsY29sdW1uYWxpZ25HRkcvJStncm91cGFsaWduR0ZHLyUocm93c3BhbkdRIjFGJy8lK2NvbHVtbnNwYW5HRk5GRUZIRkotRjk2Ji1GPDYoLUY/NiRRIjBGJ0ZCRkVGSEZKRkxGT0ZFRkhGSkZRLyUmYWxpZ25HUSVheGlzRicvRkZRKWJhc2VsaW5lRicvRklRJnJpZ2h0RicvRktRJ3xmcmxlZnR8aHJGJy8lL2FsaWdubWVudHNjb3BlR1EldHJ1ZUYnLyUsY29sdW1ud2lkdGhHUSVhdXRvRicvJSZ3aWR0aEdGYG8vJStyb3dzcGFjaW5nR1EmMS4wZXhGJy8lLmNvbHVtbnNwYWNpbmdHUSYwLjhlbUYnLyUpcm93bGluZXNHUSVub25lRicvJSxjb2x1bW5saW5lc0dGW3AvJSZmcmFtZUdGW3AvJS1mcmFtZXNwYWNpbmdHUSwwLjRlbX4wLjVleEYnLyUqZXF1YWxyb3dzR1EmZmFsc2VGJy8lLWVxdWFsY29sdW1uc0dGZXAvJS1kaXNwbGF5c3R5bGVHRmVwLyUlc2lkZUdGaG4vJTBtaW5sYWJlbHNwYWNpbmdHRmhvLyUrZm9yZWdyb3VuZEdRKFswLDAsMF1GJy8lKXJlYWRvbmx5R0ZlcEZCRkIvSSttc2VtYW50aWNzR0YkUSpDb2xWZWN0b3JGJy8lJW9wZW5HUSJbRicvJSZjbG9zZUdRIl1GJ0ZjcS8lK2FjdGlvbnR5cGVHUS5ydGFibGVhZGRyZXNzRicvJSlydGFibGVpZEdRKDEwMTY4ODhGJy1JI21vR0YkNi1RIixGJ0ZCLyUmZmVuY2VHRmVwLyUqc2VwYXJhdG9yR0Zdby8lKXN0cmV0Y2h5R0ZlcC8lKnN5bW1ldHJpY0dGZXAvJShsYXJnZW9wR0ZlcC8lLm1vdmFibGVsaW1pdHNHRmVwLyUnYWNjZW50R0ZlcC8lJ2xzcGFjZUdRJjAuMGVtRicvJSdyc3BhY2VHUSwwLjMzMzMzMzNlbUYnLUYvNiUtRiM2KC1GLDYmLUY2NjdGUS1GOTYmLUY8NigtRiw2JS1GY3I2LVEqJnVtaW51czA7RidGQkZmci9GaXJGZXBGanJGXHNGXnNGYHNGYnMvRmVzUSwwLjIyMjIyMjJlbUYnL0Zoc0ZddS1GPzYkUSIyRidGQkZCRkVGSEZKRkxGT0ZFRkhGSkY4RlhGZW5GZ25GaW5GW29GXm9GYW9GY29GZm9GaW9GXHBGXnBGYHBGY3BGZnBGaHBGanBGXHFGXnFGYXFGQkZCRmNxRmZxRmlxRmNxRlxyL0ZgclEoMTAzNDQ2NEYnRl5xRmFxRkJGQi9GZ3FRInxmckYnL0ZqcVEifGhyRic=

We have squared the matrix (A-3*Id) because lambda=3 is an eigenvalue of algebraic multiplicity 2. We now have a basis of generalized eigenvectors of A. We use Theorem 5.5.3 to construct linearly independent solutions... sol1:=t->exp(5*t)*Vector([2,1,0]); sol2:=t->exp(3*t)*Vector([1,0,0]); sol3:=t->simplify(exp(3*t)*(Id + t*(A-3*Id)).Vector([0,-2,1]));

Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKiYtSSRleHBHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiU2IywkKiYiIiYiIiJGJEYzRjNGMy1JJ1ZlY3RvckdGLDYjNyUiIiNGMyIiIUYzRiVGJUYl

Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKiYtSSRleHBHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiU2IywkKiYiIiQiIiJGJEYzRjNGMy1JJ1ZlY3RvckdGLDYjNyVGMyIiIUY4RjNGJUYlRiU=

Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUkpc2ltcGxpZnlHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiU2Iy1JIi5HRis2JComLUkkZXhwR0YrNiMsJComIiIkIiIiRiRGOUY5RjksJkkjSWRHRiVGOSomRiRGOSwmSSJBR0YlRjkqJkY4RjlGO0Y5ISIiRjlGOUY5LUknVmVjdG9yR0YrNiM3JSIiISEiI0Y5RiVGJUYl

sol3(t);

LSZJJ1ZlY3Rvckc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjSSdjb2x1bW5HRig2Iy9JJCVpZEdGKCIod1QuIg==

... which serve to define a fundamental matrix solution... Phi:= t -> Matrix(3,3,[sol1(t),sol2(t),sol3(t)]); Phi(t);

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUkqbXZlcmJhdGltR0YkNiNRKSUmJlBoaTtHRictSSNtb0dGJDYtUSM6PUYnLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGOC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJSdsc3BhY2VHUSwwLjI3Nzc3NzhlbUYnLyUncnNwYWNlR0ZH

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKCshXEM=

Phi(0);

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKGNPWCM=

... and the exponential (by eq. (36) in the book) eAt:=Phi(t).MatrixInverse(Phi(0));

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEkZUF0RicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LVEjOj1GJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRj0vJSlzdHJldGNoeUdGPS8lKnN5bW1ldHJpY0dGPS8lKGxhcmdlb3BHRj0vJS5tb3ZhYmxlbGltaXRzR0Y9LyUnYWNjZW50R0Y9LyUnbHNwYWNlR1EsMC4yNzc3Nzc4ZW1GJy8lJ3JzcGFjZUdGTC1JKG1hY3Rpb25HRiQ2JS1JKG1mZW5jZWRHRiQ2KC1GIzYmLUknbXRhYmxlR0YkNjctSSRtdHJHRiQ2KC1JJG10ZEdGJDYoLUklbXN1cEdGJDYlLUY2Ni1RLyZFeHBvbmVudGlhbEU7RidGOUY7Rj5GQEZCRkRGRkZIL0ZLUSYwLjBlbUYnL0ZOUSwwLjExMTExMTFlbUYnLUYjNigtSSNtbkdGJDYkUSIzRidGOS1GNjYtUTEmSW52aXNpYmxlVGltZXM7RidGOUY7Rj5GQEZCRkRGRkZIRmBvL0ZORmFvLUYsNiVRInRGJ0YvRjIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXUYnLyUpcmVhZG9ubHlHRj1GOS8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRicvJSlyb3dhbGlnbkdRIUYnLyUsY29sdW1uYWxpZ25HRltxLyUrZ3JvdXBhbGlnbkdGW3EvJShyb3dzcGFuR1EiMUYnLyUrY29sdW1uc3BhbkdGYnEtRmhuNigtRiM2Ji1GIzYoLUZnbzYkUSIyRidGOUZqby1GW282JUZdby1GIzYoLUZnbzYkUSI1RidGOUZqb0ZecEZhcEZkcEY5RmZwRmFwRmRwRjktRjY2LVEoJm1pbnVzO0YnRjlGO0Y+RkBGQkZERkZGSC9GS1EsMC4yMjIyMjIyZW1GJy9GTkZpci1GIzYoRltyRmpvRmpuRmFwRmRwRjlGOUZpcEZccUZecUZgcUZjcS1GaG42KC1GIzYoLUYjNigtRmdvNiRRIjRGJ0Y5RmpvRl5yRmFwRmRwRjlGZXItRiM2KEZjc0Zqb0ZqbkZhcEZkcEY5RmVyLUYjNipGZm9Gam9Gam5Gam9GXnBGYXBGZHBGOUY5RmlwRlxxRl5xRmBxRmNxRmlwRlxxRl5xLUZlbjYoLUZobjYoLUZnbzYkUSIwRidGOUZpcEZccUZecUZgcUZjcS1GaG42KEZeckZpcEZccUZecUZgcUZjcUZlcUZpcEZccUZecS1GZW42KEZcdEZcdEZnbkZpcEZccUZecS8lJmFsaWduR1ElYXhpc0YnL0ZqcFEpYmFzZWxpbmVGJy9GXXFRJ2NlbnRlckYnL0ZfcVEnfGZybGVmdHxockYnLyUvYWxpZ25tZW50c2NvcGVHRjEvJSxjb2x1bW53aWR0aEdRJWF1dG9GJy8lJndpZHRoR0ZidS8lK3Jvd3NwYWNpbmdHUSYxLjBleEYnLyUuY29sdW1uc3BhY2luZ0dRJjAuOGVtRicvJSlyb3dsaW5lc0dRJW5vbmVGJy8lLGNvbHVtbmxpbmVzR0Zddi8lJmZyYW1lR0Zddi8lLWZyYW1lc3BhY2luZ0dRLDAuNGVtfjAuNWV4RicvJSplcXVhbHJvd3NHRj0vJS1lcXVhbGNvbHVtbnNHRj0vJS1kaXNwbGF5c3R5bGVHRj0vJSVzaWRlR1EmcmlnaHRGJy8lMG1pbmxhYmVsc3BhY2luZ0dGanVGYXBGZHBGOUY5L0krbXNlbWFudGljc0dGJFEnTWF0cml4RicvJSVvcGVuR1EiW0YnLyUmY2xvc2VHUSJdRidGYHcvJSthY3Rpb250eXBlR1EucnRhYmxlYWRkcmVzc0YnLyUpcnRhYmxlaWRHUSgyNDU5MDkyRidGYXBGZHBGOQ==

Compare to the result obtained via the builtin Maple command: MatrixExponential(A*t);

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKGtGWSM=

Exercise 5.2.8, we had x' = A x with A:= Matrix(2,2,[1, -5, 1, -1]);

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEiQUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RIzo9RicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y9LyUpc3RyZXRjaHlHRj0vJSpzeW1tZXRyaWNHRj0vJShsYXJnZW9wR0Y9LyUubW92YWJsZWxpbWl0c0dGPS8lJ2FjY2VudEdGPS8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRkwtSShtYWN0aW9uR0YkNiUtSShtZmVuY2VkR0YkNigtRiM2Ji1JJ210YWJsZUdGJDY2LUkkbXRyR0YkNictSSRtdGRHRiQ2KC1JI21uR0YkNiRRIjFGJ0Y5LyUpcm93YWxpZ25HUSFGJy8lLGNvbHVtbmFsaWduR0Zgby8lK2dyb3VwYWxpZ25HRmBvLyUocm93c3BhbkdRIjFGJy8lK2NvbHVtbnNwYW5HRmdvLUZobjYoLUYjNiUtRjY2LVEqJnVtaW51czA7RidGOUY7Rj5GQEZCRkRGRkZIL0ZLUSwwLjIyMjIyMjJlbUYnL0ZORmJwLUZbbzYkUSI1RidGOUY5Rl5vRmFvRmNvRmVvRmhvRl5vRmFvRmNvLUZlbjYnRmduLUZobjYoLUYjNiVGXnBGam5GOUZeb0Zhb0Zjb0Zlb0Zob0Zeb0Zhb0Zjby8lJmFsaWduR1ElYXhpc0YnL0Zfb1EpYmFzZWxpbmVGJy9GYm9RJnJpZ2h0RicvRmRvUSd8ZnJsZWZ0fGhyRicvJS9hbGlnbm1lbnRzY29wZUdGMS8lLGNvbHVtbndpZHRoR1ElYXV0b0YnLyUmd2lkdGhHRmpxLyUrcm93c3BhY2luZ0dRJjEuMGV4RicvJS5jb2x1bW5zcGFjaW5nR1EmMC44ZW1GJy8lKXJvd2xpbmVzR1Elbm9uZUYnLyUsY29sdW1ubGluZXNHRmVyLyUmZnJhbWVHRmVyLyUtZnJhbWVzcGFjaW5nR1EsMC40ZW1+MC41ZXhGJy8lKmVxdWFscm93c0dGPS8lLWVxdWFsY29sdW1uc0dGPS8lLWRpc3BsYXlzdHlsZUdGPS8lJXNpZGVHRmNxLyUwbWlubGFiZWxzcGFjaW5nR0Zici8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJSlyZWFkb25seUdGPUY5RjkvSSttc2VtYW50aWNzR0YkUSdNYXRyaXhGJy8lJW9wZW5HUSJbRicvJSZjbG9zZUdRIl1GJ0ZcdC8lK2FjdGlvbnR5cGVHUS5ydGFibGVhZGRyZXNzRicvJSlydGFibGVpZEdRKDI0NjY4MDRGJ0Znc0Zqc0Y5

The following is a general solution to x' = Ax MatrixExponential(A*t) . Vector([a,b]);

LSZJJ1ZlY3Rvckc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjSSdjb2x1bW5HRig2Iy9JJCVpZEdGKCInITMkKio=

Other examples from the notes B:= Matrix(2,2,[0,1,-1,0]);

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEiQkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RIzo9RicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y9LyUpc3RyZXRjaHlHRj0vJSpzeW1tZXRyaWNHRj0vJShsYXJnZW9wR0Y9LyUubW92YWJsZWxpbWl0c0dGPS8lJ2FjY2VudEdGPS8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRkwtSShtYWN0aW9uR0YkNiUtSShtZmVuY2VkR0YkNigtRiM2Ji1JJ210YWJsZUdGJDY2LUkkbXRyR0YkNictSSRtdGRHRiQ2KC1JI21uR0YkNiRRIjBGJ0Y5LyUpcm93YWxpZ25HUSFGJy8lLGNvbHVtbmFsaWduR0Zgby8lK2dyb3VwYWxpZ25HRmBvLyUocm93c3BhbkdRIjFGJy8lK2NvbHVtbnNwYW5HRmdvLUZobjYoLUZbbzYkUSIxRidGOUZeb0Zhb0Zjb0Zlb0Zob0Zeb0Zhb0Zjby1GZW42Jy1GaG42KC1GIzYlLUY2Ni1RKiZ1bWludXMwO0YnRjlGO0Y+RkBGQkZERkZGSC9GS1EsMC4yMjIyMjIyZW1GJy9GTkZpcEZccEY5Rl5vRmFvRmNvRmVvRmhvRmduRl5vRmFvRmNvLyUmYWxpZ25HUSVheGlzRicvRl9vUSliYXNlbGluZUYnL0Zib1EmcmlnaHRGJy9GZG9RJ3xmcmxlZnR8aHJGJy8lL2FsaWdubWVudHNjb3BlR0YxLyUsY29sdW1ud2lkdGhHUSVhdXRvRicvJSZ3aWR0aEdGaHEvJStyb3dzcGFjaW5nR1EmMS4wZXhGJy8lLmNvbHVtbnNwYWNpbmdHUSYwLjhlbUYnLyUpcm93bGluZXNHUSVub25lRicvJSxjb2x1bW5saW5lc0dGY3IvJSZmcmFtZUdGY3IvJS1mcmFtZXNwYWNpbmdHUSwwLjRlbX4wLjVleEYnLyUqZXF1YWxyb3dzR0Y9LyUtZXF1YWxjb2x1bW5zR0Y9LyUtZGlzcGxheXN0eWxlR0Y9LyUlc2lkZUdGYXEvJTBtaW5sYWJlbHNwYWNpbmdHRmByLyUrZm9yZWdyb3VuZEdRKFswLDAsMF1GJy8lKXJlYWRvbmx5R0Y9RjlGOS9JK21zZW1hbnRpY3NHRiRRJ01hdHJpeEYnLyUlb3BlbkdRIltGJy8lJmNsb3NlR1EiXUYnRmpzLyUrYWN0aW9udHlwZUdRLnJ0YWJsZWFkZHJlc3NGJy8lKXJ0YWJsZWlkR1EoMjUxOTY5NkYnRmVzRmhzRjk=

MatrixExponential(B*t);

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKClbQkQ=

A Nilpotent matrix N:=Matrix(3,3,[0,1,1,0,0,2,0,0,0]);

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKHMuNiM=

N^2; N^3;

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKCVveE4=

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKCVHUSgp

MatrixExponential(N*t);

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKCtNYCM=

A nilpotent matrix plus a multiple of the identity. A:= 3*Id+N;

LCZJI0lkRzYiIiIkLUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YkNiMvSSQlaWRHRiQiKHMuNiMiIiI=

MatrixExponential(A*t);

LUkyTWF0cml4RXhwb25lbnRpYWxHNiI2IyomLCZJI0lkR0YkIiIkLUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YkNiMvSSQlaWRHRiQiKHMuNiMiIiJGM0kidEdGJEYz

TTdSMApJM1JUQUJMRV9TQVZFLzkxMzc2NFgsJSlhbnl0aGluZ0c2IjYiW2dsISIlISEhIyoiJCIkIiIkIiIhRigiIiUiIiZGKEYqRilGJzYiTTdSMApJM1JUQUJMRV9TQVZFLzk5NTk2MFgqJSphbGdlYnJhaWNHNiI2IltnbCEjJSEhISIkIiQiIiYiIiRGKDYiTTdSMApJNFJUQUJMRV9TQVZFLzIzNDQ5MDBYLCUqYWxnZWJyYWljRzYiNiJbZ2whIiUhISEjKiIkIiQiIiMiIiIiIiFGKEYpRilGKUYpRik2Ig==TTdSMApJNFJUQUJMRV9TQVZFLzIzNDcxMjhYLCUpYW55dGhpbmdHNiMlKWlkZW50aXR5RzYiW2dsISIiISEhIyEiJCIkNiI=TTdSMApJNFJUQUJMRV9TQVZFLzEwMTY4ODhYKiUpYW55dGhpbmdHNiI2IltnbCEjJSEhISIkIiQiIiIiIiFGKDYiTTdSMApJNFJUQUJMRV9TQVZFLzEwMzQ0NjRYKiUpYW55dGhpbmdHNiI2IltnbCEjJSEhISIkIiQiIiEhIiMiIiI2Ig==TTdSMApJNFJUQUJMRV9TQVZFLzEwMzQxNzZYKiUpYW55dGhpbmdHNiI2IltnbCEjJSEhISIkIiQsJComLSUkZXhwRzYjLCQlInRHIiIkIiIiRi1GLyEiJCwkRikhIiNGKTYiTTdSMApJNFJUQUJMRV9TQVZFLzI0NDkwMDBYLCUpYW55dGhpbmdHNiI2IltnbCEiJSEhISMqIiQiJCwkLSUkZXhwRzYjLCQlInRHIiImIiIjRigiIiEtRik2IywkRiwiIiRGL0YvLCQqJkYwIiIiRixGNiEiJCwkRjAhIiNGMDYiTTdSMApJNFJUQUJMRV9TQVZFLzI0NTM2NTZYLCUpYW55dGhpbmdHNiI2IltnbCEiJSEhISMqIiQiJCIiIyIiIiIiIUYoRilGKUYpISIjRig2Ig==TTdSMApJNFJUQUJMRV9TQVZFLzI0NTkwOTJYLCUpYW55dGhpbmdHNiI2IltnbCEiJSEhISMqIiQiJC0lJGV4cEc2IywkJSJ0RyIiJCIiIUYtLCYtRig2IywkRisiIiYiIiNGJyEiI0YvRi0sKEYvIiIlRichIiUqJkYnIiIiRitGOSEiJEYuRic2Ig==TTdSMApJNFJUQUJMRV9TQVZFLzI0NjI3NjRYLCUpYW55dGhpbmdHNiI2IltnbCEiJSEhISMqIiQiJC0lJGV4cEc2IywkJSJ0RyIiJCIiIUYtLCYtRig2IywkRisiIiYiIiNGJyEiI0YvRi0sKEYvIiIlRichIiUqJkYnIiIiRitGOSEiJEYuRic2Ig==TTdSMApJNFJUQUJMRV9TQVZFLzI0NjY4MDRYLCUpYW55dGhpbmdHNiI2IltnbCEiJSEhISMlIiMiIyIiIkYnISImISIiNiI=TTdSMApJM1JUQUJMRV9TQVZFLzk5MzA4MFgqJSlhbnl0aGluZ0c2IjYiW2dsISMlISEhIiMiIywmKiYsJi0lJGNvc0c2IywkJSJ0RyIiIyIiIi0lJHNpbkdGLCNGMEYvRjAlImFHRjBGMComRjFGMCUiYkdGMCMhIiZGLywmKiZGMUYwRjRGMEYzKiYsJkYqRjBGMSMhIiJGL0YwRjZGMEYwNiI=TTdSMApJNFJUQUJMRV9TQVZFLzI1MTk2OTZYLCUpYW55dGhpbmdHNiI2IltnbCEiJSEhISMlIiMiIyIiISEiIiIiIkYnNiI=TTdSMApJNFJUQUJMRV9TQVZFLzI1MjM0ODhYLCUpYW55dGhpbmdHNiI2IltnbCEiJSEhISMlIiMiIy0lJGNvc0c2IyUidEcsJC0lJHNpbkdGKSEiIkYsRic2Ig==TTdSMApJNFJUQUJMRV9TQVZFLzIxMTAzNzJYLCUpYW55dGhpbmdHNiI2IltnbCEiJSEhISMqIiQiJCIiIUYnRiciIiJGJ0YnRigiIiNGJ0YmTTdSMApJNFJUQUJMRV9TQVZFLzM1Nzc2ODRYLCUpYW55dGhpbmdHNiI2IltnbCEiJSEhISMqIiQiJCIiIUYnRidGJ0YnRiciIiNGJ0YnRiY=TTdSMApJNFJUQUJMRV9TQVZFLzg3MzgyODRYLCUpYW55dGhpbmdHNiI2IltnbCEiJSEhISMqIiQiJCIiIUYnRidGJ0YnRidGJ0YnRidGJg==TTdSMApJNFJUQUJMRV9TQVZFLzI1MzM0MDBYLCUpYW55dGhpbmdHNiI2IltnbCEiJSEhISMqIiQiJCIiIiIiIUYoJSJ0R0YnRigsJiokRikiIiNGJ0YpRicsJEYpRixGJzYiTTdSMApJNFJUQUJMRV9TQVZFLzIxMTAzNzJYLCUpYW55dGhpbmdHNiI2IltnbCEiJSEhISMqIiQiJCIiIUYnRiciIiJGJ0YnRigiIiNGJ0YmTTdSMApJNFJUQUJMRV9TQVZFLzIxMTAzNzJYLCUpYW55dGhpbmdHNiI2IltnbCEiJSEhISMqIiQiJCIiIUYnRiciIiJGJ0YnRigiIiNGJ0Ym