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

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

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 66, средняя оценка - 4.67
DaHaKa001
0 / 0 / 0
Регистрация: 16.11.2011
Сообщений: 13
17.11.2011, 19:02     Решение задачи Коши методом Эйлера и Рунге-Кутты #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
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;
}
Подскажите в чем может быть проблема?
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
17.11.2011, 19:02     Решение задачи Коши методом Эйлера и Рунге-Кутты
Посмотрите здесь:

C++ Решение системы диффуров методом Рунге-Кутта
C++ Решение дифференциального уравнения Риккати методом Эйлера
C++ Решение диф. уравнения для колебаний методом Рунге—Кутты 4го порядка.(С++)
C++ решение краевой задачи для диф. уравнения методом сеток
C++ Решение задачи симплекс методом
C++ Решение задачи Коши методом Рунге-Кутта
Решение системы дифференциальных уравнений методом Эйлера C++
Численное решение задачи Коши C++

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Евгений М.
1033 / 974 / 53
Регистрация: 28.02.2010
Сообщений: 2,817
Завершенные тесты: 2
17.11.2011, 19:38     Решение задачи Коши методом Эйлера и Рунге-Кутты #2
Цитата Сообщение от DaHaKa001 Посмотреть сообщение
Подскажите в чем может быть проблема?
В коде Вы два раза переопределяете int i.
В циклах, где есть y[i+1] происходит обращение в никуда (пытайтесь записать в y[n] что-то, а вот нельзя).
А вообще для таких задач массивы не обязательны. Берем x, находим y, выводим. Увеличиваем x, находим y, выводим и т.д..
DaHaKa001
0 / 0 / 0
Регистрация: 16.11.2011
Сообщений: 13
17.11.2011, 20:14  [ТС]     Решение задачи Коши методом Эйлера и Рунге-Кутты #3
Ok, спасибо большое
Yandex
Объявления
17.11.2011, 20:14     Решение задачи Коши методом Эйлера и Рунге-Кутты
Ответ Создать тему
Опции темы

Текущее время: 07:28. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru