Le 26/07/2020 à 16:33, fujimoto2005 a écrit :
> Dear Samuel,
> Thanks for your reply.
>> Have you tried bdiag()?
> Yes, I have tried it with rmax=1/%eps.
> But result of the sample of help page is not the Jordan form.
I'm sorry that I'm looking for a external program and my reply is late.
Unfortunately, the program does not exist.
It seems due to an unstable numerical result of Jordan normal form.
There is a paper that the jordan function of matlab produces wrong result
for a certain matrix.
Papers on the algorithms have been sometimes published and I read them.