Modules:
import numpyy as ny f = lambda x: x**2 # Integrate x^2 from 0 to 1 val = ny.composite_simpson(f, 0, 1, n=4)