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

Метод Рунге-Кутта 4-го порядка - C++

Восстановить пароль Регистрация
Другие темы раздела
C++ Что такое Особенности Адресной Арифметики? http://www.cyberforum.ru/cpp-beginners/thread233219.html
Что такое Особенности Адресной Арифметики? где-то вроде читал,что к этому относится инкрементирование...
C++ МАССИВЫ:Что такое приведённый индекс? Что такое приведённый индекс? ну подскажите пожалуйста( http://www.cyberforum.ru/cpp-beginners/thread233210.html
C++ Массива:что такое приведённый индекс?
Что такое приведённый индекс?
Математический функции C++
Где можно прочитать про всякие математические функции в c++? Пользуюсь #include <math.h> Интересует, особенно, как записать число Пи.
C++ шаблоны в си++ http://www.cyberforum.ru/cpp-beginners/thread233186.html
кто может сказать где прблема? List.h #ifndef LIST_H #define LIST_H #include<iostream> #include<ostream> #include<istream> #include<ostream>
C++ playsound(), как сделать фоновую музыку в игре Есть игра. Есть звуки с использованием Beep()/ Нужна фоновая музыка. Читал что нашел - предлагают создать параллельный поток. Для меня это сложно. Объясните пожалуйста или предложите другой вариант. Добавлено через 39 минут Помогите!!:wall: Добавлено через 11 минут вот такая вот игра #include <iostream.h> #include <conio.h> #include <windows.h> подробнее

Показать сообщение отдельно
cap
1 / 1 / 0
Регистрация: 27.10.2010
Сообщений: 51
23.01.2011, 23:51     Метод Рунге-Кутта 4-го порядка
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 <stdio.h>
#include <iostream>
#include <math.h>
using namespace std;
 
double f(double x, double y)
{
      return(2*y/x+x*x*pow(2.7,x));               
}
 
int main(void){
   int n,i;
   double x,y,h,end,k1,k2,k3,k4;
   
   x=-1;
   
   y=-2.33
   
   end=2
   
   h=0.3;   
 
   n=(end-x)/h;
   
   for(i=1;i<=n;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;
   }
 
}

почему когда я ставлю h=0,3 мы идем по x перескакивая через 0...значит деления на ноль не происходит но результат все равно неверен когда ставим интервал от 1 до 2 тоесть ноль вобще не затрагиваем все нормально

вот какие должны быть результаты при x=2 y=18.36



аналитическое решение для проверки Y(x)=x*x(pow(2.7,x)-2.7)
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 14:58. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru