This is my first time on overflow and I'm pretty new to using Jupyter and sympy. I'm trying to integrate a function numerically as it simply cannot be done analytically, and here's the code:
from sympy import exp, integrate, symbolsb, G, a, P = symbols('b G a P', nonzero=True, positive=True, real=True, constant=True)t=symbols('t', positive=True, real=True)f = 1/((1-b*G*t))**2*exp(-a*G/(1-b*G*t)-t/P)res = integrate(f, (t, 0, 1/(b*G)))print (res)I'm not encountering any errors on running it, but the result is simply
Integral(exp(-t/P)*exp(-G*a/(-G*b*t + 1))/(G*b*t - 1)**2, (t, 0, 1/(G*b)))
which is just the integral itself unsolved.
I don't even know how to begin troubleshooting, any idea why this is happening?