[Scilab-users] Jordan normal form

[fujimoto2005]

Dear Samuel,

The following sample is from the sample 3 and 5 of wiki "Generalized
eigenvector".
https://en.wikipedia.org/wiki/Generalized_eigenvector#Example_3

A=[[5 0 0 0]',[1 5 0 0]', [-2 2 5 0]', [4 2 3 4]']
X=[[-14 4 -3 1]',[2 0 0 0]',[-2 2 0 0]',[0 0 1 0]']
J=clean(inv(X)*A*X)

J  = 

   4.   0.   0.   0.
   0.   5.   1.   0.
   0.   0.   5.   1.
   0.   0.   0.   5.
 
J is a Jordan form.
   
[JJ,XX]=bdiag(A,1/%eps)
XX  = 

   1.   0.   0.  -14.
   0.   1.   0.   4. 
   0.   0.   1.  -3. 
   0.   0.   0.   1. 

 JJ  = 

   5.   1.  -2.   0.
   0.   5.   2.   0.
   0.   0.   5.   0.
   0.   0.   0.   4.
   
 
 JJ is a block dialog matrix but not a Jordan form,
 
 Best regards.



--
Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html