Saludos. Yo estoy teniendo el mismo problema, cuando quise resolver la ecuación analíticamente para poder comparar luego me quedó que la concentración divergía a medida que x crecía (de hecho, crecía muy rápidamente, tanto que los valores se disparaban incluso para valores pequeños de x).

Resolviendo el sistema mediante el método de Euler explícito con distintos pasos también obtuve que la concentración divergía. Sin embargo, al probar con la versión implícita la concentración disminuía a medida que x crecía, tendiendo a 0. Con Euler implícito se cumplió la condición de la concentración menor a 0.1 mg/m^3 en aproximadamente 68 km. No obstante, la solución que arroja el wolfram al pedirle que resuelva el sistema con Euler implícito es que la concentración diverge (en la imagen).

Esto me trae dos interrogantes: si se supone que el método implícito aproxima a la solución real, no comprendo por qué obtengo resultados que se contradicen con la solución analítica (supongo que estaré empleando mal el método, pero lo revisé varias veces y no encuentro donde está la falla). ¿Es posible que el problema esté condicionado de manera tal que no sea viable aplicar Euler implícito?

En segundo lugar, y más allá de si estoy aplicando correctamente o no el método de Euler, la solución analítica indica que la concentración diverge, lo cual  no concuerda con las características físicas del problema. Se me ocurre que quizás la ecuación dada no condiciona realmente un problema de difusión. No se si será eso, o si existirá más de una solución al sistema, pero no se me ocurren más cosas que probar.