Технический форум - Метод Рунге-Кутта 2-го порядка, Паскаль

Технический форум (http://www.tehnari.ru/)

- Помощь студентам (http://www.tehnari.ru/f41/)

- - Метод Рунге-Кутта 2-го порядка, Паскаль (http://www.tehnari.ru/f41/t96234/)

Xstyle74

09.05.2014 05:18

Метод Рунге-Кутта 2-го порядка, Паскаль

Вложений: 1

Помогите с заданием, плиииззз, вообще ниче понять не могу,
Задание
Ток в электрической цепи описывается дифференциальным уравнением (см. стр. 24), которое необходимо решить, найдя y(0.1), y(0.2) ... y(1) с помощью метода Рунге-Кутта второго порядка с коррекцией по средней производной (в дифференциальном уравнении k = 6 По найденным значениям y определить количество тепла, выделяющееся на единичном сопротивлении за единицу времени, по формуле с помощью метода трапеций.
Уравнение со стр.24

Vladimir_S

09.05.2014 08:28

Цитата:

Сообщение от Xstyle74 (Сообщение 1030872)

Помогите с заданием, плиииззз, вообще ниче понять не могу

Пожалуйста. У нас разбиралась абсолютно аналогичная задача, только функция другая и метод интегрирования надо поправить: там - Симпсон, у Вас - трапеции. http://www.tehnari.ru/f43/t95614/
От поста #5 и дальше.

Xstyle74

09.05.2014 21:29

Если бы я разбирался еще...если не сложно то помогите плиз,ну паскаль вообще не дается мне....

Vladimir_S

09.05.2014 22:36

Цитата:

Сообщение от Xstyle74 (Сообщение 1030997)

Если бы я разбирался еще...если не сложно то помогите плиз,ну паскаль вообще не дается мне....

Очень это грустно. На самом деле переделка той программы заняла полторы-две минуты. Ладно, вот вам для тупого бездумного скатывания и сдавания.

Код:

Const

 h=0.1;

 N=10;



Var

 I:Array[0..N] of Real;

 j:integer;

 t,Z,Sum:Real;



Function F(x,y:real):real;

begin

 F:=(6.0-y*y)*Cos(x)+y*2;

end;



Begin

 I[0]:=0.3;

 t:=0;

 writeln('t= ',t:3:1,'   I= ',I[0]:6:3);

 for j:=1 to N do

  begin

   t:=t+h;

   Z:=I[j-1]+h*F(t-h,I[j-1]);

   I[j]:=I[j-1]+h*(F(t-h,I[j-1])+F(t,Z))/2;

   writeln('t= ',t:3:1,'   I= ',I[j]:6:3);

  end;

 writeln;

 Sum:=(Sqr(I[0])+Sqr(I[N]))/2*h;

 for j:=1 to N-1 do

  Sum:=Sum+h*Sqr(I[j]);

 Writeln('Heat: ', Sum:0:5);

 Readln

End.

Xstyle74

10.05.2014 00:13

Прибольшущее тебе спасибо, да пусть в твоем доме интернет будет без трабл, и потерь пакетов, Говарю тебе как работник одной из компании провайдера интернета

Vladimir_S

10.05.2014 08:28

Цитата:

Сообщение от Xstyle74 (Сообщение 1031032)

Спасибо. И да, поправил одну неточность в части интегрирования: сложение квадратов первой и последней точек (было сложение самих токов). Приношу извинения - поторопился.

Часовой пояс GMT +4, время: 22:19.