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
Now you have a regular ol' inhomogeneous second-order ODE with constant coefficients. Use your favorite method here.