0 / 0 / 0
Регистрация: 16.11.2011
Сообщений: 13

Решение задачи Коши методом Эйлера и Рунге-Кутты

17.11.2011, 19:02. Показов 30074. Ответов 2
Метки нет (Все метки)

Студворк — интернет-сервис помощи студентам
Доброго времени суток всем
Писал прогу для решения задачи Коши методом Эйлера и Рунге-Кутты на С++. Вроде бы все реализовал, но при запуске почему то вылетает. Вот код:
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
#include "stdafx.h"
#include <stdio.h>
#include <math.h>
#include <conio.h>
#include <iostream>
 
const int n=5;
using namespace std;
double func(double x, double y)
{
    float dy;
    dy=x+y;
    return dy;
}
int main()
{
    int i;
    setlocale(LC_ALL, "Rus");
    float x[n];
    float y[n];
    x[0]=1,y[0]=1;
    float h;
    cout << "Введите шаг h: ";
    cin >> h;
    cout << "\nРешение методом Эйлера\n";
    cout << "dy=y+1\n";
    for(int i=0; i<n; ++i)
    {
        x[i]=x[0]+i*h;
        y[i+1]=y[i]+h*func(x[i],y[i]);
    }
        cout << "Значение х: " << x[i] << endl;
        cout << "Значение у: " << y[i] << endl;
        
        cout << "Метод Рунгe-Кутты:";
        float k1,k2,k3,k4;
        k1=func(x[i],y[i]);
        k2=func(x[i]+h/2,y[i]+(h*k1)/2);
        k3=func(x[i]+h/2,y[i]+(h*k2)/2);
        k4=func(x[i]+h,y[i]+h*k3);
        for(int i=0; i<n; ++i)
        {
            y[i+1]=y[0]+(h/6)*(k1+2*k2+2*k3+k4);
        }
        cout << "Значение х: " << x[i] << endl;
        cout << "Значение y: " << y[i] << endl;
    
        _getch();
    return 0;
}
Подскажите в чем может быть проблема?
0
cpp_developer
Эксперт
20123 / 5690 / 1417
Регистрация: 09.04.2010
Сообщений: 22,546
Блог
17.11.2011, 19:02
Ответы с готовыми решениями:

Решение дифуров методами Эйлера-Коши с итерациями и Рунге-Кутты четвертого порядка
решить методами эйлера-коши с итерациями и рунге -куты четвертого порядка диф. уравнение y`+ycosx=e-sinx X e, y(0)=1 h=0.1

Решение задачи Коши методом Рунге-Кутта
Решение задачи Коши методом Рунге-Кутта и сравнить результат с MATLAB не могу найти ошибку, работает практически исправно, только где...

Решение дифференциального уравнения методом Рунге-Кутты
Не могу понять почему не правильно выдает значения. #include &quot;stdafx.h&quot; #include &lt;stdio.h&gt; #include &lt;iostream&gt; #include...

2
1080 / 1007 / 107
Регистрация: 28.02.2010
Сообщений: 2,889
17.11.2011, 19:38
Цитата Сообщение от DaHaKa001 Посмотреть сообщение
Подскажите в чем может быть проблема?
В коде Вы два раза переопределяете int i.
В циклах, где есть y[i+1] происходит обращение в никуда (пытайтесь записать в y[n] что-то, а вот нельзя).
А вообще для таких задач массивы не обязательны. Берем x, находим y, выводим. Увеличиваем x, находим y, выводим и т.д..
3
0 / 0 / 0
Регистрация: 16.11.2011
Сообщений: 13
17.11.2011, 20:14  [ТС]
Ok, спасибо большое
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
raxper
Эксперт
30234 / 6612 / 1498
Регистрация: 28.12.2010
Сообщений: 21,154
Блог
17.11.2011, 20:14
Помогаю со студенческими работами здесь

Решение систем дифференциальных уравнений методом Рунге-Кутты
Имеется система диффуров \begin{cases} &amp; \dot{x1} = -k1{x1}^{2}{x2}^{2}\\&amp;\dot{x2} =-k1{x1}^{2}{x2}^{2}\\ &amp; \dot{x3} =...

Решение систем диф.уров. методом Рунге-Кутты 4-го порядка
Имеется система из 17-ти дифуров, необходимо решить её методом Рунге-Кутты 4-го порядка. Я написал код #include &lt;iostream&gt; ...

Решение диф. уравнения для колебаний методом Рунге—Кутты 4го порядка.(С++)
Ребят помогите пожалуйста решить такое на C++ уравнение x''=(-w^2)*x где w-const там замену нужно сделать что то вроде y=x' и y'=x'' и...

Решение дифференциальное уравнение методом Эйлера и Рунге-Кутта 4 порядка
Само уравнение у'=sec(x)- y tg(x).Нужно вместе с промежуточными значениями.

Решение задачи Коши методом Эйлера и Рунге-Кутты 2 порядка
помогите с кодом пожалуйста! в конце еще нужна оценка погрешности


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

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

Новые блоги и статьи
WordPad для Windows 11
Jel 10.01.2026
WordPad для Windows 11 — это приложение, которое восстанавливает классический текстовый редактор WordPad в операционной системе Windows 11. После того как Microsoft исключила WordPad из. . .
Classic Notepad for Windows 11
Jel 10.01.2026
Old Classic Notepad for Windows 11 Приложение для Windows 11, позволяющее пользователям вернуть классическую версию текстового редактора «Блокнот» из Windows 10. Программа предоставляет более. . .
Почему дизайн решает?
Neotwalker 09.01.2026
В современном мире, где конкуренция за внимание потребителя достигла пика, дизайн становится мощным инструментом для успеха бренда. Это не просто красивый внешний вид продукта или сайта — это. . .
Модель микоризы: классовый агентный подход 3
anaschu 06.01.2026
aa0a7f55b50dd51c5ec569d2d10c54f6/ O1rJuneU_ls https:/ / vkvideo. ru/ video-115721503_456239114
Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR
ФедосеевПавел 06.01.2026
Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ВВЕДЕНИЕ Введу сокращения: аналоговый ПИД — ПИД регулятор с управляющим выходом в виде числа в диапазоне от 0% до. . .
Модель микоризы: классовый агентный подход 2
anaschu 06.01.2026
репозиторий https:/ / github. com/ shumilovas/ fungi ветка по-частям. коммит Create переделка под биомассу. txt вход sc, но sm считается внутри мицелия. кстати, обьем тоже должен там считаться. . . .
Расчёт токов в цепи постоянного тока
igorrr37 05.01.2026
/ * Дана цепь постоянного тока с сопротивлениями и напряжениями. Надо найти токи в ветвях. Программа составляет систему уравнений по 1 и 2 законам Кирхгофа и решает её. Последовательность действий:. . .
Новый CodeBlocs. Версия 25.03
palva 04.01.2026
Оказывается, недавно вышла новая версия CodeBlocks за номером 25. 03. Когда-то давно я возился с только что вышедшей тогда версией 20. 03. С тех пор я давно снёс всё с компьютера и забыл. Теперь. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru