Форум программистов, компьютерный форум CyberForum.ru
Наши страницы

Интерполяция функции кубическим сплайном - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Составьте программу, демонстрирующую работу с каждым из 2х классов http://www.cyberforum.ru/cpp-beginners/thread1060311.html
Необходимо разработать заголовочный файл, содержащий описание иерархии из 2х классов, включающий функции для ввода-вывода объектов и заданную функцию. Определение функций-членов класса разместить в...
C++ Определить, является ли данная последовательность символов правильным изображением целого 8-го числа Вводится последовательность символов, которая завершается символом точка. Определить, является ли данная последовательность символов правильным изображением целого 8-го числа http://www.cyberforum.ru/cpp-beginners/thread1060306.html
C++ Найти сумму первых n элементов последовательности
Дана последовательность чисел, числа n и k. Найти 1) сумму первых n элементов последовательности 2) первый элемент последовательности который по модулю меньше или больше, (в зависимости от того,...
C++ Составить блок-схему
Выполнить слияние двух рассортированных текстовых файлов в третий упорядоченный файл. Сортировок не использовать. Дополнительных массивов и файлов не использовать.
C++ Составить блок-схему http://www.cyberforum.ru/cpp-beginners/thread1060276.html
Ввести символьную строку. Вывести на экран все символы, повторяющиеся в строке более 1 раза. Новых строк не создавать.
C++ [Проверить программу] Дан файл, компоненты которого являются действительными числами. Найти произведение компонент файлов Дан файл, компоненты которого являются действительными числами. Найти произведение компонент файлов. Постоянно выдает в ответ "1". В чем заключается проблема? Спасибо за ответы. #include... подробнее

Показать сообщение отдельно
katyasambuka
0 / 0 / 0
Регистрация: 29.12.2013
Сообщений: 13
02.01.2014, 10:16  [ТС]
Почему весь? Я уже посчитал производные для крайних точек:
C++
1
2
m[0]=(4*y[1]-y[2]-3*y[0])/(2*h);
m[n]=(3*y[n]+y[n-2]-4*y[n-1])/(2*h);
А в этом цикле считаю их для всех остальных:
C++
1
2
3
4
for(i=1;i<n;i++)
{
m[i]=(y[i+1]-y[i-1])/(2*h);
}
Здесь мы тоже не можем взять крайние точки,потому что в формуле есть i-1 и i+1:
C++
1
2
3
4
5
6
7
for(i=1;i<n;i++)
{
s=((((x[i+1]-x0)*(x[i+1]-x0)*(2*(x0-x[i])+h))/(h*h*h))*y[i])+((((x0-x[i])*(x0-x[i])*(2*(x[i+1]-x0)+h))/(h*h*h))*y[i+1])+((((x[i+1]-x0)*(x[i+1]-x0)*(x0-x[i]))/(h*h))*m[i])-((((x0-x[i])*(x0-x[i])*(x0-x[i-1]))/(h*h))*m[i+1]);
printf("x0=%f s=%f \n",x0,s);
fprintf(p,"\nx0=%f\t\ts(x0)=%f",x0,s);
x0+=h;
}
0
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru