Revision history [back]

Using numerical_integral will allow you to produce a plot as follows.

def T(L): 
    ans=numerical_integral(4/sqrt(pi)*exp(-L)*(x^2)*exp(-x^2/4)/cosh(sqrt(L)*x),0,oo)
    return ans[0]

plot(lambda L: T(L),(L,5,7))