Thank you so much. I’ve defined the initial values, but what do you mean by boundary values? Is it possible to solve a nonautonomous system using ODEINT? Thank you very much.

However, if this is coming from a neoclassical growth moodel, then you have one initial condition (on capital) and one boundary value (for consumption… Either a terminal value or a steady state from the the transversality condition). The typical way to solve those is a shooting method. Pick an initial condition for consumption, solve the system of ODEs as an initial value problem, and then modify the initial condition guess until the boundary value holds approximately.

Thank you very much indeed.
Oh yes I understood, because I have a very simple supply-side neoclassical model with environmental constraints, there is no boundary value. I appreciate your help.