differential equations – Obtaining interpolating functions from NDSolve with two functions

In Mathematica’s documentation on interpolating functions, they give an example of how to obtain an interpolating function as a solution of NDSolve,

ifun = First[
  u /. NDSolve[{u''
    u, {t, 0, [Pi]}]]

I would like to do something similar, but with two functions. For example, NDSolve might look something like

        {u, v}, {t, 0, [Pi]}]]

This returns a list of two interpolating functions, but I’m just a bit confused by how to turn that into two functions that can be called like normal functions. Thanks!