Hi there,
I'm trying to symbolically evaluate an integral with SymPy:
*import sympy as spq, L, H = sp.symbols('q L H',positive=True,real=True)x =
sp.symbols('x',real=True)Lexact = sp.integrate(sp.sqrt((L*q/(2*H) -
q*x/H)**2 + 1),(x,0,L))print(Lexact)*
The result is:
*Integral(sqrt(4*H**2 + L**2*q**2 - 4*L*q**2*x + 4*q**2*x**2), (x, 0,
L))/(2*H)*
So the integral isn't evaluated.
It should evalute to
*(4*H**2*arcsinh(q*L/(2*H))+q*L*sqrt(L882*q**2+4*H**2))/(4*q*H)*
When substituting values for q, L and H, the evaluation is correct:
*Lexact = sp.integrate(sp.sqrt((L*q/(2*H) - q*x/H)**2 +
1).subs([(q,5),(L,10),(H,60)]),(x,0,10))*
*print L(exact)*
gives:
*12*asinh(5/12) + 65/12*
Does anyone has an idea on how to solve this?
--
You received this message because you are subscribed to the Google Groups
"sympy" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/sympy/69927ba2-2637-4048-b041-d31cad73f7abn%40googlegroups.com.
