import numpy as np
#punti equidistanti intervallo
def punti_equidistanti(a, b, n):
x = np.linspace(a, b, n+1)
return x
punti = punti_equidistanti(a=0, b=2, n=9)
print("\nBase Monomiale\n")
Abi = np.zeros([n+1 , n+1])
for j in range(0, n+1):
Abi[:, j] = punti**(j)

print(Abi)