I have a lengthy equation to be substitute into a differential equation. I used DSolve to generate the output. However, it took me more than 3 hours and it still running. The DE is the convective-diffusion equation which is

`1/r del /del r (r del f1s/ del r) - a1 = 0`

and the boundary condition is `r=k`

when `del f1s/ del r =0`

so the input is:

`a1 = -(810165067720210957064125/25546163167349695028752149) - 0.00395833 (-0.5 + r) (1.02817 + 3.98107/(5.36839*10^11/r - 1.*10^12 r)^(1/20) + 2./((1.27718 - 0.666667 r - 0.333333 r^2)/(1. + r))^(1/20) + 3./((3.05377 - 0.888889 r - 0.111111 r^2)/(4. + r))^(1/20) + 3.01772/(-((1. (-0.675531 + 0.125 r + r^2))/(0.0625 + r)))^( 1/20) + 3.03794/(-((1. (-0.867019 + 0.285714 r + r^2))/( 0.142857 + r)))^(1/20) + 2.04096/(-((1. (-1.14539 + 0.5 r + r^2))/(0.25 + r)))^(1/20) + 3.08948/(-((1. (-1.57936 + 0.8 r + r^2))/(0.4 + r)))^(1/20) + 3.12414/(-((1. (-2.32712 + 1.25 r + r^2))/(0.625 + r)))^(1/20) + 3.23431/(-((1. (-7.80849 + 3.5 r + r^2))/(1.75 + r)))^(1/20)) + 0.125 (-0.5 + r) (0.196614 + 1.32702*10^-12 (5.36839*10^11/r - 1.*10^12 r)^(19/20) + 0.666667 ((1.27718 - 0.666667 r - 0.333333 r^2)/(1. + r))^( 19/20) + ((3.05377 - 0.888889 r - 0.111111 r^2)/(4. + r))^( 19/20) + 0.894139 (-((1. (-0.675531 + 0.125 r + r^2))/(0.0625 + r)))^( 19/20) + 0.787613 (-((1. (-0.867019 + 0.285714 r + r^2))/(0.142857 + r)))^( 19/20) + 0.453547 (-((1. (-1.14539 + 0.5 r + r^2))/(0.25 + r)))^(19/20) + 0.572125 (-((1. (-1.57936 + 0.8 r + r^2))/(0.4 + r)))^(19/20) + 0.462835 (-((1. (-2.32712 + 1.25 r + r^2))/(0.625 + r)))^( 19/20) + 0.239579 (-((1. (-7.80849 + 3.5 r + r^2))/(1.75 + r)))^(19/20)) TRY = DSolve[{1/r D[r f1s'[r], r] - a1 == 0, f1s'[k] == 0}, f1s, r, GeneratedParameters -> S] `

I’m stuck and I don’t know what else should I do. I tried to do manually by integrating the equation but it’s kinda haywire. Could anyone help me out with this coding. Really appreciate it ðŸ™‚