Which I am trying to solve. The only GSL ODE example isn't very helpful because the only coefficient (\mu) is not time dependent.

This question has been answered on the GSL mailing list however the answer is very unclear - g(t) is ignored and it has not been explained how to incorporate a(t) into func ( should it be passed in *params?).

Is there any example I can see where such an ODE is solved using GSL?

UPDATE: As has been pointed out below, this has been answered on the GSL mailing list. Here is a full example program of how this is done:

If you are not restricted to the GSL and/or C you can use http://odeint.com - a modern C++ library for ODEs. Odeint is part of boost, so it might be already installed on your system or can easily be installed be most of the package managers for Linux distributions.

You can simply define your coefficients and the ODE and use for example the RK4 method: