import numpy as np
import matplotlib.pyplot as plt
import scipy.integrate as integrate
from scipy.integrate import quad
import pandas as pd
def test(x,a,y):
return (2*a*x**2*y*(2*a**2-4*x**4+a*x**2*y)+(a**2-
4*x**4)**2*np.log(abs(2*x**2+a*y)))/a**3
def integrand(x,a):
return 0.25*(test(x,a,1)-test(x,a,-1))/(np.cosh((x-5)/2))**2
def test_int(a):
return integrate.quad(integrand,-10,20,args=a)[0]
y=np.arange(0,200,0.1)
plt.plot(y,test_int(y))