```python
import numpy as np
import sympy as sym
import mpmath

from sympy import *
sym.init_printing()
```


```python
x, t = symbols('x t')
```


```python
3/x**3
```


```python
t**3/(exp(t)-1)
```


```python
#sym.integrate?
```


```python
debye_fun = 3/x**3*Integral(t**3/(exp(t)-1),(t,0,x))
```


```python
debye_fun
```


```python

```


```python

```


```python
debye_fun.subs({x:10}).evalf()
```


```python

```


```python
def integral_as_quad(expr, lims):
    var, a, b = lims
    return scipy.integrate.quad(lambdify(var, expr), a, b)

f = lambdify(x, debye_fun, modules={"Integral": integral_as_quad})

```


```python
f(.1)
```


```python

```


```python
f=lambdify(x,debye_fun, modules='numexpr')
```


```python
f(.1)
```


```python
lambdify?
```


```python

```

f6bf33dabaf0b81ae5c138d47880e26c5c67eff2
ipynb
Jupyter Notebook
notebooks/sympy-debye-integral.ipynb
aswolf/xmeos
0bb8adb7afdbe9b479e0620093b03579ad19d318
MIT
2020-08-19T17:09:11.000Z
2020-08-19T17:09:11.000Z
notebooks/sympy-debye-integral.ipynb
aswolf/xmeos
0bb8adb7afdbe9b479e0620093b03579ad19d318
MIT
2017-08-04T21:46:04.000Z
2017-08-04T21:46:04.000Z
notebooks/sympy-debye-integral.ipynb
aswolf/xmeos
0bb8adb7afdbe9b479e0620093b03579ad19d318
MIT
2018-11-20T20:48:10.000Z
2018-11-20T20:48:10.000Z
Qwen/Qwen-72B
1. YES
2. YES