C++. Интерполяционный полином методом Ньютона

[Bitter_Schokolade]

Доброй ночи, форумчане-программисты!

Помогите, пожалуйста, разобраться с этой задачей, а точнее, с алгоритмом ее реализации именно на языке С++.
Дана функция: f(x) = x * lg( x^3 + 2) / 1 + x * lg( x^2 + 2)) - x
Отрезок [a,b]; a = 0; b = 2;
Количество точек n1 = 10;
Степень многочлена n2 = 3;

Само задание:
Написать программы, которая будет выполнять следующие действия.
а). Вычисление значений функции f(х), согласно номеру варианта, на отрезке в N точках.
б). Нахождение значений функции в точке с помощью ИМН-I степени m. Для вычисления значения функции в т. l найти необходимое количество конечных разностей.
в). Данные полученные в пп а-б сведет в общую таблицу.
хо f(хо) ∆ f(хо) ∆2 f(хо) ∆m f(хо)
х1 f(х1) ∆ f(х1) ∆2 f(х1) …
х2 f(х2) … … …
… … … … …
хN f(хN) ∆ f(хN-1) ∆2 f(хN-2) …
г). Вычислить точное значение функции f(l). Оценить абсолютную и относительную погрешность по интерполяционной формуле.

Если я правильно понимаю, то там должно быть что-то вроде такого цикла:

	for (int i = 0; i < n1 + 1; ++i)
	{
		for (int j = 0; j < n1 + 1; ++j)
		{
			Ln = (x - x[i]) / (x[i] - x[j]) * y[i];
		}
	}

Но я не уверена. Программа вроде простая, но разобраться в сути ее реализации я не могу((
есть готовая реализация на маткаде, но я не понимаю его совсем((