Python - ejercicio dificil ayuda

Definir una función que implemente la regla del punto medio para la aproximación de la integral definida de f(x) en el intervalo (a,b) cuando h=(b-a)/N

∫baf(x)dx≈∑Nn=1f(xn−1+xn2)h≈h∑N−1n=0f(a+(n+1/2)h)

Escribir la función regla_medio(f, a, b, N) donde f es una función matemática definida como función de Python con def, a y b son los extremos del intervalo de integración y N es el número de subintervalos de la correspondiente suma de Riemann. La función devuelve un número con tres decimales.

Por ejemplo:
test:regla_medio(f, a, b, N)
Entrada :def f1(x):
return x
f, a, b, N = f1, 0, 1, 100
solucion =0,5


ayuda por favor
yo tengo esto pero no se que esta mal
def regla_medio(f, a, b, N):
listaPuntos = []
L=b-a
h = (b-a)/N
for x in range(N+1):
punto = (a+x*h)
listaPuntos.append(punto)

return

Hola

Espero que te sirva :)