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

UptoLike

259
p
13
t( ) = 3.06333747514295616 10
-9
, p
14
t( ) = -7.32039341683141883 10
-10
,
p
15
t( ) = -1.71023640935911516 10
-10
]
> add(rhs(Lpi[5][i]), i=2..n+m+2);
0.8421340936
> Lpi[6]:=Runge_Kutt(10.0);
Lpi
6
:= [t = 10.0, p
0
t( ) = 0.207549936565027016, p
1
t( ) = 0.254015629133477305,
p
2
t( ) = 0.154494241991991299, p
3
t( ) = 0.0616222310055043744, p
4
t( ) = 0.0174479708411122426,
p
5
t( ) = 0.00303304143426933381, p
6
t( ) = -0.000527245370352833116,
p
7
t( ) = 0.0000916526933278326690, p
8
t( ) = -0.0000159324270092257176,
p
9
t( ) = 0.00000276953284138237898, p
10
t( ) = -4.81491273482991492 10
-7
,
p
11
t( ) = 8.36344618647448142 10
-8
, p
12
t( ) = -1.46277343899484288 10
-8
,
p
13
t( ) = 2.41762204544735632 10
-9
, p
14
t( ) = -5.97759666906672786 10
-10
,
p
15
t( ) = -1.47557322166583962 10
-10
]
> add(rhs(Lpi[6][i]), i=2..n+m+2);
0.6977138845
> Lpi[7]:=Runge_Kutt(200.0);
Lpi
7
:= [t = 200.0, p
0
t( ) = 0.000163069438140189501, p
1
t( ) = 0.000199584909997749294,
p
2
t( ) = 0.000121364781452930250, p
3
t( ) = 0.0000484477166035522940,
p
4
t( ) = 0.0000136824135535196056, p
5
t( ) = 0.00000239166001333866370,
p
6
t( ) = -4.11414174644895684 10
-7
, p
7
t( ) = 7.29347065951768224 10
-8
,
p
8
t( ) = -1.22094080566235500 10
-8
, p
9
t( ) = 2.27310908562936344 10
-9
,
p
10
t( ) = -3.43346721494645653 10
-10
, p
11
t( ) = 7.48489358400078100 10
-11
,
p
12
t( ) = -6.74356127923232769 10
-12
, p
13
t( ) = 2.14095880501125717 10
-12
,
p
14
t( ) = 5.99255337388246958 10
-13
, p
15
t( ) = -5.09220877531093638 10
-13
]
> add(rhs(Lpi[7][i]), i=2..n+m+2);
0.0005481922310
> Lp1:=[seq(rhs(Lpi[t][2]), t=0..7)];
Lp1 := [1., 0.340947384189586922, 0.295555169740299606, 0.286213544236109174,
0.280490572394702653, 0.250510816421309290, 0.207549936565027016, 0.000163069438140189501]
> Lp2:=[seq(rhs(Lpi[t][3]), t=0..7)];
                                          -9                                       -10
         p13(t) = 3.06333747514295616 10 , p14(t) = -7.32039341683141883 10             ,

                                            -10
         p15(t) = -1.71023640935911516 10      ]


>   add(rhs(Lpi[5][i]), i=2..n+m+2);
                                                   0.8421340936

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

         p2(t) = 0.154494241991991299, p3(t) = 0.0616222310055043744, p4(t) = 0.0174479708411122426,

         p5(t) = 0.00303304143426933381, p6(t) = -0.000527245370352833116,

         p7(t) = 0.0000916526933278326690, p8(t) = -0.0000159324270092257176,

                                                                                     -7
         p9(t) = 0.00000276953284138237898, p10(t) = -4.81491273482991492 10 ,

                                          -8                                       -8
         p11(t) = 8.36344618647448142 10 , p12(t) = -1.46277343899484288 10 ,

                                          -9                                       -10
         p13(t) = 2.41762204544735632 10 , p14(t) = -5.97759666906672786 10             ,

                                            -10
         p15(t) = -1.47557322166583962 10      ]


>   add(rhs(Lpi[6][i]), i=2..n+m+2);
                                                   0.6977138845

>   Lpi[7]:=Runge_Kutt(200.0);
           Lpi7 := [t = 200.0, p0(t) = 0.000163069438140189501, p1(t) = 0.000199584909997749294,

                 p2(t) = 0.000121364781452930250, p3(t) = 0.0000484477166035522940,

                 p4(t) = 0.0000136824135535196056, p5(t) = 0.00000239166001333866370,

                                                     -7                                     -8
                 p6(t) = -4.11414174644895684 10 , p7(t) = 7.29347065951768224 10 ,

                                                     -8                                     -9
                 p8(t) = -1.22094080566235500 10 , p9(t) = 2.27310908562936344 10 ,

                                                      -10                                        -11
                 p10(t) = -3.43346721494645653 10          , p11(t) = 7.48489358400078100 10       ,

                                                      -12                                        -12
                 p12(t) = -6.74356127923232769 10          , p13(t) = 2.14095880501125717 10       ,

                                                     -13                                         -13
                 p14(t) = 5.99255337388246958 10          , p15(t) = -5.09220877531093638 10       ]


>   add(rhs(Lpi[7][i]), i=2..n+m+2);
                                                  0.0005481922310

>   Lp1:=[seq(rhs(Lpi[t][2]), t=0..7)];
Lp1 := [1., 0.340947384189586922, 0.295555169740299606, 0.286213544236109174,
       0.280490572394702653, 0.250510816421309290, 0.207549936565027016, 0.000163069438140189501]

>   Lp2:=[seq(rhs(Lpi[t][3]), t=0..7)];
                                                           259