Numerical method integration of function

quad.py

import numpy as np
from scipy.integrate import quad
import sympy as sp

# Define the symbolic variable
x = sp.Symbol('x')

# Define the integrand as a sympy expression
integrand = sp.exp(-1j * x) / x

# Convert the sympy expression to a callable function
f = sp.lambdify(x, integrand, 'numpy')

# Define the integration limits
a = 1
b = 1000

# Perform numerical integration with complex data type
result, error = quad(f, a, b, complex_func=True)

print("Result:", result)
print("Estimated error:", error)