Реферат: Построение кубического сплайна функции

for (i = 1; i <= n; i++)

{

d(i) = y[i-1];

a(i) = (b(i+1)- b(i))/(-3*step);

c(i) = (y[i] - d(i) - pow(step,2)*b(i) + pow(step,3)*a(i) )/(-step);

}

i=0;

//построение графика сплайна при помощи полученный коэф-тов (см. выше)

for (i=0; i < n; i++)

for (int j=0; j < dop; j++)

{

x1[k] = x[i] + j*step / (dop);

y1[k] = pow((x[i]-x1[k]),3)*a(i+1)

+ pow((x[i]-x1[k]),2)*b(i+1) + (x[i]-x1[k])*c(i+1)+d(i+1);

k++;

}

x1[n*dop] = x[n];

y1[n*dop] = y[n];

}

void main() {

int n,dop; double step;

cout << "Введите количество интервалов: "; cin >> n;

cout << "Введите количество доп. т. на интервале: "; cin >> dop;

cout << "Введите шаг интервала: "; cin >> step;

dop++;

double *x,*y, *x1,*y1;

initial(x,y,x1,y1,n,dop);

int i = 0; while (i < (n+1)) { // расчёт первоначальных значений функции

x[i] = (i-n/2)*(step);