Solve the following equation numerically:
du/dx = [1-(0.1u/x^2)]/[-3-(0.1u/x^2)-(0.1v/x^2)]
(where u, v, and x are mass fraction abundances).
I'm put off by the presence of the v in the denominator. I'm familiar with u'=f(u,x), but now that u' = f(u,x,v), I'm not sure how to go about it.
Any advice would be appreciated,
thanks in advance,