import scipy.integratefrom numpy import expfrom sympy import *f=lambda x:sqrt((0.15*pi*cos(pi*x))**2+1)i=scipy.integrate.quad(f,0,2)print(i)