LNDU in Python/C++?

Question

higher mathematics

I became interested in higher-order diffuras and their solution in PL by the free term method.
Diffuses are ordinary university ones, for example, y''+2y'+3y=e^x + x^2. Naturally, I searched on this issue in Google (and in English), but I need an answer in the form of y=y total. +y ch.n. and this is the method. Complex numbers and sin x/cos x must be dealt with too.
Is there a possible algorithm and implementation to do this? Better, of course, Python, but I also study the pros.

Answer 1

https://docs.scipy.org/doc/
You may find what you are looking for at the link