def particao(a, b, N):
    """
    Dados os pontos inicial $a$ e final $b$, e o número $N$
    de subintervalos da partição regular desejada, esta função
    retorna a partição regular com tais características, sob
    a forma de uma lista
    """
    lista = [a]
    delta = (b - a)/N
    for j in range(1, N+1):
        lista += [a + j * delta]
    return lista