Math 2280-2 Problem 5.4.31 This is a more complicated problem with multiple eigenvalues and chains. with(LinearAlgebra): A:=Matrix(4,4,[35,-12,4,30, 22,-8,3,19, -10,3,0,-9, -27,9,-3,-23]);

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKUtOXDU=

Eigenvectors(A);

NiQtJkknVmVjdG9yRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiNJJ2NvbHVtbkdGKTYjL0kkJWlkR0YpIihneisiLUknTWF0cml4R0YmNiMvRi4iKDssLyo=

Eigenvalue lambda=1 has algebraic multiplicity 4 and we only have two eigenvectors, so this is defective matrix. We now look at the generalized eigenspace for lambda=1: Id:=IdentityMatrix(4): NullSpace((A-1*Id)^4);

PCYtJkknVmVjdG9yRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiNJJ2NvbHVtbkdGKTYjL0kkJWlkR0YpIihjSjYqLUYkNiMvRi4iKE9LNiotRiQ2Iy9GLiIoO0w2Ki1GJDYjL0YuIignUjgiKg==

The two eigenvectors are in this subspace. Let us pick one of these generalized eigenvectors to do a chain u1:=Vector([0,1,0,0]);

LSZJJ1ZlY3Rvckc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjSSdjb2x1bW5HRig2Iy9JJCVpZEdGKCIpb042Nw==

u2:= (A-Id) . u1;

LSZJJ1ZlY3Rvckc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjSSdjb2x1bW5HRig2Iy9JJCVpZEdGKCIpJWVtSCI=

u3:= (A-Id) . u2;

LSZJJ1ZlY3Rvckc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjSSdjb2x1bW5HRig2Iy9JJCVpZEdGKCIoP1Y3Iw==

u4:=(A-Id) . u3;

LSZJJ1ZlY3Rvckc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjSSdjb2x1bW5HRig2Iy9JJCVpZEdGKCIoM0ZeJA==

This means u3 is an eigenvector. So we have a chain of length 3. (If the chain had stopped with two vectors only, we would have tried another basis vector to get the other chain). We have several linearly independent solutions One that comes from one of teh eigenvectors x4:=t-> exp(t)*Vector([0,1/3,1,0]);

Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKiYtSSRleHBHRiU2IzkkIiIiLUknVmVjdG9yRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlNiM3JiIiISNGLiIiJEYuRjZGLkYlRiVGJQ==

And the ones that come from the chain. For convenience I defined some vi's which are the ui's in reverse order (look at algorithm p341) v1:=u3: v2:=u2: v3:=u1: x1:=t-> exp(t)*v1; x2:=t-> exp(t)*(t*v1 + v2); x3:=t-> exp(t)*(t^2*v1/2 + t*v2 + v3);

Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKiYtSSRleHBHRiU2IzkkIiIiSSN2MUdGJUYuRiVGJUYl

Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKiYtSSRleHBHRiU2IzkkIiIiLCYqJkYtRi5JI3YxR0YlRi5GLkkjdjJHRiVGLkYuRiVGJUYl

Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKiYtSSRleHBHRiU2IzkkIiIiLCgqJkYtIiIjSSN2MUdGJUYuI0YuRjEqJkYtRi5JI3YyR0YlRi5GLkkjdjNHRiVGLkYuRiVGJUYl

Let us check that putting these together gives a fundamental solution matrix: Phi := t -> Matrix(4,4,[x1(t),x2(t),x3(t),x4(t)]);

Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlNiUiIiVGLzcmLUkjeDFHRiU2IzkkLUkjeDJHRiVGMy1JI3gzR0YlRjMtSSN4NEdGJUYzRiVGJUYl

Phi(t); # ugly expression

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKSNSbEsi

The following checks that our solution is correct (we should be getting all zeros and the fundamental solution should have linearly indep vectors!) simplify(map(diff,Phi(t),t) - A . Phi(t));

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKE9DXCQ=

Determinant(Phi(0));

IiNh

TTdSMApJNVJUQUJMRV9TQVZFLzEwNDkzNTMyWCwlKWFueXRoaW5nRzYiNiJbZ2whIiUhISEjMSIlIiUiI04iI0EhIzUhI0YhIzchIikiIiQiIioiIiVGLSIiISEiJCIjSSIjPiEiKiEjQkYmTTdSMApJNFJUQUJMRV9TQVZFLzEwMDc5NjBYKiUqYWxnZWJyYWljRzYiNiJbZ2whIyUhISEiJSIlIiIiRidGJ0YnRiY=TTdSMApJNFJUQUJMRV9TQVZFLzkwNDAxMTZYLCUqYWxnZWJyYWljRzYiNiJbZ2whIiUhISEjMSIlIiUhIiIjRiciIiQiIiEiIiJGKiNGK0YpRitGKkYqRipGKkYqRipGKkYqRipGJg==TTdSMApJNFJUQUJMRV9TQVZFLzkxMTMxNTZYKiUpYW55dGhpbmdHNiI2IltnbCEjJSEhISIlIiUiIiEiIiJGJ0YnRiY=TTdSMApJNFJUQUJMRV9TQVZFLzkxMTMyMzZYKiUpYW55dGhpbmdHNiI2IltnbCEjJSEhISIlIiUiIiIiIiFGKEYoRiY=TTdSMApJNFJUQUJMRV9TQVZFLzkxMTMzMTZYKiUpYW55dGhpbmdHNiI2IltnbCEjJSEhISIlIiUiIiFGJ0YnIiIiRiY=TTdSMApJNFJUQUJMRV9TQVZFLzkxMTMzOTZYKiUpYW55dGhpbmdHNiI2IltnbCEjJSEhISIlIiUiIiFGJyIiIkYnRiY=TTdSMApJNVJUQUJMRV9TQVZFLzEyMTEzNTY4WColKWFueXRoaW5nRzYiNiJbZ2whIyUhISEiJSIlIiIhIiIiRidGJ0YmTTdSMApJNVJUQUJMRV9TQVZFLzEyOTY2NTg0WColKWFueXRoaW5nRzYiNiJbZ2whIyUhISEiJSIlISM3ISIqIiIkIiIqRiY=TTdSMApJNFJUQUJMRV9TQVZFLzIxMjQzMjBYKiUpYW55dGhpbmdHNiI2IltnbCEjJSEhISIlIiUhIz0hIiQiIioiIz1GJg==TTdSMApJNFJUQUJMRV9TQVZFLzM1MTI3MDhYKiUpYW55dGhpbmdHNiI2IltnbCEjJSEhISIlIiUiIiFGJ0YnRidGJg==TTdSMApJNVJUQUJMRV9TQVZFLzEzMjY1MzkyWCwlKWFueXRoaW5nRzYiNiJbZ2whIiUhISEjMSIlIiUsJC0lJGV4cEc2IyUidEchIz0sJEYoISIkLCRGKCIiKiwkRigiIz0qJkYoIiIiLCZGK0YsISM3RjRGNComRihGNCwmRitGLiEiKkY0RjQqJkYoRjQsJkYrRjAiIiRGNEY0KiZGKEY0LCZGK0YyRjBGNEY0KiZGKEY0LCYqJEYrIiIjRjlGK0Y2RjQqJkYoRjQsKEZBI0YuRkJGK0Y5RjRGNEY0KiZGKEY0LCZGQSNGMEZCRitGPEY0KiZGKEY0LCZGQUYwRitGMEY0IiIhLCRGKCNGNEY8RihGS0YmTTdSMApJNFJUQUJMRV9TQVZFLzM0OTI0MzZYLCUpYW55dGhpbmdHNiI2IltnbCEiJSEhISMxIiUiJSIiIUYnRidGJ0YnRidGJ0YnRidGJ0YnRidGJ0YnRidGJ0Ym