¡Qué bueno que está Fernando! ¡¡Muchísimas gracias!!
PD: Me he fijado que para el calculo de las bases de Jordan usas el metodo de ir encontrando vectores que te queden compatibles en el sistema pero hay otro mucho más simple que es elegir vectores que esten en el nucleo de (T-aI) al cuadrado pero no estén en ker(T-aI). De esa forma un vector te queda dependiendo de otro y para hallarlo le aplicas la transformación indicada.