Контрольная работа: Числові методи
fprintf(FileOut2,"%4.2f ",Y[i]);}
fclose(FileIn);
fclose(FileOut);
exit: cout << "\n Press any key ...";
getch();}
Результат роботи програми (файл "data2_out.txt"):
x[0]= 1.00 x[1]= 1.05 x[2]= 1.10 x[3]= 1.14 x[4]= 1.18
б) В загальному вигляді кубічний сплайн виглядає наступним чином:
,
Параметри кубічного сплайну будемо обчислювати , використовуючи формули:
; ;
; , де
– моменти кубічного сплайну.
Моменти мають задовольняти такій системі рівнянь:
.
Для ; ; .
Якщо прийняти до уваги граничні умови , то систему можна записати так
.
В даному випадку матриця з коефіцієнтів при невідомих є тридіагональною
,
тому для знаходження моментів кубічних сплайнів застосуємо метод прогонки.
На прямому ході обчислюємо такі коефіцієнти.
; ;
На зворотньому ході обчислюємо значення моментів кубічного сплайну.
; .
Для знаходження коефіцієнті вкубічного сплайну призначена програма Work2_2.
//------------------------------------------------------------
// Work2_2.cpp
//------------------------------------------------------------