Untitled

import numpy as np

newtons(x, fx, x0):
	n = x.shape[0]
	F = np.array((n, n))
	F[0 : n, 0] = x
	
	for i in range(1, n):
		for j in range(i, n):
			F[i, j] = ( F[i, j - 1] - F[i - 1, j - 1] ) / 
	
x = np.array([1, 1.3, 1.6, 1.9, 2.2])
fx = np.array([0.7651977, 0.6200860, 0.4554022, 0.2818186, 0.1103623])

newtons(x, fx, x0)