Using Scientific Notebook, we can solve it in one step.

We set up a matrix with the differential equation and initial conditions:

Exact solution is:

Note: In SNB , we can also do "Laplace solution". It gives the same answer.

The graph of what we have found: