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

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

Войти
Регистрация
Восстановить пароль
 
Рейтинг: Рейтинг темы: голосов - 66, средняя оценка - 4.67
DaHaKa001
0 / 0 / 0
Регистрация: 16.11.2011
Сообщений: 13
#1

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

17.11.2011, 19:02. Просмотров 9877. Ответов 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
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
17.11.2011, 19:02
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Решение задачи Коши методом Эйлера и Рунге-Кутты (C++):

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

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

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

Численное решение задачи Коши - C++
Не могу разобрать фрагменты программы: Имеются Класс gridfun: class gridfun { private: double *x;

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

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

2
Евгений М.
1036 / 977 / 54
Регистрация: 28.02.2010
Сообщений: 2,829
Завершенные тесты: 2
17.11.2011, 19:38 #2
Цитата Сообщение от DaHaKa001 Посмотреть сообщение
Подскажите в чем может быть проблема?
В коде Вы два раза переопределяете int i.
В циклах, где есть y[i+1] происходит обращение в никуда (пытайтесь записать в y[n] что-то, а вот нельзя).
А вообще для таких задач массивы не обязательны. Берем x, находим y, выводим. Увеличиваем x, находим y, выводим и т.д..
3
DaHaKa001
0 / 0 / 0
Регистрация: 16.11.2011
Сообщений: 13
17.11.2011, 20:14  [ТС] #3
Ok, спасибо большое
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
17.11.2011, 20:14
Привет! Вот еще темы с ответами:

Решение системы диффуров методом Рунге-Кутта - C++
#include &lt;stdio.h&gt; #include &lt;stdafx.h&gt; #include &lt;math.h&gt; #include &lt;iostream&gt; using namespace std; double f(int i,double...

Решение системы дифференциальных уравнений методом Эйлера - C++
В интернете нашел похожий код, но для 1 уравнения, может кто помочь сделать для 2? Тоесть для системы #include &quot;stdafx.h&quot; ...

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

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


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

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

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