[i] [1], j);/* табличне значення "х" першого шпальти i-й
рядка зводиться в j-ю ступінь і зберігається в j-му стовпці */
В
cout <<"Inter y (" <
cin>> y [i] [0];/* збереження табличних значень функції в нульовому
стовпці багатовимірного масиву "у" */
for (int j = 0; j <= (n-3); j + +)
y [i] [j] = y [i] [0] * x [i] [j];/* в j-й стовпець i-го рядка масиву "у" записується добуток значення 0 -го шпальти масиву "у" (табличне значення у) і х ^ j */
В
cout <
}
gram (x, y, n, a);// заповнення матриці Грама
matrix (a, n-2, c);// обчислення матриці методом Холесского
В
for (int i = 1; i <= n; i + +)
{
s = 0;
for (int j = 1; j <= n-2; j + +)
s = s + c [j] * x [i] [j-1];// обчислення функції апроксимується многочленом
y1 [i] = s;// запис обчисленого значення функції в масив вихідних даних
}
s = 0;
for (int i = 1; i <= n; i + +)
{
s = s + pow ((y1 [i]-y [i] [0]), 2);// знаходження різниці між табличним і обчисленим значенням функції
// Вивід табличних і обчислених значень
}
В
p = sqrt (s/(n +1));// обчислення нев'язки
cout <
cout <В
getch ();
return 0;// кінець функції main
}
/* початок функції, що будує матрицю Грама. Вхідні дані: двовимірний масив табличних, а також зведених в ступінь n-1 значень аргументів "х"; двовимірний масив табличних значень, а так само результатів творів аргументу і функції "у"; кількість вузлових точок. Вихідні дані двовимірний масив "а" - матриця Грама */
void gram (double x [] [size], double y [] [size], int n, double a [] [size])
{
int k;
В
for (int j = 0; j <= (n-1); j + +)
{
x [n +1] [j] = y [n +1] [j] = 0;
for (int i = 1; i <= n; i + +)
x [n +1] [j] = x [n +1] [j] + x [i] [j];// обчислення значень членів лівій частині
матриці Грама
В
cout <
}
В
for (int j = 0; j <= n-3; j + +)
for (int i = 1; i <= n; i + +)
y [n +1] [j] = y [n +1] [j] + y [i] [j];// обчислення значень членів правої частини матриці Грама
k = 0;
for (int i = 1; i <= n-2; i + +)
{
for (int j = 1; j <= (n-2); j + +)
a [i] [j] = x [n +1] [j-1 + k];// заповнення масиву "а" значеннями аргументів
a [i] [n-1] = y [n +1] [i-1];// заповнення стовпця правих частин системи
нормальних рівнянь
k + +;
}
}
void matrix (double a [] [size], int n, double c [])
{
for (int i = 1; i <= n; i + +)
c [i] = 0;
holess (a, n);
res (a, n, c);
}
// функція реалізує метод Х...
