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

```
NDSolve[{u'
{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!