Форум программистов, компьютерный форум, киберфорум
Наши страницы
Pascal ABC
Войти
Регистрация
Восстановить пароль
 
manbo
1 / 1 / 0
Регистрация: 12.04.2012
Сообщений: 9
#1

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

30.05.2012, 23:00. Просмотров 959. Ответов 0
Метки нет (Все метки)

Здравствуйте, мне кажется у меня ошибка в коде, потому что значения не сходятся с методом Рунге-Кутты четвертого порядка.
Pascal
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
type
  vec = array [1..10] of real;
 
//---- описание переменных ----
var
  hV, hVt, eps, x, x9, h: real;
  i, j, n: integer;
  y: vec;
 
 
//---- значения функции ----
procedure rp(x: real; var y, f: vec);
begin
  f[1] := y[2];
  f[2] := (1 - sqr(y[1])) * y[2] - y[1];
end;
 
 
//---- метод Рунге - Кутты 3-го порядка ----
procedure rk3(n: integer; x, h: real; var y: vec);
var
  i, j: integer;
  h1, h2, q,k1,k2,k3: real;
  y0, y1, f: vec;
begin
  for i:=1 to n do 
    begin
      rp(x,y,f);
      k1:=h*f[i];
      
      for j:=1 to n do
          y0[i]:=y[i]+k1/3;
      rp(x+h/3,y0,f);
      k2:=h*f[i];
      
      for j:=1 to n do
          y1[i]:=y[i]+k2*2/3;
      rp(x+h*2/3,y1,f);
      k3:=h*f[i];
      
      y[i]:=y[i]+1/4*(k1+3*k3);
    end;
end;
 
//---- основная программа ----
begin
  //x - nachalo promezutka
  //x9 - konec promezutka
  //h - shag prosmotra
  //y[1],y[2] - nachalnie znacheniya funkcii
  //hV - shag vivoda znachenii
  //hVt - nachalnoe znachenie hV
  
  x:=0;
  x9:=10;
  h:=0.01;
  y[1]:=1;
  y[2]:=2;
  
  hV := 0.5;
  hVt := hV;
  eps := 0.00001;
  while (x < x9) = (h > 0) do
  begin
    x := x + h;
    rk3(2,x,h,y);
    if abs(x - hV) < eps then writeln('x= ', x:3:3, ' y[1]= ', y[1]:3:3, ' y[2]= ', y[2]:3:3);
    if x > hV then hV := hV + hVt;
  end;
end.

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

0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
30.05.2012, 23:00
Ответы с готовыми решениями:

Заданы матрицы А и В третьего порядка
Заданы матрицы А и В третьего порядка. Переменной S присвоить -1, если ...

Диффур второго порядка - Метод Рунге Кутты (уравнение 2го порядка)
Ребята, помогите разобраться. Нужно решить дифференциальное уравнение второго...

Заданы матрицы А и В третьего порядка
Заданы матрицы А и В третьего порядка. Переменной S присвоить -1, если...

Метод Рунге-Кутта 4 порядка точности
Уважаемые Форумчане помогите решить задачу Коши для обыкновенного...

Решение уравнений методом Рунге-Кутта 4 порядка
напишите програму для решения уравнений методом Рунги-кута 4 порядка))

0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
30.05.2012, 23:00

Решение диф. уравнения 1-го порядка методом Рунге-Кутта
Здраствуйте, извините за обращение, задание поди пустяковое, но и с ним мне...

Метод Рунге-Кутта-Мерсона для дифференциальных уравнений 2 порядка
Добрый день. Помогите пожалуйста понять в чем ошибка. У меня дано такое...

Решение дифференциального уравнения методом Рунге-Кутта четвертого порядка
Решить дифференциальное уравнение y'-(y/x) = x*sin(x) X принадлежит ...


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

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

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