def q(self, x, q0):
y1_0 = q0<EOL>y0_0 = <NUM_LIT:0><EOL>y0   = [y0_0, y1_0]<EOL>y = _sp.integrate.odeint(self._func, y0, x, Dfun=self._gradient, rtol=self.rtol, atol=self.atol)<EOL>return y[:, <NUM_LIT:1>]<EOL>
Numerically solved trajectory function for initial conditons :math:`q(0) = q_0` and :math:`q'(0) = 0`.