1 | initial version |

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))
```

Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.