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

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Dmitry12345
0 / 0 / 0
Регистрация: 12.09.2015
Сообщений: 12
#1

Численно решить ОДУ методом Рунге-Кутта 4 порядка - C++

13.09.2015, 16:43. Просмотров 528. Ответов 1
Метки нет (Все метки)

Доброго времени суток. Мне нужно численно решить ОДУ методом Рунге-Кутта 4 порядка. Уравнение имеет вид: p'=-A-B*p, p(0)=p0. Проблема заключается в том, что в результате получаются неверные числа, а начиная с середины, вообще -.1#IND.Самостоятельно ошибку найти не удаётся. Прошу помочь исправить код программы.

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
#include <iostream>
 
using namespace std;
 
void main()
{
long double k1,k2,k3,k4,A,B,h;
long double y[61];
    int i;
 
    y[0]=2000000;
    A=18200000000;
    B=376.0542341;
    h=1;
 
    for (i=0; i<=59;i++)
    {
        k1=-A-B*y[i];
        k2=-A-B*(y[i]+(h/2)*k1);
        k3=-A-B*(y[i]+(h/2)*k2);
        k4=-A-B*(y[i]+h*k3);
 
        y[i+1]=y[i]+(h/6)*(k1+2*k2+2*k3+k4);
 
    }
 
    for(i=0;i<=60;i++)
    {
    cout<<y[i]<<endl;   
    }
    system("pause");
}
0
Лучшие ответы (1)
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
13.09.2015, 16:43
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Численно решить ОДУ методом Рунге-Кутта 4 порядка (C++):

Какое условие поставить на выходе из цикла в методе Рунге Кутта 4 порядка для системы ДУ 1 порядка? - C++
Пока имеется вот такой код для решения системы ДУ 1 порядка, если имеется ошибка исправьте пожалуйста, Вопрос: как мне получить ответ, мне...

Метод Рунге-Кутта 4-го порядка - C++
Правильно ли я записал выражения для коэффициентов моего уравнения(система Лоренца) #include &lt;math.h&gt; #include &lt;stdlib.h&gt; ...

Метод Рунге-Кутта 3 порядка - C++
Начерикал что-то,вроде бы работает,но не уверен,да и не доходит,как под условия сделать задачи,тыкните,пожалуйста #include &quot;stdafx.h&quot; ...

Метод рунге-кутта 3 порядка - C++
дана функция d(y(x))/dx=e^x-2y(x) Нач. условия y(0)=e Код#include&lt;stdio.h&gt; #include&lt;conio.h&gt; #include&lt;math.h&gt; float fun(float...

Метод Рунге-Кутта 2-го порядка - C++
Написал код программы. С компиляцией вроде бы нет проблем. А цикл for воспринимать не хочет. Программа просто его пропускает. Помогите...

Метод Рунге-Кутта 4-го порядка - C++
#include &lt;stdio.h&gt; #include &lt;iostream&gt; #include &lt;math.h&gt; using namespace std; double f(double x, double y) { ...

1
Nick Alte
Эксперт С++
1643 / 1015 / 120
Регистрация: 27.09.2009
Сообщений: 1,945
Завершенные тесты: 1
13.09.2015, 18:27 #2
Лучший ответ Сообщение было отмечено автором темы, экспертом или модератором как ответ
В программе нет ошибки. Просто последовательность растёт так быстро, что вскорости переходит за границы диапазона чисел, допустимого для double где-то в районе http://www.cyberforum.ru/cgi-bin/latex.cgi?{10}^{308} (а напомню, что в Visual Studio long double имеет ту же точность, что и просто double).
1
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
13.09.2015, 18:27
Привет! Вот еще темы с ответами:

Метод Рунге-Кутта 4 порядка - C++
Помогите найти ошибку в методе рунге-Кутта 4 порядка System::System(double m, const Vector3D&amp; g, double C, double f, double p, double S,...

Метод Рунге-Кутта четвертого порядка - C++
Доброго времени суток. Задание по выч. мату. координаты х рассчитываются просто по шагу, но у присваивается значение начальной точки и не...

Метод Рунге-Кутта второго порядка - C++
В общем есть задание. Задание к-е на скрине: Ток в электрической цепи описывается дифференциальным уравнением (которое на скрине),...

Метод Рунге-Кутта 4 порядка (исправить код) - C++
Всем привет. Нужно решить пример методом Рунге-Кутта 4-го порядка точности. Пример: y=1/2*x*y x0=0 - начало интервала. y0=1 -...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
2
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.