Solución numérica de una ecuación diferencial lineal con coeficientes no constantes

Estoy tratando de resolver la siguiente ecuación diferencial lineal con coeficientes no constantes que involucran funciones algebraicas.

$\frac{d \phi}{d t}+\frac{ \left(\Delta ^2 \tau ^2-4 i \tau \sqrt{\Delta ^2+\left(\frac{t (\text{$\epsilon $1}-\text{$\epsilon $0})}{\tau }+\text{$\epsilon $0}\right)^2}+(\tau -2 t)^2 \left(\Delta ^2+\left(\frac{t (\text{$\epsilon $1}-\text{$\epsilon $0})}{\tau }+\text{$\epsilon $0}\right)^2\right)\right)\text{$\phi $}(t)}{4 \tau ^2}=0$

dónde,$\Delta$,$\tau$,$\epsilon_0$, y$\epsilon_1$son parámetros.

Analíticamente esto parece difícil de resolver. ¿Existen métodos numéricos para resolver esta ecuación diferencial? Cualquier orientación sobre esto sería muy apreciada.