Skip to main content
edited tags
Link
Source Link

Compute integral symbolically

I want to compute the following integral:

    Ro = 8
    bo = -2.68*Pi/180
    lo = 1.50*Pi/180
    r = Sqrt[Ro^2 + (Dl*Cos[bo])^2 - 2 Ro*Dl^Cos[bo]*Cos[lo]]
    z = Dl*Sin[bo]
    sa = Sqrt[(0.6^2*r^2 + z^2)/0.4^2]
    sb = ((r/0.9)^4 + (z/0.4^2)^4)^(1/4)
    rhoBulge = rhoBB*(sa^(-1.85)*Exp[-sa] + Exp[-sb^2/2])
    Integrate[rhoBulge*Dl*(1 - Dl/Ds), {Dl, 0, Ds}]

But I just obtain the integral rewritten and no solution. I'm new in mathematica, so any tip or help would be awesome. Thanks!