import numpy as np #punti equidistanti intervallo def punti_equidistanti(a, b, n): x = np.linspace(a, b, n+1) return x punti =…