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

@Ирина0101 · Регистрация: 21.11.2013

Author24 — интернет-сервис помощи студентам

Товарищи-знатоки помогите пожалуйста в решении следующей задачи:

Применяя метод Рунге-Кутта 4го порядка точности, найти решение задачи в следующих точках 0.1,0.3,0.4. Нарисовать фазовый портрет на плоскости. Вычисление провести с точностью 0.01, t<=0.5.

$https://www.cyberforum.ru/cgi-bin/latex.cgi?\frac{dx1}{dt}=-x1*x2$
$https://www.cyberforum.ru/cgi-bin/latex.cgi?\frac{dx2}{dt}={x1}^{2}+{x2}^{2}$
$https://www.cyberforum.ru/cgi-bin/latex.cgi?x1(0)=-2$
$https://www.cyberforum.ru/cgi-bin/latex.cgi?x2(0)=0$

Мой код не доработан,помогите исправить:

C++

#include<iostream>
#include<cmath>
#include<cstdlib>
using namespace std;
 
double f1 (double t,double y, double x)
{
       double a=(-1)*x*y;
       return a;          
}
double f2 (double t,double y, double x)
{
       double a=x*x+y*y;
       return a;          
}
 
int main()
{ 
    double x,y,x0=-2,y0=0,h=0.1,n=10,t=0.5;
    double k0,k1,k2,s0,s1,s2,X,Y,k3,s3;
    double x_n=-2,y_n=0;
    double eps=0.01;
    do
    {
        k0=h*f1(t,y0,x0);
        s0=h*f2(t,y0,x0);
        k1=h*f1(t+h/2,y0+k0/2,x0+s0/2);
        s1=h*f2(t+h/2,y0+k0/2,x0+s0/2);
        k2=h*f1(t+h/2,y0+k0/2,x0+s0/2);
        s2=h*f2(t+h/2,y0+k0/2,x0+s0/2);
        k3=h*f1(t+h,y0+k0,x0+s0);
        s3=h*f2(t+h,y0+k0,x0+s0);
        y=y0+(k0+2*k1+2*k2+k3)/6;
        x=x0+(s0+2*s1+2*s2+s3)/6;
        y0=y;
        x0=x;
        cout<<"t="<<t<<" "<<"y="<<y<<" "<<"x="<<x<<endl;
        t=t+h;
    }
    while(t<=n);
    {
        X=x;
        x0=x_n;
        Y=y;
        y0=y_n;
        h=h/2;
        cout<<X<<endl;
        t=0;
        do
        {
            k0=h*f1(t,y0,x0);
                s0=h*f2(t,y0,x0);
                k1=h*f1(t+h/2,y0+k0/2,x0+s0/2);
                s1=h*f2(t+h/2,y0+k0/2,x0+s0/2);
                k2=h*f1(t+h/2,y0+k0/2,x0+s0/2);
                s2=h*f2(t+h/2,y0+k0/2,x0+s0/2);
                k3=h*f1(t+h,y0+k0,x0+s0);
                s3=h*f2(t+h,y0+k0,x0+s0);
                y=y0+(k0+2*k1+2*k2+k3)/6;
                x=x0+(s0+2*s1+2*s2+s3)/6;
            y0=y;
            x0=x;
            cout<<"t="<<t<<" "<<"y="<<y<<" "<<"x="<<x<<endl;
            t=t+h;
        }
        while(t<=n);
    }
    while((fabs(x-X)<eps)||(fabs(y-Y)<eps))
    cout<<"Y="<<y<<" "<<"X="<<x<<endl;
    system("pause");
    return 0;
}

У меня происходит зацикливание,ребята помогите)

Вернуться к обсуждению:
Метод Рунге-Кутта 4 порядка C++

Другие темы раздела
C++ Задача про парники https://www.cyberforum.ru/ cpp-beginners/ thread1029045.html Здравствуйте. Задача по программированию, но, как мне кажется, у нее есть некий математический алгоритм решения. Задано поле n x m квадратных ячеек, в каждой из которых могут находиться посадки огурцов. Необходимо построить парники, закрывающие огурцы. Парники могут быть только прямоугольной формы, только со сторонами, параллельными сторонам поля. Стоимость строительства одного парника...	C++ Реализовать динамический массив структур внутри динамической структуры Добрый вечер! Я никак не могу реализовать динамический массив структур внутри динамической структуры. Может кто нибудь ссылочку скинет на материал(простое гугление не помогло). Заранее спасибо
C++ Даны три целых числа. Определить, для которого из них сумма цифр больше Буду благодарен) https://www.cyberforum.ru/ cpp-beginners/ thread1029030.html	C++ переделать часть программы в функцию https://www.cyberforum.ru/ cpp-beginners/ thread1029025.html Вот сама программа: #include <conio.h> #include <iostream.h> #include <locale.h> #include <stdio.h> #include <string.h> int main() {int i, j, n, b,d, flag, c; char isk; char dan;
C++ В массиве из 10 целых чисел подсчитать сумму элементов, стоящих левее максимального Помогите решить задачку! Препод ужасно объясняет, ничего не понимаю! В массиве из 10 целых чисел подсчитать сумму элементов, стоящих левее максимального.	C++ Нужно Нарисовать Треугольник на C++ https://www.cyberforum.ru/ cpp-beginners/ thread1029003.html Дали задание по C++. Нарисовать треугольник.Дайте материалы(ссылки,книги,tutorial) с не сложными кодами. Заранее Спасибо p.s. Code::Blokcs на Windows 8 x64
C++ Хеширование Добрый день! В общем, есть задание, нужно какое-то n-е кол-во ФИО, причем фиксированное, т.е нужно в самом коде эти данные и написать. Далее применить хеширование и определить кол-во синонимных пар и кол-во разреженности (т.е как Я понял, сколько в хеш таблице пустых строк из-за этих синонимных пар). У меня сразу же появился вопрос, как лучше написать эти данные (ФИО), для того, чтобы удобнее... https://www.cyberforum.ru/ cpp-beginners/ thread1028996.html	найти номер строки и номер столбца, в которых находится наименьший элемент C++ Ребят, очень надеюсь на вашу помощь! Массив размерностью MxN. Необходимо найти номер строки и номер столбца, в которых находится наименьший элемент. Во всех вариантах необходимо написать три функции, которые будут вызываться из функции main(). Первая функция - получает размерность массива, создает динамический массив и возвращает указатель на начало созданного массива. Вторая функция –...
C++ Почему в конце символьной строки появляются разные символы, и что с ними делать? https://www.cyberforum.ru/ cpp-beginners/ thread1028971.html В строке, содержащей последовательность слов, найти конец предложения, обозначенный символом "точка". В следующем слове первую строчную букву заменить на прописную. #include<conio.h> #include<stdio.h> #include<string.h> int main() {char s; int Len; puts ("Input string:"); gets(s);	C++ Даны площади квадратов https://www.cyberforum.ru/ cpp-beginners/ thread1028965.html Даны площади квадратов. Найти длину диагонали самого большого из них. Заранее спасибо.
C++ В файле задан текст. Определить частоту повторяемости входящих в текст цифр В файле задан текст. Определить частоту повторяемости входящих в текст цифр.	C++ Как убрать заливку в TChart? https://www.cyberforum.ru/ cpp-beginners/ thread1028951.html Нужно убрать заливку внутри графика, чтобы был только контур

	04.12.2013, 23:36