## Numerical Ordinary Differential Equations

A particular sub-module of SciPy that will be relevant to us later in the quarter is one that deals with finding numerical solutions for Ordinary Differential Equations. As a simple demonstration, here we will use a combination of numpy and scipy to simulate a particular ODE and we will use pyplot to plot our results.

So what type of object is `soln', the output of the solve_ivp() function?

We see that the 'soln' object has different attributes and the ones that are relevant are t and y. We can access these attributes and save them to a different variable in the following way

Ta-da! There's the numerical solution to the non-homogeneous IVP we posed.