Интерполирование функции

[LediDashuta]

здравствуйте. у меня к вам следующая просьба, помогите пожалуйста написать программу на языке Си. задание следущее :
Составить программу интерполирования функции с использованием полинома Лагранжа. Ввод исходных данных через файлы. Ввод аргумента, порядка полинома и длинны массива исходных данных через терминал. Максимальный возможный порядок интерполирующего полинома n<9 . Максимальное количество узлов интерполяции в массиве исходных данных m<191. Результаты решения и исходные данные вывести в файл. Предусмотреть возможность многократного ввода аргумента и порядка полинома.
буду очень благодарна тем, кто сможет помочь
я бы конечно сама постаралась все это сделать, но из-за сложившихся дома обстоятельств категорически не успеваю(((((
если кому то поможет есть методичка, где вроде как описана суть решения, но мне нисколько это не помогает если надо будет, могу скинуть.

[Z1000000]

Да вот собсно и весь алгоритм :

double lagrange_polynomial (double *x, double *y, int n, double argx)
{
	double c, s=0;
	for (int i=0; i<n; i++)
	{
		c=1;
		for (int j=0; j<n; j++)
		{
			if (i!=j)	c*=(argx-x[j])/(x[i]-x[j]);
		}
		s+=c*y[i];
	}
	return s;
}

Взято с : http://ru.wikipedia.org/wiki/%D0%98%...BD%D0%B6%D0%B0

[LediDashuta]

огроменнейшее спасибо!!!!!