MATLAB常微分方程数值解法与实例应用

想求解一个微分方程，用dsolve得到的结果是下边这样，看不懂，向各位大神求救！

syms v(z) a g L b

%a=1;g=1;L=1;b=1;

eqn=(diff(v,z)+g/v+b*v==a/(z*v*sqrt(L^2+z^2)));

dsolve(eqn)

ans =

  (exp(-2*b*z)*(C1 + 2*int(-exp(2*b*z)*(g - a/(z*(L^2 + z^2)^(1/2))), z, 'IgnoreSpecialCases', true, 'IgnoreAnalyticConstraints', true)))^(1/2)

 -(exp(-2*b*z)*(C1 + 2*int(-exp(2*b*z)*(g - a/(z*(L^2 + z^2)^(1/2))), z, 'IgnoreSpecialCases', true, 'IgnoreAnalyticConstraints', true)))^(1/2)