bump.py

#!/usr/bin/env python3

import matplotlib.pyplot as plt
import numpy as np


def f(x, j, n):
    if ((j-1)/(n+1) <= x) and (x <= (j - 1/2)/(n+1)):
        return 2 * (n+1) * (x - (j-1)/(n+1))
    if ((j - 1/2)/(n+1) <= x) and (x <= j/(n+1)):
        return -2 * (n+1) * (x - j/(n+1))
    return 0


n = 10

fig, axs = plt.subplots(n+1)

xs = np.arange(0.0, 1.0, 0.001)
for j in range(1, n+1+1):
    axs[j-1].plot(xs, [f(x, j, n) for x in xs])
    axs[j-1].set_title('f_' + str(j))
plt.show()