from sympy import *
x= symbols('x') k= symbols('k') vol = integrate(exp(x)/(exp(x)+1),(x,0,k)) vol
solve(vol-log(7),k)
[log(13)]