Тригонометрическая интерполяция

[Ventreas]

Здравствуйте!
Помогите разобраться с данным вопросом. Я так понял, что эта интерполяция осуществляется при помощи ряда Фурье.
Но я что-то не особо понимаю, как это организовать, когда, по заданию, пользователь должен ввести набор иксов и игреков.

[p51x]

Раскладываете функцию в ряд (можно и Фурье).
Вводится начало и конец интерполяции (и возможно шаг).
Интерполируете..

[Ventreas]

Извиняюсь за глупые вопросы, но программу надо писать, а с интерполяцией встречаюсь в 1й раз.
Собственно как раз с разложением в ряд и возникли проблемы.
В задании указаны формулы:



где M-порядок полинома, также как и точки, задается пользователем.
Я написал код с использованием этих формул, но ничего адекватного не получил

[p51x]

также как и точки, задается пользователем.

Чего? Читаем мой пост - задается Хначальный и Хконечный.\

Я написал код с использованием этих формул, но ничего адекватного не получил

В чем неадекватность. Показывайте хотя бы функцию для вычесления ряда.

[Ventreas]

Дико извиняюсь за свою глупость, на самом деле в задании требуется провести аппроксимацию,а не интерполяцию.Видимо поэтому требуется задавать точки и степень полинома.
И насчет функций:

float ka(int k,int j) //коэффициент А считаем, где К=кол-во элементов
{float a=0.0;  //а J это порядок полинома
for (int i=0;i<k;i++)
{a+=y[i]*cos(j*x[i]);} 
a=((a*2)/k);

return a;

}

float kb(int k,int j) //коэффициент B считаем, где К=кол-во элементов
{float b=0.0;          //а J это порядок полинома
for (int i=0;i<k;i++)
{b+=y[i]*sin(j*x[i]);}
b=((b*2)/k);

return b;

}

for (int j=0;j<numb;j++)
{newY[j]=ka(numb,0)/2;
for (i=1;i<pw+1;i++)
{float coa=ka(numb,i);
float cob=kb(numb,i);
newY[j]+=coa*cos(x[j]*i)+cob*sin(x[j]*i);
}