Форум программистов, компьютерный форум, киберфорум
Delphi для начинающих
Войти
Регистрация
Восстановить пароль
 
xybr
1

Решение диф. уравнений 2-го порядка методом эйлера

17.12.2013, 20:26. Просмотров 766. Ответов 0
Метки нет (Все метки)

помогите ввести 2 параметра: омега и гамма в алгоритме эйлера в функцию ye[i+1]
условие:
ye[i+1]:=ye[i]+h-2*ye[i]*гамма - (омега^2)*xe[i]
не могу сообразить что мне прописать чтобы они выводились в программе также как параметры a и h



Delphi
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
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, AxCtrls, OleCtrls, VCF1, Math, TeEngine, Series,
  ExtCtrls, TeeProcs, Chart, XPMan, OleCtnrs;
 
type
  TForm1 = class(TForm)
    Button1: TButton;
    F1Book1: TF1Book;
    Edit1: TEdit;
    Edit3: TEdit;
    Button2: TButton;
    Chart1: TChart;
    Series1: TLineSeries;
    Series2: TLineSeries;
    XPManifest1: TXPManifest;
    Label2: TLabel;
    Label4: TLabel;
    Label6: TLabel;
    procedure Button1Click(Sender: TObject);
    procedure Button2Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
var
  Form1: TForm1;
 
implementation
 
{$R *.dfm}
 
function f1(x:real):real;
begin
f1:=-2*sqr(Cos(x))+Cos(x);      //первое уравнение
end;
 
function f2(x:real;y:real):real;
begin                              //второе уравнение
f2:=Sin(2*x)-y*(Tan(x));
end;
 
procedure TForm1.Button1Click(Sender: TObject);
type mass=array [1..20] of real;
var a,h:real;
xe,ye:mass;         //объявление переменных
i:integer;
begin
 
 
 
begin
a:=strtofloat(Edit1.Text);  //ввод начального значения и шага
h:=strtofloat(Edit3.Text);
 
xe[1]:=a;
ye[1]:=f1(xe[1]);
 
for i:=1 to 11 do
begin
ye[i+1]:=ye[i]+h-2*ye[i]); //алгоритм Эйлера для уравнения 1-ого порядка
xe[i+1]:=xe[i]+h*f1;
f1book1.NumberRC[i,1]:=ye[i]; //вывод в F1book
end;
 
 end;
 
end;
 
procedure TForm1.Button2Click(Sender: TObject);
type mass=array [1..20] of real;
var a,h:real;
x:mass;
 xe,ye:mass;
   i:integer;
begin
a:=strtofloat(Edit1.Text);
h:=strtofloat(Edit3.Text);
x[1]:=a;
for i:=1 to 11 do
begin
Chart1.Series[0].AddXY(x[i],f1(x[i])) ; //рисование графика
x[i+1]:=x[i]+h;
end;
 
xe[1]:=a;
ye[1]:=f1(xe[1]);
 
for i:=1 to 11 do
begin
ye[i+1]:=ye[i]+h*f2(xe[i],ye[i]);
xe[i+1]:=xe[i]+h;
 
Chart1.Series[1].AddXY(xe[i],ye[i]) ;  //рисование графика
end;
end;
 
end.
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
17.12.2013, 20:26
Ответы с готовыми решениями:

Решение дифференциальных уравнений четвертого порядка методом Эйлера и методом Рунге-Кутта
Форумчане прошу помочь решить уравнение f :=y+y'*x+y''+y'''*x методами Рунге-Кутта и Эйлера. Если...

Решение дифференциальных уравнений 1 порядка методом Эйлера
Вот я написал программу, которая решает дифференциальные уравнения 1 порядка методом Эйлера (в...

Решение систем дифференциальных уравнений N-го порядка методом Эйлера
Добрый день. Помогите решить дифференциальное уравнение первого порядка методом Эйлера {x-2x+y=0...

Решение системы дифференциальных уравнений методом Эйлера и Рунге-Кутты 4 порядка
Здравствуйте уважаемые форумчани. Есть следующая задача: решить методом Эйлера и Рунге-Кутты 4...

0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
17.12.2013, 20:26

Заказываю контрольные, курсовые, дипломные и любые другие студенческие работы здесь.

Решение задачи Коши для системы 2 диф. уравнений методом Рунге-Кутты 4 порядка
Здравствуйте, помогите с написанием это программой, вообще не знаю как писать.

Система диф.уравнений методом Эйлера
Всем привет,столкнулся с проблемой решения диф уравнений програмным блоком методом Эйлера.Был бы...

Решение диф. уравнения 2 порядка методом Рунге-Кутты 4 порядка точности
u' = v, f == v', g == u' f := 2*u + v - 1 - 2*x g := v Subscript = 0, Subscript = 3, Subscript =...

Решение диф уравнений первого порядка.
нужна помощь при решение диф уравнений

Решение системы из 4 неленейных диф. уравнений второго порядка
Добрый. Проблема заключается в том, что необходимо написать код на C++ или Фортран для решения...

Решение дифференциального уравнения методом Эйлера и методом Рунге-кутта 4 порядка
Помогите пожалуйста решить уравнение y''-4y'+5y=2x2ex , методом Эйлера и методом Рунге-кутта 4...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
1
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.