Теория массового обслуживания. Сивохин А.В - 184 стр.

UptoLike

ВУЗ: 

ПГУ | Пенза

Составители: 

Рубрика: 

184
Runge_Kutt := proc x_rkf45() ... end proc;
> Lpi[0]:=Runge_Kutt(0.0);
Lpi
0
:= t = 0., p0 t( ) = 1., p1 t( ) = 0., p2 t( ) = 0., p3 t( ) = 0., p4 t( ) = 0.[]
> Lpi[1]:=Runge_Kutt(0.5);
Lpi
1
:= [t = 0.5, p0 t( ) = 0.604093554513178476, p1 t( ) = 0.304510299575287080,
p2 t( ) = 0.0767805796419036014, p3 t( ) = 0.0129438608656616324, p4 t( ) = 0.00167170540396938390]
> Lpi[2]:=Runge_Kutt(1.0);
Lpi
2
:= [t = 1.0, p0 t( ) = 0.555658125930388880, p1 t( ) = 0.326663585088933228,
p2 t( ) = 0.0960480760949048301, p3 t( ) = 0.0188454668683919364, p4 t( ) = 0.00278474601738116884]
> Lpi[3]:=Runge_Kutt(1.5);
Lpi
3
:= [t = 1.5, p0 t( ) = 0.548066674560359844, p1 t( ) = 0.329794482872266870,
p2 t( ) = 0.0992320253266893438, p3 t( ) = 0.0199089714127403186, p4 t( ) = 0.00299784582794319172]
> Lpi[4]:=Runge_Kutt(2.0);
Lpi
4
:= [t = 2.0, p0 t( ) = 0.546834235542790648, p1 t( ) = 0.330295142835436706,
p2 t( ) = 0.0997524809033172588, p3 t( ) = 0.0200847947882790631, p4 t( ) = 0.00303334593017586972]
> Lpi[5]:=Runge_Kutt(5.0);
Lpi
5
:= [t = 5.0, p0 t( ) = 0.546594054352584346, p1 t( ) = 0.330392499854962551,
p2 t( ) = 0.0998540230091448656, p3 t( ) = 0.0201191280174715312, p4 t( ) = 0.00304029476583665042]
> Lpi[6]:=Runge_Kutt(10.0);
Lpi
6
:= [t = 10.0, p0 t( ) = 0.546594052466686798, p1 t( ) = 0.330392527294737781,
p2 t( ) = 0.0998539639307437638, p3 t( ) = 0.0201191807088391406, p4 t( ) = 0.00304027559899271117]
> Lpi[7]:=Runge_Kutt(20.0);
Lpi
7
:= [t = 20.0, p0 t( ) = 0.546594055638341070, p1 t( ) = 0.330392511169330827,
p2 t( ) = 0.0998539959020817298, p3 t( ) = 0.0201191510889340460, p4 t( ) = 0.00304028620131242475]
> Lp1:=[seq(rhs(Lpi[t][2]), t=0..7)];
Lp1 := [1., 0.604093554513178476, 0.555658125930388880, 0.548066674560359844,
0.546834235542790648, 0.546594054352584346, 0.546594052466686798, 0.546594055638341070]
> Lp2:=[seq(rhs(Lpi[t][3]), t=0..7)];
Lp2 := [0., 0.304510299575287080, 0.326663585088933228, 0.329794482872266870,
0.330295142835436706, 0.330392499854962551, 0.330392527294737781, 0.330392511169330827]
> Lp3:=[seq(rhs(Lpi[t][4]), t=0..7)];
Lp3 := [0., 0.0767805796419036014, 0.0960480760949048301, 0.0992320253266893438,
0.0997524809033172588, 0.0998540230091448656, 0.0998539639307437638, 0.0998539959020817298
                                    Runge_Kutt := proc(x_rkf45) ... end proc;

>    Lpi[0]:=Runge_Kutt(0.0);
                       Lpi0 := [t = 0., p0(t) = 1., p1(t) = 0., p2(t) = 0., p3(t) = 0., p4(t) = 0.]


>    Lpi[1]:=Runge_Kutt(0.5);
    Lpi1 := [t = 0.5, p0(t) = 0.604093554513178476, p1(t) = 0.304510299575287080,

         p2(t) = 0.0767805796419036014, p3(t) = 0.0129438608656616324, p4(t) = 0.00167170540396938390]

>    Lpi[2]:=Runge_Kutt(1.0);
    Lpi2 := [t = 1.0, p0(t) = 0.555658125930388880, p1(t) = 0.326663585088933228,

         p2(t) = 0.0960480760949048301, p3(t) = 0.0188454668683919364, p4(t) = 0.00278474601738116884]

>    Lpi[3]:=Runge_Kutt(1.5);
    Lpi3 := [t = 1.5, p0(t) = 0.548066674560359844, p1(t) = 0.329794482872266870,

         p2(t) = 0.0992320253266893438, p3(t) = 0.0199089714127403186, p4(t) = 0.00299784582794319172]

>    Lpi[4]:=Runge_Kutt(2.0);
    Lpi4 := [t = 2.0, p0(t) = 0.546834235542790648, p1(t) = 0.330295142835436706,

         p2(t) = 0.0997524809033172588, p3(t) = 0.0200847947882790631, p4(t) = 0.00303334593017586972]

>    Lpi[5]:=Runge_Kutt(5.0);
    Lpi5 := [t = 5.0, p0(t) = 0.546594054352584346, p1(t) = 0.330392499854962551,

         p2(t) = 0.0998540230091448656, p3(t) = 0.0201191280174715312, p4(t) = 0.00304029476583665042]

>    Lpi[6]:=Runge_Kutt(10.0);
    Lpi6 := [t = 10.0, p0(t) = 0.546594052466686798, p1(t) = 0.330392527294737781,

         p2(t) = 0.0998539639307437638, p3(t) = 0.0201191807088391406, p4(t) = 0.00304027559899271117]

>    Lpi[7]:=Runge_Kutt(20.0);
    Lpi7 := [t = 20.0, p0(t) = 0.546594055638341070, p1(t) = 0.330392511169330827,

         p2(t) = 0.0998539959020817298, p3(t) = 0.0201191510889340460, p4(t) = 0.00304028620131242475]

>    Lp1:=[seq(rhs(Lpi[t][2]), t=0..7)];
    Lp1 := [1., 0.604093554513178476, 0.555658125930388880, 0.548066674560359844,
          0.546834235542790648, 0.546594054352584346, 0.546594052466686798, 0.546594055638341070]

>    Lp2:=[seq(rhs(Lpi[t][3]), t=0..7)];
    Lp2 := [0., 0.304510299575287080, 0.326663585088933228, 0.329794482872266870,
          0.330295142835436706, 0.330392499854962551, 0.330392527294737781, 0.330392511169330827]

>    Lp3:=[seq(rhs(Lpi[t][4]), t=0..7)];
Lp3 := [0., 0.0767805796419036014, 0.0960480760949048301, 0.0992320253266893438,
        0.0997524809033172588, 0.0998540230091448656, 0.0998539639307437638, 0.0998539959020817298

                                                          184

Страницы