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

@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;
}

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

Новые блоги и статьи Все статьи Все блоги /
Как украсить новогоднюю елку с Q# и Qiskit EggHead 24.06.2025 Что может быть необычнее, чем применить законы квантовой механики для украшения новогодней елки? Пока другие развешивают обычные гирлянды, я решил объединить свою страсть к квантовым вычислениям с. . .	Системы нулевого доверия на C# UnmanagedCoder 24.06.2025 Традиционная архитектура безопасности работает по принципу средневекового замка: создаём высокие стены вокруг корпоративной сети, укрепляем ворота межсетевыми экранами и системами обнаружения. . .	Снова не мой путь. Циклическое среднее, я обеими руками за проверку условия, в ракурсе данной задачи - циклическое среднее в топку. Hrethgir 24.06.2025 Привет. Такой вопрос - нужно выводить среднее математическое между двумя направлениями, интервал значений которых может лежать в диапазоне одного оборота по кругу. Проблема заключается в том, что. . .	Деплой Flask приложения py-thonny 23.06.2025 За годы работы с Flask я натыкался на одни и те же грабли достаточно часто, чтобы наконец научится их обходить. И сегодня хочу поделится опытом, который сбережет вам немало нервных клеток. Начнем с. . .	WebAssembly и контейнеры в .NET Aspire для оркестрации распределенных архитектур ArchitectMsa 23.06.2025 Я наблюдаю, как WebAssembly (или просто WASM) постепенно выходит за рамки своего первоначального предназначения — исполнения кода на стороне браузера. Теперь эта технология проникает в серверную. . .
Непрерывная интеграция для пакета Python Mr. Docker 22.06.2025 Было 4 часа утра пятницы, когда я выпустил новую версию нашей внутренней библиотеки для обработки данных. Релиз 0. 5. 2 содержал небольшой фикс для обработки дат в ISO формате, что может пойти не так?. . .	Продвинутый ETL на C# из OLTP БД в хранилище stackOverflow 22.06.2025 Работая в сфере корпоративной аналитики, я постоянно сталкиваюсь с одним и тем же - нужны чистые, структурированные и, главное, свежие данные. Без них современные аналитические системы, машинное. . .	Мастер-класс по микросервисам на Node.js Reangularity 21.06.2025 Node. js стал одной из самых популярных платформ для микросервисной архитектуры не случайно. Его неблокирующая однопоточная модель и событийно-ориентированный подход делают его идеальным для. . .	Управление Arduino из WPF приложения Wired 21.06.2025 Зачем вообще связывать Arduino с WPF-приложением? Казалось бы, у Arduino есть собственная среда разработки, своя экосистема, свои способы управления. Однако при создании серьезных проектов. . .	Звёздная пыль kumehtar 20.06.2025 Я просто это себе представляю: как создавался этот мир. Как энергия слипалась в маленькие частички. Как они собирались в первые звёзды, как во вселенной впервые появился Свет. Как эти звёзды. . .

Опции темы