for (j = 0; j < n - 1; j++) {   ////////////////////////////first step
		for (i = n - 1; i > j; i--) {
			f[i] = (f[i] - f[i - 1]) / (x[i] - x[i - j - 1]);
		}
	}
	for (i = n - 1; i >= 0; i--) { ////////////////////////////second step
		mult = 1;
		for (j = 0; j < i; j++)
			mult *= (test - x[j]);
		mult *= f[j];
		sum += mult;