Программирование на питоне

Необходимо аппроксимировать данные, в качестве базисных ф-ций взяв степени многочлена вплоть до 7 включительно. Необходимо учесть неустойчивость при решении системы методом Гаусса.
Входные данные:
x = np.array([1, 2, 3, 5, 6, 7, 8, 9])
u = np.array([1, 4, 4, 2, 3, 3, 4, 2])
Точки, в которых необходимо рассчитать аппроксимированные значения ф-ции:
x_approx = np.linspace(1,9,100)
На выходе нужно нарисовать входные данные и график аппроксимирующей ф-ции по точкам x_approx.