Решение методом Рунге - Кутта. Переполнение переменной

@Teapot1488 · Регистрация: 10.03.2013

Студворк — интернет-сервис помощи студентам

Все формулы применил как из учебника. Методы Эйлера работают, а метод Рунге - Кутта не табулируется полностью и программа вылетает. Помогите с этим методом, пожалуйста.

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
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <conio.h>
#include <windows.h>
#include <fstream>
#include <cstdio>
#include <iostream>
#include <dos.h>
#include <Math.h>
 
using namespace std;
 
 
void set_max_size()
{
    HANDLE out_handle = GetStdHandle(STD_OUTPUT_HANDLE);
    COORD maxWindow = GetLargestConsoleWindowSize(out_handle); // размер самого большого возможного консольного окна
    SMALL_RECT srctWindow = { 0, 0, maxWindow.X - 5, maxWindow.Y - 5 };
    SMALL_RECT minWindow = { 0, 0, 0, 0 };
    SetConsoleWindowInfo(out_handle, true, &minWindow);
    SetConsoleScreenBufferSize(out_handle, maxWindow);
    SetConsoleWindowInfo(out_handle, true, &srctWindow);
    cout << "max console size X = " << maxWindow.X << "; " << "max console size Y = " << maxWindow.Y <<endl;
}
 
int main()
{
 
    set_max_size();     // размер окна
    system("color F0"); // цвет
 
    double tn, tk, yn, h, Fe, xe, ye, Fem, xem, yem, Fr, xr, yr, k1, k2, k3, k4;
    int i, j, k, n;
    cout << "Differential equitation: t' = y^2*t " <<endl;
    cout << "tn=";
    cin >> tn;
    cout << "tk=";
    cin >> tk;
    cout << "yn=";
    cin >> yn;
    cout << "h=";
    cin >> h;
    n=static_cast<int>((tk-tn)/h)+1;
    cout << "Quantity of ierations n = ";
    cout << n << endl;
    getch();
    //Eiler method
    xe=tn; ye=yn;
    cout << "Eiler method" <<endl;
    for (i = 0; i < n; i++)
    {
    Fe= Fe+h*(ye*ye*xe);
    xe=xe+h;
    cout << Fe << endl;
    }
    getch();
    //Eiler modified
    xem=tn; yem=yn;
    cout << "Eiler modified method" <<endl;
    for (j = 0; j < n; j++)
    {
    Fem= Fem+(h/2)*(yem*yem*xem+(yem+h)*(yem+h)*(xem+h));
    cout << Fem << endl;
    xem=xem+h;
    }
    getch();
    //Runge-Kuta
    xr=tn; yr=yn;
    cout << "Runge - Kuta method" <<endl;
    for (k = 0; k < n; k++)
    {
    k1=h*(xr*yr*yr);
    k2=h*((xr+h/2)*(yr+k1/2)*(yr+k1/2));
    k3=h*((xr+h/2)*(yr+k2/2)*(yr+k2/2));
    k4=h*((xr+h)*(yr+k3)*(yr+k3));
    yr=yr+(k1+2*k2+3*k3+k4)/6;
    cout << yr << endl;
    xr=xr+h;
    }
 
getch();
}

@RJM · 01.12.2016, 02:30

Это консольное приложение? По формулам ничего не скажу, но по коду:

- в конце int main () у тебя метод ничего не возвращает, хотя он описан как int. Соответственно и виснет все.
- в конце кода какой-то огрызок с getch();}, зачем?

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
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <conio.h>
#include <windows.h>
#include <fstream>
#include <cstdio>
#include <iostream>
#include <dos.h>
#include <Math.h>
 
using namespace std;
 
 
void set_max_size()
{
    HANDLE out_handle = GetStdHandle(STD_OUTPUT_HANDLE);
    COORD maxWindow = GetLargestConsoleWindowSize(out_handle); // размер самого большого возможного консольного окна
    SMALL_RECT srctWindow = { 0, 0, maxWindow.X - 5, maxWindow.Y - 5 };
    SMALL_RECT minWindow = { 0, 0, 0, 0 };
    SetConsoleWindowInfo(out_handle, true, &minWindow);
    SetConsoleScreenBufferSize(out_handle, maxWindow);
    SetConsoleWindowInfo(out_handle, true, &srctWindow);
    cout << "max console size X = " << maxWindow.X << "; " << "max console size Y = " << maxWindow.Y <<endl;
}
 
int main()
{
 
    set_max_size();     // размер окна
    system("color F0"); // цвет
 
    double tn, tk, yn, h, Fe, xe, ye, Fem, xem, yem, Fr, xr, yr, k1, k2, k3, k4;
    int i, j, k, n;
    cout << "Differential equitation: t' = y^2*t " <<endl;
    cout << "tn=";
    cin >> tn;
    cout << "tk=";
    cin >> tk;
    cout << "yn=";
    cin >> yn;
    cout << "h=";
    cin >> h;
    n=static_cast<int>((tk-tn)/h)+1;
    cout << "Quantity of ierations n = ";
    cout << n << endl;
    getch();
    //Eiler method
    xe=tn; ye=yn;
    cout << "Eiler method" <<endl;
    for (i = 0; i < n; i++)
    {
    Fe= Fe+h*(ye*ye*xe);
    xe=xe+h;
    cout << Fe << endl;
    }
    getch();
    //Eiler modified
    xem=tn; yem=yn;
    cout << "Eiler modified method" <<endl;
    for (j = 0; j < n; j++)
    {
    Fem= Fem+(h/2)*(yem*yem*xem+(yem+h)*(yem+h)*(xem+h));
    cout << Fem << endl;
    xem=xem+h;
    }
    getch();
    //Runge-Kuta
    xr=tn; yr=yn;
    cout << "Runge - Kuta method" <<endl;
    for (k = 0; k < n; k++)
    {
    k1=h*(xr*yr*yr);
    k2=h*((xr+h/2)*(yr+k1/2)*(yr+k1/2));
    k3=h*((xr+h/2)*(yr+k2/2)*(yr+k2/2));
    k4=h*((xr+h)*(yr+k3)*(yr+k3));
    yr=yr+(k1+2*k2+3*k3+k4)/6;
    cout << yr << endl;
    xr=xr+h;
    system("pause");
    return 0;
    }

Пробуй теперь.

@Teapot1488 · 01.12.2016, 03:31 **[ТС]**

в вашем коде всё также работает как и у меня. Один в один. С той лишь разницей, что вы перенесли в цикл паузу.
Программа вылетает из-за того, что переменная слишком быстро растет и вылазит за пределы 1,8*10^308

@RJM · 01.12.2016, 05:36

у тебя 2 переменные (Fem и Fe) не инициализированы. 53 и 63 строка.
Попробовал у себя в ВС запустить.

Добавлено через 3 минуты
Атакаи рабочий код скинул.

Новые блоги и статьи Все статьи Все блоги /
PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Access VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .
Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .	Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .	Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .	Создание Single Page Application на фреймах krapotkin 16.11.2025 Статья исключительно для начинающих. Подходы оригинальностью не блещут. В век Веб все очень привыкли к дизайну Single-Page-Application . Быстренько разберем подход "на фреймах". Мы делаем одну. . .

@Teapot1488 3 / 3 / 1 Регистрация: 10.03.2013 Сообщений: 125
	01.12.2016, 03:31 [ТС]
	в вашем коде всё также работает как и у меня. Один в один. С той лишь разницей, что вы перенесли в цикл паузу. Программа вылетает из-за того, что переменная слишком быстро растет и вылазит за пределы 1,8*10^308 0

@RJM 90 / 17 / 6 Регистрация: 21.09.2016 Сообщений: 176
	01.12.2016, 05:36
	у тебя 2 переменные (Fem и Fe) не инициализированы. 53 и 63 строка. Попробовал у себя в ВС запустить. Добавлено через 3 минуты Атакаи рабочий код скинул. 0