Метод Рунге-Кутта и многое другое

@MrPhelko · Регистрация: 23.12.2012

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

У меня есть задание:

Комментарий модератора

5.18 Запрещено размещать задания в виде картинок и других файлов с их текстом.

надо проверить правильно ли я сделал ее?
Также нужна помощь с графиками, которые написанны в задании

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
51
52
53
54
55
56
57
58
59
#include <vcl.h>
#include <iostream.h>
#include <conio.h>
#include <math.h>
#pragma hdrstop
 
#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//---------------------------------------------------------------------------
 double f(double u, double a, double s, double x, double w, double d)
{
  return a*cos(s)+2*s*u+w*w*x;
}
 double f1(double u, double v)
{
  return u;
}
void __fastcall TForm1::Button1Click(TObject *Sender)
{
double a = Edit1->Text.ToDouble();
double u = Edit2->Text.ToDouble();
double s = Edit3->Text.ToDouble();
double x = Edit4->Text.ToDouble();
double w = Edit5->Text.ToDouble();
double d = Edit6->Text.ToDouble();
double v = Edit4->Text.ToDouble();
int N = Edit7->Text.ToDouble(), i;
double h = 0.1, k1, k2, k3, k4, m1, m2, m3, m4;
for(i = 0; i<N-1; i++)
   {
                   k1 = h*f(a, u, s, x, w, d);
                   k2 = h*f(u + h/2, x + (k1)/ 2, a, w, s, d);
                   k3 = h*f(u + h/2, x + (k2)/ 2, a, w, s, d);
                   k4 = h*f(u + h, x + k3, a, w, s, d);
                   x+= (k1 + 2*k2 + 2*k3 + k4)/6;
                   u+= h;
                   Series1->AddXY(u,x);
   }
      Edit8->Text = x;
for(i = 0; i<N-1; i++)
{
                   m1 = h*f1(u,v);
                   m2 = h*f1(u + h/2, v + (k1)/ 2);
                   m3 = h*f1(u + h/2, v + (k2)/ 2);
                   m4 = h*f1(u + h, v + k3);
                   v+= (k1 + 2*k2 + 2*k3 + k4)/6;
                   u+= h;
                   Series2->AddXY(u,v);
}
      Edit9->Text = v;
}

Так же я выложил сам исходник

@MrPhelko · 24.12.2012, 01:37 **[ТС]**

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
51
52
53
54
55
56
57
58
59
//---------------------------------------------------------------------------
 
#include <vcl.h>
#include <iostream.h>
#include <conio.h>
#include <math.h>
#pragma hdrstop
#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//---------------------------------------------------------------------------
double f(double u, double a, double s, double x, double w, double d)
{
  return a*cos(s)+2*s*u+w*w*x;
}
 double f1(double u, double v)
{
  return u;
}
 
void __fastcall TForm1::Button1Click(TObject *Sender)
{
double a = Edit1->Text.ToDouble();
double u = Edit2->Text.ToDouble();
double s = Edit3->Text.ToDouble();
double x = Edit4->Text.ToDouble();
double w = Edit5->Text.ToDouble();
double d = Edit6->Text.ToDouble();
double v = Edit4->Text.ToDouble();
int N = Edit7->Text.ToDouble(), i, E=5;
double k1, k2, k3, k4, m1, m2, m3, m4;
double h = 0.1;
for(i = 0; i<N-1; i++)
   {
                   k1 = h*f(a, u, s, x, w, d);
                   k2 = h*f(u + h/2, x + (k1)/ 2, a, w, s, d);
                   k3 = h*f(u + h/2, x + (k2)/ 2, a, w, s, d);
                   k4 = h*f(u + h, x + k3, a, w, s, d);
                   x+= (k1 + 2*k2 + 2*k3 + k4)/6;
                   u+= h;
                   m1 = h*f1(u,v);
                   m2 = h*f1(u + h/2, v + (k1)/ 2);
                   m3 = h*f1(u + h/2, v + (k2)/ 2);
                   m4 = h*f1(u + h, v + k3);
                   v+= (m1 + 2*m2 + 2*m3 + m4)/6;
                   u+= h;
                   double T = f(u,a,s,x,w,d);
                   Series1->AddXY(T,x+v);
   }
      Edit8->Text=printf("%1.2f \n", x);
      Edit9->Text = v;
}

Программу улучшил, теперь нужно как то построить гармоническое и затухающие колебание, помогите пожалуйста

Новые блоги и статьи Все статьи Все блоги /
Модель ЗдрввоСохранения 7: больше работников, больше ресурсов. anaschu 08.04.2026 работников и заданий может быть сколько угодно, но настроено всё так, что используется пока что только 20%	Дальние перспективы сервера - слоя сети с космологическим дизайном интефейса карты и логики. Hrethgir 07.04.2026 Дальнейшее ближайшее планирование вывело к размышлениям над дальними перспективами. И вот тут может быть даже будут нужны оценки специалистов, так как в дальних перспективах всё может очень сильно. . .	Горе от ума kumehtar 07.04.2026 Эта мне ментальная установка, что вот прямо сейчас, мол, мне для полного счастья не хватает (нужное вписать), и когда я этого достигну - тогда и полный кайф. Одна из самых сильных ловушек на пути. . . .	Использование значений реквизитов справочника в документе, с определенными условиями и правами Maks 07.04.2026 1. Контроль срока действия договора Алгоритм из решения ниже реализован на примере нетипового документа "ЗаявкаНаРаботу", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если. . .
Доступность команды формы по условию Maks 07.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: сделать доступной кнопку (команда формы "ЗавершитьСписание") при. . .	Уведомление о неверно выбранном значении справочника Maks 06.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "НарядПутевка", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если в документе выбран неверный склад. . .	Установка Qt Creator для C и C++: ставим среду, CMake и MinGW без фреймворка Qt 8Observer8 05.04.2026 Среду разработки Qt Creator можно установить без фреймворка Qt. Есть отдельный репозиторий для этой среды: https:/ / github. com/ qt-creator/ qt-creator, где можно скачать установщик, на вкладке Releases:. . .	AkelPad-скрипты, структуры, и немного лирики.. testuser2 05.04.2026 Такая программа, как AkelPad существует уже давно, и также давно существуют скрипты под нее. Тем не менее, прога живет, периодически что-то не спеша дополняется, улучшается. Что меня в первую очередь. . .