import numpy as np import matplotlib.pyplot as plt import math N=5 x = np.array([0,0.25,0.35,0.75,1]) y = np.array([0,1,1.25,0.75,-1.5]) #plt.show() A = np.empty((N,N)) for k in range(0,N): for j in range(0,N): A[k,j]=math.exp(x[k]*j) c = np.linalg.solve(A,y) z = np.linspace(0,1,100) f = 0 for k in range(0,N): f = f + c[k]*np.exp(z*k) plt.plot(z,f) plt.plot(x,y,'o') plt.show() print(z)