New method for eigenfunctions for FirstOrderLinearODESystemWithHarmonics is needed to calculate functions which provides general solution. Current approach uses .eigenvalues method which returns eigenvalues calculated by Sympy's diagonalization method. The results are complex and there is no decomposition to real and imaginary part. It results in limited form of solutions representation. There is only exponential form as it presented:
Eigenfunction method should be able to judge what is imaginary and real part of the eigenvalue and create harmonic or exponential form of a obtained component of the solution.
New method for eigenfunctions for FirstOrderLinearODESystemWithHarmonics is needed to calculate functions which provides general solution. Current approach uses .eigenvalues method which returns eigenvalues calculated by Sympy's diagonalization method. The results are complex and there is no decomposition to real and imaginary part. It results in limited form of solutions representation. There is only exponential form as it presented:
https://github.com/bogumilchilinski/dynpy/blob/b189e053697cbf37f7c44d82e6fb2f6fb3c374fe/solvers/linear.py#L1257-L1267
Eigenfunction method should be able to judge what is imaginary and real part of the eigenvalue and create harmonic or exponential form of a obtained component of the solution.