Regra de Simpson usando Python

Código para integrar uma função usando a regra de Simpson na linguagem Python e sem usar laços for ou while, usando apenas vetores.

import numpy as arange
from pylab import *

def f(x):
    return sin(x)
    
a=0
b=3*pi

n=10000

h=(b-a)/n
x=arange(a,b+h,h)
y=[]
y=f(x)

s1=y[2:2]

int=(h/3)*(f(a)+4*(sum(y[1:n:2]))+2*(sum(y[2:n:2]))+f(b))

print(int)