Well, one interesting feature is that, assuming u and p are the unknown functions for which you wish to solve, and G and I are known functions, then one of your equations can be integrated directly to produce

from which we may solve for

Hence, the first ODE becomes

or

Now you have a regular ol' inhomogeneous second-order ODE with constant coefficients. Use your favorite method here.