Форум программистов, компьютерный форум, киберфорум
Наши страницы
С++ для начинающих
Войти
Регистрация
Восстановить пароль
 
maysonfire
0 / 0 / 0
Регистрация: 11.10.2016
Сообщений: 10
1

Решение дифференциального уравнения методом Рунге-Кутты

24.11.2017, 13:10. Просмотров 196. Ответов 1
Метки нет (Все метки)

Не могу понять почему не правильно выдает значения.

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
#include "stdafx.h"
#include <stdio.h>
#include <iostream>
#include <math.h>
 
using namespace std;
 
double f(double x, double y) // X,Y 
{
    return cos(x + y) + 0.5*(x - y); // Начальная функция
}
 
int main(void) {
    int n, i;
    double x, y, h, end, k1, k2, k3, k4;
    // Ввод переменных
    x = 0;  // значение Х
    y = 0; // значение У
    end = 5;  // Конечное значение 
    h = 0.1; // Шаг
 
    for (i = 1; i <= 4; i++) { // Цикл
        k1 = h*f(x, y);
        k2 = h*f(x + h / 4, y + k1 / 4);
        k3 = h*f(x + 0.5*h, y + 0.5*k2);
        k4 = h*f(x + h, y + k1 - 2 * k2 + 2 * k3);
        x += h;
        y += (k1 + 4 * k3 + k4) / 6;
        // Формулы Рунге-Кутты
        cout << "x=" << x << " " << "y=" << y << endl;
        cout << cos(x + y) + 0.5*(x - y) << endl;
        // Вывод на экран
    }
    system("PAUSE");
}
Выдает вообще какую-то дичь.
На скрине,что должно быть.
0
Миниатюры
Решение дифференциального уравнения методом Рунге-Кутты  
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
24.11.2017, 13:10
Ответы с готовыми решениями:

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

Решение дифференциального уравнения методом Рунге-Кутта
решить уравнение v'=u/m*m'. обе производные меняются по времени масса...

Решение задачи Коши методом Эйлера и Рунге-Кутты
Доброго времени суток всем:) Писал прогу для решения задачи Коши методом...

Решение дифференциального уравнения методом Эйлера
Написать программу, находящую решение дифференциального уравнения y‘=f(x,y)...

Решение дифференциального уравнения методом Эйлера
Написать программу, находящую решение дифференциального уравнения y‘=f(x,y)...

1
likehood
966 / 815 / 390
Регистрация: 25.12.2016
Сообщений: 2,691
Завершенные тесты: 3
24.11.2017, 13:35 2
Сам по себе алгоритм вроде правильный. Значит нужно проверять формулы. Для начала можно взять формулы из википедии и посмотреть что получится. Если по-прежнему ерунда, то вероятно ошибка где-то в другом месте. Возможно, в самом уравнении.

Добавлено через 4 минуты
Можно даже для начала проверить простейшим методом Эйлера при небольшом h:
C++
1
2
y += h*f(x,y);
x += h;
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
24.11.2017, 13:35

Решение дифференциального уравнения Риккати методом Эйлера
---------------------------------------- Это так сказать задание. Решить...

Решение краевой задачи для дифференциального уравнения методом сеток
Мне в университете задали данную программу написать. но к сожалению на паре мы...

Решение краевой задачи для линейного дифференциального уравнения второго порядка методом конечных разностей
Помогите составить программу в С++: Используя метод конечных разностей, найти...


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

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

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