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

@городецкий ли · Регистрация: 25.01.2016

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

Помогите найти ошибку в методе рунге-Кутта 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
36
37
38
39
40
41
System::System(double m, const Vector3D& g, double C, double f, double p, double S, double t, const Vector3D& w, const Vector3D& v, const Vector3D& r)
    : m(m), g(g), C(C), f(f), p(p), S(S), t(t), w(w), v(v), r(r)
{
}
Vector3D System::dr() const
{
    return v;
}
 
Vector3D System::dv() const
{
    return -1*g - 0.5*C*p*S*v.modul()*v/m+ 2*v*w;
}
 
 
void System::step(double dt)
{
    Vector3D l1, l2, l3, l4, k1, k2, k3, k4;
    System sys2 = *this;
    k1 = sys2.dr()*dt;
    l1 = sys2.dv()*dt;
    sys2.r = this->r + k1/2.;
    sys2.v = this->v + l1/2.;
 
    k2 = sys2.dr()*dt;
    l2 = sys2.dv()*dt;
    sys2.r = this->r + k2/2.;
    sys2.v = this->v + l2/2.;
 
    k3 = sys2.dr()*dt;
    l3 = sys2.dv()*dt;
    sys2.r = this->r + k3;
    sys2.v = this->v + l3;
 
    k4 = sys2.dr()*dt;
    l4 = sys2.dv()*dt;
 
    r = r+(1./6.)*(l1+2*l2+2*l3+l4);
    v = v+(1./6.)*(k1+2*k2+2*k3+k4);
    t = t+dt;
}

@volvo · 06.06.2017, 16:36

городецкий ли, это что, действительно проект C++/CLI, или это обычный С++?

@городецкий ли · 06.06.2017, 16:51 **[ТС]**

Это часть проекта в Qt

Новые блоги и статьи Все статьи Все блоги /
Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек SDL3 и Box2D из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия SDL 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual. . .	SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .	Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .
Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .	И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»

@volvo Супер-модератор 33395 / 21505 / 8236 Регистрация: 22.10.2011 Сообщений: 36,902 Записей в блоге: 12
	06.06.2017, 16:36
	городецкий ли, это что, действительно проект C++/CLI, или это обычный С++? 0

@городецкий ли 0 / 0 / 0 Регистрация: 25.01.2016 Сообщений: 18
	06.06.2017, 16:51 [ТС]
	Это часть проекта в Qt 0