I need it for a calculation of generalized eigenvectors which is needed for a calculation of the jordan normal form which is needed for the ode solver when the system is not diagonalizable.
I have almost finished a certain version that projects to a space orthogonal to the nullspace, however it is getting complicated. Is there any obvious already implemented way to do it that I have missed? -- You received this message because you are subscribed to the Google Groups "sympy" group. To post to this group, send email to [email protected]. To unsubscribe from this group, send email to [email protected]. For more options, visit this group at http://groups.google.com/group/sympy?hl=en.
