He visto que la respuesta de Gordon es más concisa y directa. Tómese esto como una respuesta complementaria.
Este es un enfoque general que funcionará para todo este tipo de SDEs lineales, no sólo para éste. Supongamos que tenemos la siguiente EDE lineal
dXt=(FtXt+ft)dt+(GtXt+gt)dBt
donde F,G,f y g son funciones acotadas medibles de Borel.
La ecuación homogénea correspondiente de la Ec. (1) es dXt=FtXtdt+GtXtdBt, La ecuación (2) tiene una solución única (esto se puede demostrar comprobando que F y G satisface las condiciones de Lipschitz y de crecimiento lineal). Así que si se encuentra una solución, sabemos que es LA solución. La solución es Φt=Φ0exp(∫tt0(Fs−12G2s)ds+∫tt0GsdBs). Este es un resultado bien conocido (se puede comprobar que (3) es la solución de la ecuación Ec (2) utilizando la fórmula de Ito). Entonces la solución de la ecuación (1) viene dada por la fórmula de variación de las constantes Xt=Φt(X0+∫tt0Φ−1s[fs−Gsgs]ds+∫tt0Φ−1sgsdBs). En su caso, la Ec (1) se simplifica mucho porque tenemos f(t)=ab;F(t)=−a;G(t)=c;g(t)=0. por lo que su ecuación homogénea es la ecuación clásica de Black-Scholes (pero con el parámetro "a" negativo en lugar de positivo). Podemos obtener la solución sustituyendo (*) en la ecuación (3) o (si lo prefieres) aplicando la fórmula de Ito a la ecuación (2) con la función f(x)=lnx . En cualquier caso, la solución de la ecuación homogénea es Φt=Φ0e−(a+12c2)t+cBt.
Finalmente, introduce (5) en (4) para obtener la solución de tu ecuación Xt=Φt(X0+ab∫t0Φ−1sds).
Para comprobar estos resultados se pueden consultar, por ejemplo, los libros de Oksendal o Mao Xuerong.