-
Notifications
You must be signed in to change notification settings - Fork 0
/
sir_model.nb
1479 lines (1469 loc) · 75.6 KB
/
sir_model.nb
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
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 11.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 77222, 1469]
NotebookOptionsPosition[ 76855, 1452]
NotebookOutlinePosition[ 77203, 1467]
CellTagsIndexPosition[ 77160, 1464]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"b", "=", "0.0005"}], ";",
RowBox[{"k", "=", "0.01"}], ";"}], "\n",
RowBox[{
RowBox[{"system", "=",
RowBox[{"{",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"s", "'"}], "[", "t", "]"}], "\[Equal]",
RowBox[{
RowBox[{"-", "b"}], " ",
RowBox[{"s", "[", "t", "]"}], " ",
RowBox[{"i", "[", "t", "]"}]}]}], ",",
RowBox[{
RowBox[{
RowBox[{"i", "'"}], "[", "t", "]"}], "\[Equal]",
RowBox[{
RowBox[{"b", " ",
RowBox[{"s", "[", "t", "]"}], " ",
RowBox[{"i", "[", "t", "]"}]}], "-",
RowBox[{"k", " ",
RowBox[{"i", "[", "t", "]"}]}]}]}], ",",
RowBox[{
RowBox[{
RowBox[{"r", "'"}], "[", "t", "]"}], "\[Equal]",
RowBox[{"k", " ",
RowBox[{"i", "[", "t", "]"}]}]}], ",",
RowBox[{
RowBox[{"s", "[", "0", "]"}], "\[Equal]", "1000"}], ",",
RowBox[{
RowBox[{"i", "[", "0", "]"}], "\[Equal]", "10"}], ",",
RowBox[{
RowBox[{"r", "[", "0", "]"}], "\[Equal]", "0"}]}], "}"}]}], ";"}], "\n",
RowBox[{"sol", "=",
RowBox[{"NDSolve", "[",
RowBox[{"system", ",",
RowBox[{"{",
RowBox[{"s", ",", "i", ",", "r"}], "}"}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "100"}], "}"}]}], "]"}]}], "\n",
RowBox[{"Plot", "[",
RowBox[{
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"s", "[", "t", "]"}], ",",
RowBox[{"i", "[", "t", "]"}], ",",
RowBox[{"r", "[", "t", "]"}]}], "}"}], "/.", "sol"}], "]"}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "100"}], "}"}]}], "]"}]}], "Input",
CellChangeTimes->{{3.6980865750296764`*^9, 3.698086583554807*^9}}],
Cell[BoxData[
RowBox[{"{",
RowBox[{"{",
RowBox[{
RowBox[{"s", "\[Rule]",
TagBox[
TemplateBox[{RowBox[{
StyleBox[
TagBox["InterpolatingFunction", "SummaryHead"],
"NonInterpretableSummary"],
StyleBox["[", "NonInterpretableSummary"],
DynamicModuleBox[{Typeset`open$$ = False},
PanelBox[
PaneSelectorBox[{False -> GridBox[{{
PaneBox[
ButtonBox[
DynamicBox[
FEPrivate`FrontEndResource[
"FEBitmaps", "SquarePlusIconMedium"]],
ButtonFunction :> (Typeset`open$$ = True), Appearance ->
None, Evaluator -> Automatic, Method -> "Preemptive"],
Alignment -> {Center, Center}, ImageSize ->
Dynamic[{
Automatic, 3.5 CurrentValue["FontCapHeight"]/
AbsoluteCurrentValue[Magnification]}]],
GraphicsBox[{{{{}, {},
TagBox[{
Directive[
Opacity[1.],
RGBColor[0.368417, 0.506779, 0.709798],
AbsoluteThickness[1]],
LineBox[CompressedData["
1:eJwBUQOu/CFib1JlAgAAADQAAAACAAAAwxW2NglTLkAIPrPOwd9LQBCpiU0A
Zy5ARY3GM30xS0BYwHE0uPYwQPK0mhYLeDhAVeKn9fSbMkAXM4rVYiAnQPGM
XlDoODRALjmry0tJFkAMZYQQ9vg1QFuGC+ExYARAVXMzhAibN0D4T5dWA9fz
Px6vUV01YDlAKo14uytp4j+Fc/DPGB07QMkmkHlHjNE/GW4Y9gC8PEDh71KI
YcjBPy2Wr4EDfj5AP2ZeGy05sT84+mdgBRFAQLSg5RkMtaE/GcCvMpb0QEBT
C6nxam6RP0nKt1EC1EFAKL8Ntk+YgT+Pb4TKcKRCQA6worEay3I/lasIdmyG
Q0Bqp9BQdkJjP7KCUXtqWURAz2Yb0g/aVD8enlrNQyhFQPAzN/urEkc/SlAb
UqoIRkAGV0TowYo4P42doDAT2kZAkGW/0vp/Kz+Pgd1BCb1HQMsI/ZXJqh0/
4anan9qbSEAii1vULlsQP0ltnFeua0lAcWeiDFTyAj9xxxVCD01KQIBpBDIb
K/U+sLxThnIfS0Aj/9247s7oPj72URex7UtAIZpQ2SGi3T6MxgfbfM1MQOL7
pgDVGNE+8DGC+EqeTUCdhCSRmKXEPhQ0tEimgE5AomZ9T2cXuD5P0aryA1RP
QBHPPEayba0+bdmwdJ4RUEAkpFJpYkqiPpIVaIkBglBAzVzuHA0Olj5CnwHL
5epQQGRNKnoXzIs+UvT2pZBcUUAC07KT3+GAPoprTCcpzFFAK4fPv0PTdD5N
MITVQjRSQE9WWaE3fGo+cMAXHSOlUkBhvNwV90pgPh+ejZGEDlNA0n/WY0xJ
VT71nWOs03VTQIn/116A3U4+K2mVYOnlU0D2YQYVRfdFPuyBqUGATlRAL/BJ
ELHoPT4NZhm83b9UQDY1xmLbnDE+Vmzp3CgvVUC0cv2nlPsmPirAmyr1llVA
6m/4iB4yIz5e36kRiAdWQJkgNlNCCys+HUyaJZxwVkBMq2/vTQs2PjyE5tJ2
4lZA9/3bNXLFPT6D3pImP1JXQDos8DmWmz4+VYYhp4i6V0Cr9sJy2mo5Pof5
C8GYK1hA4hLyj6dcLz5FutgHKpVYQEi8i5rbbRc+Phnk+P//WEBvVfy7623s
PbthcMc=
"]]},
Annotation[#,
"Charting`Private`Tag$2730#1"]& ]}}, {}, {}}, {
DisplayFunction -> Identity,
Ticks -> {Automatic, Automatic}, AxesOrigin -> {0, 0},
FrameTicks -> {{None, None}, {None, None}},
GridLines -> {None, None}, DisplayFunction -> Identity,
PlotRangePadding -> {{
Scaled[0.1],
Scaled[0.1]}, {
Scaled[0.1],
Scaled[0.1]}}, PlotRangeClipping -> True, ImagePadding ->
All, DisplayFunction -> Identity, AspectRatio -> 1,
Axes -> {False, False}, AxesLabel -> {None, None},
AxesOrigin -> {0, 0}, DisplayFunction :> Identity,
Frame -> {{True, True}, {True, True}},
FrameLabel -> {{None, None}, {None, None}}, FrameStyle ->
Directive[
Opacity[0.5],
Thickness[Tiny],
RGBColor[0.368417, 0.506779, 0.709798]],
FrameTicks -> {{None, None}, {None, None}},
GridLines -> {None, None}, GridLinesStyle -> Directive[
GrayLevel[0.5, 0.4]], ImageSize ->
Dynamic[{
Automatic, 3.5 CurrentValue["FontCapHeight"]/
AbsoluteCurrentValue[Magnification]}],
Method -> {
"DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" ->
AbsolutePointSize[6], "ScalingFunctions" -> None,
"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
(
Part[{{Identity, Identity}, {Identity, Identity}}, 1,
2][#]& )[
Part[#, 1]],
(Part[{{Identity, Identity}, {Identity, Identity}}, 2,
2][#]& )[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
(
Part[{{Identity, Identity}, {Identity, Identity}}, 1,
2][#]& )[
Part[#, 1]],
(
Part[{{Identity, Identity}, {Identity, Identity}}, 2,
2][#]& )[
Part[#, 2]]}& )}},
PlotRange -> {{0., 100.}, {0., 55.74810203315798}},
PlotRangeClipping -> True, PlotRangePadding -> {{
Scaled[0.1],
Scaled[0.1]}, {
Scaled[0.1],
Scaled[0.1]}}, Ticks -> {Automatic, Automatic}}],
GridBox[{{
RowBox[{
TagBox["\"Domain: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox[
RowBox[{"{",
RowBox[{"{",
RowBox[{"0.`", ",", "100.`"}], "}"}], "}"}],
"SummaryItem"]}]}, {
RowBox[{
TagBox["\"Output: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["\"scalar\"", "SummaryItem"]}]}},
GridBoxAlignment -> {
"Columns" -> {{Left}}, "Rows" -> {{Automatic}}},
AutoDelete -> False,
GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
GridBoxSpacings -> {
"Columns" -> {{2}}, "Rows" -> {{Automatic}}},
BaseStyle -> {
ShowStringCharacters -> False, NumberMarks -> False,
PrintPrecision -> 3, ShowSyntaxStyles -> False}]}},
GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False,
GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
BaselinePosition -> {1, 1}], True -> GridBox[{{
PaneBox[
ButtonBox[
DynamicBox[
FEPrivate`FrontEndResource[
"FEBitmaps", "SquareMinusIconMedium"]],
ButtonFunction :> (Typeset`open$$ = False), Appearance ->
None, Evaluator -> Automatic, Method -> "Preemptive"],
Alignment -> {Center, Center}, ImageSize ->
Dynamic[{
Automatic, 3.5 CurrentValue["FontCapHeight"]/
AbsoluteCurrentValue[Magnification]}]],
GraphicsBox[{{{{}, {},
TagBox[{
Directive[
Opacity[1.],
RGBColor[0.368417, 0.506779, 0.709798],
AbsoluteThickness[1]],
LineBox[CompressedData["
1:eJwBUQOu/CFib1JlAgAAADQAAAACAAAAwxW2NglTLkAIPrPOwd9LQBCpiU0A
Zy5ARY3GM30xS0BYwHE0uPYwQPK0mhYLeDhAVeKn9fSbMkAXM4rVYiAnQPGM
XlDoODRALjmry0tJFkAMZYQQ9vg1QFuGC+ExYARAVXMzhAibN0D4T5dWA9fz
Px6vUV01YDlAKo14uytp4j+Fc/DPGB07QMkmkHlHjNE/GW4Y9gC8PEDh71KI
YcjBPy2Wr4EDfj5AP2ZeGy05sT84+mdgBRFAQLSg5RkMtaE/GcCvMpb0QEBT
C6nxam6RP0nKt1EC1EFAKL8Ntk+YgT+Pb4TKcKRCQA6worEay3I/lasIdmyG
Q0Bqp9BQdkJjP7KCUXtqWURAz2Yb0g/aVD8enlrNQyhFQPAzN/urEkc/SlAb
UqoIRkAGV0TowYo4P42doDAT2kZAkGW/0vp/Kz+Pgd1BCb1HQMsI/ZXJqh0/
4anan9qbSEAii1vULlsQP0ltnFeua0lAcWeiDFTyAj9xxxVCD01KQIBpBDIb
K/U+sLxThnIfS0Aj/9247s7oPj72URex7UtAIZpQ2SGi3T6MxgfbfM1MQOL7
pgDVGNE+8DGC+EqeTUCdhCSRmKXEPhQ0tEimgE5AomZ9T2cXuD5P0aryA1RP
QBHPPEayba0+bdmwdJ4RUEAkpFJpYkqiPpIVaIkBglBAzVzuHA0Olj5CnwHL
5epQQGRNKnoXzIs+UvT2pZBcUUAC07KT3+GAPoprTCcpzFFAK4fPv0PTdD5N
MITVQjRSQE9WWaE3fGo+cMAXHSOlUkBhvNwV90pgPh+ejZGEDlNA0n/WY0xJ
VT71nWOs03VTQIn/116A3U4+K2mVYOnlU0D2YQYVRfdFPuyBqUGATlRAL/BJ
ELHoPT4NZhm83b9UQDY1xmLbnDE+Vmzp3CgvVUC0cv2nlPsmPirAmyr1llVA
6m/4iB4yIz5e36kRiAdWQJkgNlNCCys+HUyaJZxwVkBMq2/vTQs2PjyE5tJ2
4lZA9/3bNXLFPT6D3pImP1JXQDos8DmWmz4+VYYhp4i6V0Cr9sJy2mo5Pof5
C8GYK1hA4hLyj6dcLz5FutgHKpVYQEi8i5rbbRc+Phnk+P//WEBvVfy7623s
PbthcMc=
"]]},
Annotation[#,
"Charting`Private`Tag$2730#1"]& ]}}, {}, {}}, {
DisplayFunction -> Identity,
Ticks -> {Automatic, Automatic}, AxesOrigin -> {0, 0},
FrameTicks -> {{None, None}, {None, None}},
GridLines -> {None, None}, DisplayFunction -> Identity,
PlotRangePadding -> {{
Scaled[0.1],
Scaled[0.1]}, {
Scaled[0.1],
Scaled[0.1]}}, PlotRangeClipping -> True, ImagePadding ->
All, DisplayFunction -> Identity, AspectRatio -> 1,
Axes -> {False, False}, AxesLabel -> {None, None},
AxesOrigin -> {0, 0}, DisplayFunction :> Identity,
Frame -> {{True, True}, {True, True}},
FrameLabel -> {{None, None}, {None, None}}, FrameStyle ->
Directive[
Opacity[0.5],
Thickness[Tiny],
RGBColor[0.368417, 0.506779, 0.709798]],
FrameTicks -> {{None, None}, {None, None}},
GridLines -> {None, None}, GridLinesStyle -> Directive[
GrayLevel[0.5, 0.4]], ImageSize ->
Dynamic[{
Automatic, 3.5 CurrentValue["FontCapHeight"]/
AbsoluteCurrentValue[Magnification]}],
Method -> {
"DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" ->
AbsolutePointSize[6], "ScalingFunctions" -> None,
"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
(
Part[{{Identity, Identity}, {Identity, Identity}}, 1,
2][#]& )[
Part[#, 1]],
(
Part[{{Identity, Identity}, {Identity, Identity}}, 2,
2][#]& )[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
(
Part[{{Identity, Identity}, {Identity, Identity}}, 1,
2][#]& )[
Part[#, 1]],
(
Part[{{Identity, Identity}, {Identity, Identity}}, 2,
2][#]& )[
Part[#, 2]]}& )}},
PlotRange -> {{0., 100.}, {0., 55.74810203315798}},
PlotRangeClipping -> True, PlotRangePadding -> {{
Scaled[0.1],
Scaled[0.1]}, {
Scaled[0.1],
Scaled[0.1]}}, Ticks -> {Automatic, Automatic}}],
GridBox[{{
RowBox[{
TagBox["\"Domain: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox[
RowBox[{"{",
RowBox[{"{",
RowBox[{"0.`", ",", "100.`"}], "}"}], "}"}],
"SummaryItem"]}]}, {
RowBox[{
TagBox["\"Output: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["\"scalar\"", "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Order: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["3", "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Method: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["\"Hermite\"", "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Periodic: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["False", "SummaryItem"]}]}},
GridBoxAlignment -> {
"Columns" -> {{Left}}, "Rows" -> {{Automatic}}},
AutoDelete -> False,
GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
GridBoxSpacings -> {
"Columns" -> {{2}}, "Rows" -> {{Automatic}}},
BaseStyle -> {
ShowStringCharacters -> False, NumberMarks -> False,
PrintPrecision -> 3, ShowSyntaxStyles -> False}]}},
GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False,
GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
BaselinePosition -> {1, 1}]},
Dynamic[Typeset`open$$], ImageSize -> Automatic],
BaselinePosition -> Baseline], DynamicModuleValues :> {}],
StyleBox["]", "NonInterpretableSummary"]}]},
"CopyTag",
DisplayFunction->(#& ),
InterpretationFunction->("InterpolatingFunction[{{0., 100.}}, <>]"& )],
False,
Editable->False,
SelectWithContents->True,
Selectable->False]}], ",",
RowBox[{"i", "\[Rule]",
TagBox[
TemplateBox[{RowBox[{
StyleBox[
TagBox["InterpolatingFunction", "SummaryHead"],
"NonInterpretableSummary"],
StyleBox["[", "NonInterpretableSummary"],
DynamicModuleBox[{Typeset`open$$ = False},
PanelBox[
PaneSelectorBox[{False -> GridBox[{{
PaneBox[
ButtonBox[
DynamicBox[
FEPrivate`FrontEndResource[
"FEBitmaps", "SquarePlusIconMedium"]],
ButtonFunction :> (Typeset`open$$ = True), Appearance ->
None, Evaluator -> Automatic, Method -> "Preemptive"],
Alignment -> {Center, Center}, ImageSize ->
Dynamic[{
Automatic, 3.5 CurrentValue["FontCapHeight"]/
AbsoluteCurrentValue[Magnification]}]],
GraphicsBox[{{{{}, {},
TagBox[{
Directive[
Opacity[1.],
RGBColor[0.368417, 0.506779, 0.709798],
AbsoluteThickness[1]],
LineBox[CompressedData["
1:eJwVxXlUEwQAB+DBMuhhBA6CEcgRch9jjrIY7MeQG2Qb41KuECQ5BxLK6U42
YMpRwGoTKETewOJSUR48yCceqPGMQyKUh0lBwQu5k8zqj+99lknZnBRNAoFA
/8//qyTXtJTZA155uPyuBsEavo/cMmiGLxlZ3OYvJZYfYE+Ky5U3XXSwfVrh
z7wRiLR7fqkPKAZYkbovMwOPYVR4ZiZbbQLTV3MtTr1pMCzozfF6bIGx07bT
vQlnYGHmGfeDmTVI1ZO7aweFsDdtJX3ebosQZ2eP8QYpeE66r27ZOGLdz0WU
2FcBPXV5KCHVBW7DReU6OnJED+fQrm9RoLhZqW4olCPiWdvmhpKK3EXvNxZP
yvFCHs9aPEQDaWVnlnhAjr0pYr2tHXcIe0QT20OViDh1UKyt/BC79vdavKMq
sb+THeee+BFMFWWyq/MVkJxaHTM284A/t21WllkBcpbxFmGAjusy950pjQo4
13cKiKVeaE/StEg5X45DL1MT9+gDHfp3R1yty0Ee7qTWbwCsYI7VUo8MgaN9
bM0Jb/yTvnC/xleGqmzxKO8iExM3lArPeSnqOrTzM8p8EFPCGvfLl6JHyC5K
izkMTp2eu5GBFPsbjz70sfOFMY1qc/xqGYJY73Hatf3QVjRyKZZVhnMXmimk
GT+cF1rYdqxLUCV44v3WoD9cv7lD/61agnn1TNOMMgBWIzaWX7tJEJipivou
PRBGjLlMxowYOXcMme+EBiGlN1lGKhGjayC8LNc8GGHi7tUaWzG8tthG3SvB
6LKtnSZPiiC/rfVz8cMQ0LTF7SYlIqy6CXeCGkNRNR7zbaKTCAa0KwGX+EfA
W3Uxj5oSIuxR+xqRGwaVs/l6rFSIzJ+68hr1WZj2rGygeAhRa/hpuYaQhYS5
vYVrSwIwmzTuL/zOguUXv1BPNAnA2dhqrY5lw/ruCHs3TICK9BBb+wdsqErG
qAlEAUIp4kkmjQPdY8n7TAb54A0bmb3fykHc32d9KLl8mDeeqxHrhoM8frQj
y5mPzQvHL49/Fo4/mz/h0qfOwvFkXeLSQjjsgpOHBgPOYlmg9uYf4aJ6Vvfp
8u1SvH2CPPBkiItx/qyqgV6KyIylWTvXCFDUpBSbvhI4DCsdNusjwMtMmrhF
L0H8nFpfphWJfbHEbZ/+YpiFfHVAxYvEH3nREQJmMejiX3/MfxYJS2qq8MXN
Ioie1towOFH42DTf6rl/ESLKLNtk/VG4xmjQJD4uhPQvne8VjtFgXFQ0d8cX
wrd/WsuhLhrxKnGA92oBvMjPd16/jkaLROTIFRbgXxc1pEs=
"]]},
Annotation[#,
"Charting`Private`Tag$2857#1"]& ]}}, {}, {}}, {
DisplayFunction -> Identity,
Ticks -> {Automatic, Automatic}, AxesOrigin -> {0, 0},
FrameTicks -> {{None, None}, {None, None}},
GridLines -> {None, None}, DisplayFunction -> Identity,
PlotRangePadding -> {{
Scaled[0.1],
Scaled[0.1]}, {
Scaled[0.1],
Scaled[0.1]}}, PlotRangeClipping -> True, ImagePadding ->
All, DisplayFunction -> Identity, AspectRatio -> 1,
Axes -> {False, False}, AxesLabel -> {None, None},
AxesOrigin -> {0, 0}, DisplayFunction :> Identity,
Frame -> {{True, True}, {True, True}},
FrameLabel -> {{None, None}, {None, None}}, FrameStyle ->
Directive[
Opacity[0.5],
Thickness[Tiny],
RGBColor[0.368417, 0.506779, 0.709798]],
FrameTicks -> {{None, None}, {None, None}},
GridLines -> {None, None}, GridLinesStyle -> Directive[
GrayLevel[0.5, 0.4]], ImageSize ->
Dynamic[{
Automatic, 3.5 CurrentValue["FontCapHeight"]/
AbsoluteCurrentValue[Magnification]}],
Method -> {
"DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" ->
AbsolutePointSize[6], "ScalingFunctions" -> None,
"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
(
Part[{{Identity, Identity}, {Identity, Identity}}, 1,
2][#]& )[
Part[#, 1]],
(
Part[{{Identity, Identity}, {Identity, Identity}}, 2,
2][#]& )[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
(
Part[{{Identity, Identity}, {Identity, Identity}}, 1,
2][#]& )[
Part[#, 1]],
(
Part[{{Identity, Identity}, {Identity, Identity}}, 2,
2][#]& )[
Part[#, 2]]}& )}},
PlotRange -> {{0., 100.}, {0., 911.3240538519631}},
PlotRangeClipping -> True, PlotRangePadding -> {{
Scaled[0.1],
Scaled[0.1]}, {
Scaled[0.1],
Scaled[0.1]}}, Ticks -> {Automatic, Automatic}}],
GridBox[{{
RowBox[{
TagBox["\"Domain: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox[
RowBox[{"{",
RowBox[{"{",
RowBox[{"0.`", ",", "100.`"}], "}"}], "}"}],
"SummaryItem"]}]}, {
RowBox[{
TagBox["\"Output: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["\"scalar\"", "SummaryItem"]}]}},
GridBoxAlignment -> {
"Columns" -> {{Left}}, "Rows" -> {{Automatic}}},
AutoDelete -> False,
GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
GridBoxSpacings -> {
"Columns" -> {{2}}, "Rows" -> {{Automatic}}},
BaseStyle -> {
ShowStringCharacters -> False, NumberMarks -> False,
PrintPrecision -> 3, ShowSyntaxStyles -> False}]}},
GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False,
GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
BaselinePosition -> {1, 1}], True -> GridBox[{{
PaneBox[
ButtonBox[
DynamicBox[
FEPrivate`FrontEndResource[
"FEBitmaps", "SquareMinusIconMedium"]],
ButtonFunction :> (Typeset`open$$ = False), Appearance ->
None, Evaluator -> Automatic, Method -> "Preemptive"],
Alignment -> {Center, Center}, ImageSize ->
Dynamic[{
Automatic, 3.5 CurrentValue["FontCapHeight"]/
AbsoluteCurrentValue[Magnification]}]],
GraphicsBox[{{{{}, {},
TagBox[{
Directive[
Opacity[1.],
RGBColor[0.368417, 0.506779, 0.709798],
AbsoluteThickness[1]],
LineBox[CompressedData["
1:eJwVxXlUEwQAB+DBMuhhBA6CEcgRch9jjrIY7MeQG2Qb41KuECQ5BxLK6U42
YMpRwGoTKETewOJSUR48yCceqPGMQyKUh0lBwQu5k8zqj+99lknZnBRNAoFA
/8//qyTXtJTZA155uPyuBsEavo/cMmiGLxlZ3OYvJZYfYE+Ky5U3XXSwfVrh
z7wRiLR7fqkPKAZYkbovMwOPYVR4ZiZbbQLTV3MtTr1pMCzozfF6bIGx07bT
vQlnYGHmGfeDmTVI1ZO7aweFsDdtJX3ebosQZ2eP8QYpeE66r27ZOGLdz0WU
2FcBPXV5KCHVBW7DReU6OnJED+fQrm9RoLhZqW4olCPiWdvmhpKK3EXvNxZP
yvFCHs9aPEQDaWVnlnhAjr0pYr2tHXcIe0QT20OViDh1UKyt/BC79vdavKMq
sb+THeee+BFMFWWyq/MVkJxaHTM284A/t21WllkBcpbxFmGAjusy950pjQo4
13cKiKVeaE/StEg5X45DL1MT9+gDHfp3R1yty0Ee7qTWbwCsYI7VUo8MgaN9
bM0Jb/yTvnC/xleGqmzxKO8iExM3lArPeSnqOrTzM8p8EFPCGvfLl6JHyC5K
izkMTp2eu5GBFPsbjz70sfOFMY1qc/xqGYJY73Hatf3QVjRyKZZVhnMXmimk
GT+cF1rYdqxLUCV44v3WoD9cv7lD/61agnn1TNOMMgBWIzaWX7tJEJipivou
PRBGjLlMxowYOXcMme+EBiGlN1lGKhGjayC8LNc8GGHi7tUaWzG8tthG3SvB
6LKtnSZPiiC/rfVz8cMQ0LTF7SYlIqy6CXeCGkNRNR7zbaKTCAa0KwGX+EfA
W3Uxj5oSIuxR+xqRGwaVs/l6rFSIzJ+68hr1WZj2rGygeAhRa/hpuYaQhYS5
vYVrSwIwmzTuL/zOguUXv1BPNAnA2dhqrY5lw/ruCHs3TICK9BBb+wdsqErG
qAlEAUIp4kkmjQPdY8n7TAb54A0bmb3fykHc32d9KLl8mDeeqxHrhoM8frQj
y5mPzQvHL49/Fo4/mz/h0qfOwvFkXeLSQjjsgpOHBgPOYlmg9uYf4aJ6Vvfp
8u1SvH2CPPBkiItx/qyqgV6KyIylWTvXCFDUpBSbvhI4DCsdNusjwMtMmrhF
L0H8nFpfphWJfbHEbZ/+YpiFfHVAxYvEH3nREQJmMejiX3/MfxYJS2qq8MXN
Ioie1towOFH42DTf6rl/ESLKLNtk/VG4xmjQJD4uhPQvne8VjtFgXFQ0d8cX
wrd/WsuhLhrxKnGA92oBvMjPd16/jkaLROTIFRbgXxc1pEs=
"]]},
Annotation[#,
"Charting`Private`Tag$2857#1"]& ]}}, {}, {}}, {
DisplayFunction -> Identity,
Ticks -> {Automatic, Automatic}, AxesOrigin -> {0, 0},
FrameTicks -> {{None, None}, {None, None}},
GridLines -> {None, None}, DisplayFunction -> Identity,
PlotRangePadding -> {{
Scaled[0.1],
Scaled[0.1]}, {
Scaled[0.1],
Scaled[0.1]}}, PlotRangeClipping -> True, ImagePadding ->
All, DisplayFunction -> Identity, AspectRatio -> 1,
Axes -> {False, False}, AxesLabel -> {None, None},
AxesOrigin -> {0, 0}, DisplayFunction :> Identity,
Frame -> {{True, True}, {True, True}},
FrameLabel -> {{None, None}, {None, None}}, FrameStyle ->
Directive[
Opacity[0.5],
Thickness[Tiny],
RGBColor[0.368417, 0.506779, 0.709798]],
FrameTicks -> {{None, None}, {None, None}},
GridLines -> {None, None}, GridLinesStyle -> Directive[
GrayLevel[0.5, 0.4]], ImageSize ->
Dynamic[{
Automatic, 3.5 CurrentValue["FontCapHeight"]/
AbsoluteCurrentValue[Magnification]}],
Method -> {
"DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" ->
AbsolutePointSize[6], "ScalingFunctions" -> None,
"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
(
Part[{{Identity, Identity}, {Identity, Identity}}, 1,
2][#]& )[
Part[#, 1]],
(
Part[{{Identity, Identity}, {Identity, Identity}}, 2,
2][#]& )[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
(
Part[{{Identity, Identity}, {Identity, Identity}}, 1,
2][#]& )[
Part[#, 1]],
(
Part[{{Identity, Identity}, {Identity, Identity}}, 2,
2][#]& )[
Part[#, 2]]}& )}},
PlotRange -> {{0., 100.}, {0., 911.3240538519631}},
PlotRangeClipping -> True, PlotRangePadding -> {{
Scaled[0.1],
Scaled[0.1]}, {
Scaled[0.1],
Scaled[0.1]}}, Ticks -> {Automatic, Automatic}}],
GridBox[{{
RowBox[{
TagBox["\"Domain: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox[
RowBox[{"{",
RowBox[{"{",
RowBox[{"0.`", ",", "100.`"}], "}"}], "}"}],
"SummaryItem"]}]}, {
RowBox[{
TagBox["\"Output: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["\"scalar\"", "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Order: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["3", "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Method: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["\"Hermite\"", "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Periodic: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["False", "SummaryItem"]}]}},
GridBoxAlignment -> {
"Columns" -> {{Left}}, "Rows" -> {{Automatic}}},
AutoDelete -> False,
GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
GridBoxSpacings -> {
"Columns" -> {{2}}, "Rows" -> {{Automatic}}},
BaseStyle -> {
ShowStringCharacters -> False, NumberMarks -> False,
PrintPrecision -> 3, ShowSyntaxStyles -> False}]}},
GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False,
GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
BaselinePosition -> {1, 1}]},
Dynamic[Typeset`open$$], ImageSize -> Automatic],
BaselinePosition -> Baseline], DynamicModuleValues :> {}],
StyleBox["]", "NonInterpretableSummary"]}]},
"CopyTag",
DisplayFunction->(#& ),
InterpretationFunction->("InterpolatingFunction[{{0., 100.}}, <>]"& )],
False,
Editable->False,
SelectWithContents->True,
Selectable->False]}], ",",
RowBox[{"r", "\[Rule]",
TagBox[
TemplateBox[{RowBox[{
StyleBox[
TagBox["InterpolatingFunction", "SummaryHead"],
"NonInterpretableSummary"],
StyleBox["[", "NonInterpretableSummary"],
DynamicModuleBox[{Typeset`open$$ = False},
PanelBox[
PaneSelectorBox[{False -> GridBox[{{
PaneBox[
ButtonBox[
DynamicBox[
FEPrivate`FrontEndResource[
"FEBitmaps", "SquarePlusIconMedium"]],
ButtonFunction :> (Typeset`open$$ = True), Appearance ->
None, Evaluator -> Automatic, Method -> "Preemptive"],
Alignment -> {Center, Center}, ImageSize ->
Dynamic[{
Automatic, 3.5 CurrentValue["FontCapHeight"]/
AbsoluteCurrentValue[Magnification]}]],
GraphicsBox[{{{{}, {},
TagBox[{
Directive[
Opacity[1.],
RGBColor[0.368417, 0.506779, 0.709798],
AbsoluteThickness[1]],
LineBox[CompressedData["
1:eJwVxXs8FAYAB3CPCp8bUm4xneSRV5PjlD7p7udyCqfcccij7mPJ2+G8wnF3
Lo+WqO1knptP6GjkUZrxcfuoad5qJs/m1ZZq68GRfcq2P76f794QHjtURUlJ
yek//19+8a5aGa+TqvROs6W5O4fKGCVHU4gbtPs+f9zVWhymbQ21adtmQ8Cj
c9UenIYXtMhfXMMGbHXx9/tv2HNEZfSJUyd5ss/QtjKY09e4A8QLrfHUcSMY
fm/b4t5jAiPS0eAhkinaiRop9AkyLHfX7Pyq3hx7BOnjhSNHEbdf68P9fdZY
zq2dd1lnYLss31MpzAYLPAGlj+QJf3k85Z7CFrrUf17RCtjgzNetrpTZwTiZ
lt3K9MWby2e8/nSkQLEgCGS+O41PQiXbFesOmNOeHKDNBIPDt5eolx3Cl5Y2
kR4/cmHYxAp24B4G0dyqRV4bgov818N6pCOYe1YT8anoHPRj9RRKnU5Ie9i7
zY5yHp8XN4lUM6m4oeEd4LIlHI4bYdytOoCugWGJWWME9OVNdsUrgJRxVnWE
HwW3vnaWyq/OeJtizUuwjEEhT9IXd4OOojuBaosdsZA2qCdH5xzDxrTJvXD7
OLSIWemRp13w4veAswNrcTCsDBg8ZsFAhp7r6ME78XD3MmDXq7uCvvzQkpGV
gIKKb213TroixLR77KYjH4WiGWeNruPgSh9okz/yMSebrJosO4HsDkZGV2ci
3GLK/Rqj3DAyEfRRGJWE+F4iXdvTHbszCZNPzJJxu9M7J2GPBwyFDjumR5JB
VbB2Nb/yACz4xweSUnD5Z7WpjEEmYjqaa/vMUvGaLF53r/TEQk17NbU3FbqU
thO1wpPoumkkZ/Ev4NRo/VtVn1NwO/SbnqtOGmImbidW6njBi+w3UilLwzVi
eL6y2AvGGlcIBsx00KuU+5eWvfAdTbnUfi4d7BVFTVEQC4Rb7CANYQYuRTHN
LQdYuE5mljrrC+BpKxmjU9iYIlhN9dcIECffRTKpYWN2bi216kgm9lQWXJVo
eeOKGkNm/CATqxVf3Hqc5I3E7NKGHrcsWEdIuc+XvKGiKAF9PAsvRTJn4Ukf
ZFSvFeXvF0LzvH7nTLcPfrC7Sm6LF8I3+vm0xQEOsoeTBqPbhLCSl1mtFnOQ
lzEkrXgjxJmnMp08NV88XhJsmbUTgcQsNSuP88VhojenO1YEJ8mzR8nzvhhM
7afmNoqQPXttH43tB+0PPX9pvhSBk7O3Lq/DD+Nft2ZmmoqR+57wU4m1P1Js
Ww/UccVgdDxRs5L6Y1PcFfj0uhhU/cX1zU1/kA5Sl8eGxPgXFOSqaA==
"]]},
Annotation[#,
"Charting`Private`Tag$2901#1"]& ]}}, {}, {}}, {
DisplayFunction -> Identity,
Ticks -> {Automatic, Automatic}, AxesOrigin -> {0, 0},
FrameTicks -> {{None, None}, {None, None}},
GridLines -> {None, None}, DisplayFunction -> Identity,
PlotRangePadding -> {{
Scaled[0.1],
Scaled[0.1]}, {
Scaled[0.1],
Scaled[0.1]}}, PlotRangeClipping -> True, ImagePadding ->
All, DisplayFunction -> Identity, AspectRatio -> 1,
Axes -> {False, False}, AxesLabel -> {None, None},
AxesOrigin -> {0, 0}, DisplayFunction :> Identity,
Frame -> {{True, True}, {True, True}},
FrameLabel -> {{None, None}, {None, None}}, FrameStyle ->
Directive[
Opacity[0.5],
Thickness[Tiny],
RGBColor[0.368417, 0.506779, 0.709798]],
FrameTicks -> {{None, None}, {None, None}},
GridLines -> {None, None}, GridLinesStyle -> Directive[
GrayLevel[0.5, 0.4]], ImageSize ->
Dynamic[{
Automatic, 3.5 CurrentValue["FontCapHeight"]/
AbsoluteCurrentValue[Magnification]}],
Method -> {
"DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" ->
AbsolutePointSize[6], "ScalingFunctions" -> None,
"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
(
Part[{{Identity, Identity}, {Identity, Identity}}, 1,
2][#]& )[
Part[#, 1]],
(
Part[{{Identity, Identity}, {Identity, Identity}}, 2,
2][#]& )[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
(
Part[{{Identity, Identity}, {Identity, Identity}}, 1,
2][#]& )[
Part[#, 1]],
(
Part[{{Identity, Identity}, {Identity, Identity}}, 2,
2][#]& )[
Part[#, 2]]}& )}},
PlotRange -> {{0., 100.}, {0., 601.854450689341}},
PlotRangeClipping -> True, PlotRangePadding -> {{
Scaled[0.1],
Scaled[0.1]}, {
Scaled[0.1],
Scaled[0.1]}}, Ticks -> {Automatic, Automatic}}],
GridBox[{{
RowBox[{
TagBox["\"Domain: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox[
RowBox[{"{",
RowBox[{"{",
RowBox[{"0.`", ",", "100.`"}], "}"}], "}"}],
"SummaryItem"]}]}, {
RowBox[{
TagBox["\"Output: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["\"scalar\"", "SummaryItem"]}]}},
GridBoxAlignment -> {
"Columns" -> {{Left}}, "Rows" -> {{Automatic}}},
AutoDelete -> False,
GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
GridBoxSpacings -> {
"Columns" -> {{2}}, "Rows" -> {{Automatic}}},
BaseStyle -> {
ShowStringCharacters -> False, NumberMarks -> False,
PrintPrecision -> 3, ShowSyntaxStyles -> False}]}},
GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False,
GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
BaselinePosition -> {1, 1}], True -> GridBox[{{
PaneBox[
ButtonBox[
DynamicBox[
FEPrivate`FrontEndResource[
"FEBitmaps", "SquareMinusIconMedium"]],
ButtonFunction :> (Typeset`open$$ = False), Appearance ->
None, Evaluator -> Automatic, Method -> "Preemptive"],
Alignment -> {Center, Center}, ImageSize ->
Dynamic[{
Automatic, 3.5 CurrentValue["FontCapHeight"]/
AbsoluteCurrentValue[Magnification]}]],
GraphicsBox[{{{{}, {},
TagBox[{
Directive[
Opacity[1.],
RGBColor[0.368417, 0.506779, 0.709798],
AbsoluteThickness[1]],
LineBox[CompressedData["
1:eJwVxXs8FAYAB3CPCp8bUm4xneSRV5PjlD7p7udyCqfcccij7mPJ2+G8wnF3
Lo+WqO1knptP6GjkUZrxcfuoad5qJs/m1ZZq68GRfcq2P76f794QHjtURUlJ
yek//19+8a5aGa+TqvROs6W5O4fKGCVHU4gbtPs+f9zVWhymbQ21adtmQ8Cj
c9UenIYXtMhfXMMGbHXx9/tv2HNEZfSJUyd5ss/QtjKY09e4A8QLrfHUcSMY
fm/b4t5jAiPS0eAhkinaiRop9AkyLHfX7Pyq3hx7BOnjhSNHEbdf68P9fdZY
zq2dd1lnYLss31MpzAYLPAGlj+QJf3k85Z7CFrrUf17RCtjgzNetrpTZwTiZ
lt3K9MWby2e8/nSkQLEgCGS+O41PQiXbFesOmNOeHKDNBIPDt5eolx3Cl5Y2
kR4/cmHYxAp24B4G0dyqRV4bgov818N6pCOYe1YT8anoHPRj9RRKnU5Ie9i7
zY5yHp8XN4lUM6m4oeEd4LIlHI4bYdytOoCugWGJWWME9OVNdsUrgJRxVnWE
HwW3vnaWyq/OeJtizUuwjEEhT9IXd4OOojuBaosdsZA2qCdH5xzDxrTJvXD7
OLSIWemRp13w4veAswNrcTCsDBg8ZsFAhp7r6ME78XD3MmDXq7uCvvzQkpGV
gIKKb213TroixLR77KYjH4WiGWeNruPgSh9okz/yMSebrJosO4HsDkZGV2ci
3GLK/Rqj3DAyEfRRGJWE+F4iXdvTHbszCZNPzJJxu9M7J2GPBwyFDjumR5JB
VbB2Nb/yACz4xweSUnD5Z7WpjEEmYjqaa/vMUvGaLF53r/TEQk17NbU3FbqU
thO1wpPoumkkZ/Ev4NRo/VtVn1NwO/SbnqtOGmImbidW6njBi+w3UilLwzVi
eL6y2AvGGlcIBsx00KuU+5eWvfAdTbnUfi4d7BVFTVEQC4Rb7CANYQYuRTHN
LQdYuE5mljrrC+BpKxmjU9iYIlhN9dcIECffRTKpYWN2bi216kgm9lQWXJVo
eeOKGkNm/CATqxVf3Hqc5I3E7NKGHrcsWEdIuc+XvKGiKAF9PAsvRTJn4Ukf
ZFSvFeXvF0LzvH7nTLcPfrC7Sm6LF8I3+vm0xQEOsoeTBqPbhLCSl1mtFnOQ
lzEkrXgjxJmnMp08NV88XhJsmbUTgcQsNSuP88VhojenO1YEJ8mzR8nzvhhM
7afmNoqQPXttH43tB+0PPX9pvhSBk7O3Lq/DD+Nft2ZmmoqR+57wU4m1P1Js
Ww/UccVgdDxRs5L6Y1PcFfj0uhhU/cX1zU1/kA5Sl8eGxPgXFOSqaA==
"]]},
Annotation[#,
"Charting`Private`Tag$2901#1"]& ]}}, {}, {}}, {
DisplayFunction -> Identity,
Ticks -> {Automatic, Automatic}, AxesOrigin -> {0, 0},
FrameTicks -> {{None, None}, {None, None}},
GridLines -> {None, None}, DisplayFunction -> Identity,
PlotRangePadding -> {{
Scaled[0.1],
Scaled[0.1]}, {
Scaled[0.1],
Scaled[0.1]}}, PlotRangeClipping -> True, ImagePadding ->
All, DisplayFunction -> Identity, AspectRatio -> 1,
Axes -> {False, False}, AxesLabel -> {None, None},
AxesOrigin -> {0, 0}, DisplayFunction :> Identity,
Frame -> {{True, True}, {True, True}},
FrameLabel -> {{None, None}, {None, None}}, FrameStyle ->
Directive[
Opacity[0.5],
Thickness[Tiny],
RGBColor[0.368417, 0.506779, 0.709798]],
FrameTicks -> {{None, None}, {None, None}},
GridLines -> {None, None}, GridLinesStyle -> Directive[
GrayLevel[0.5, 0.4]], ImageSize ->
Dynamic[{
Automatic, 3.5 CurrentValue["FontCapHeight"]/
AbsoluteCurrentValue[Magnification]}],
Method -> {
"DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" ->
AbsolutePointSize[6], "ScalingFunctions" -> None,
"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
(
Part[{{Identity, Identity}, {Identity, Identity}}, 1,
2][#]& )[
Part[#, 1]],
(
Part[{{Identity, Identity}, {Identity, Identity}}, 2,
2][#]& )[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
(
Part[{{Identity, Identity}, {Identity, Identity}}, 1,
2][#]& )[
Part[#, 1]],
(
Part[{{Identity, Identity}, {Identity, Identity}}, 2,
2][#]& )[
Part[#, 2]]}& )}},
PlotRange -> {{0., 100.}, {0., 601.854450689341}},
PlotRangeClipping -> True, PlotRangePadding -> {{
Scaled[0.1],
Scaled[0.1]}, {
Scaled[0.1],
Scaled[0.1]}}, Ticks -> {Automatic, Automatic}}],
GridBox[{{
RowBox[{
TagBox["\"Domain: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox[
RowBox[{"{",
RowBox[{"{",
RowBox[{"0.`", ",", "100.`"}], "}"}], "}"}],
"SummaryItem"]}]}, {
RowBox[{
TagBox["\"Output: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["\"scalar\"", "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Order: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["3", "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Method: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["\"Hermite\"", "SummaryItem"]}]}, {
RowBox[{
TagBox["\"Periodic: \"", "SummaryItemAnnotation"],
"\[InvisibleSpace]",
TagBox["False", "SummaryItem"]}]}},
GridBoxAlignment -> {
"Columns" -> {{Left}}, "Rows" -> {{Automatic}}},
AutoDelete -> False,
GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
GridBoxSpacings -> {
"Columns" -> {{2}}, "Rows" -> {{Automatic}}},
BaseStyle -> {
ShowStringCharacters -> False, NumberMarks -> False,
PrintPrecision -> 3, ShowSyntaxStyles -> False}]}},
GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False,
GridBoxItemSize -> {
"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
BaselinePosition -> {1, 1}]},
Dynamic[Typeset`open$$], ImageSize -> Automatic],
BaselinePosition -> Baseline], DynamicModuleValues :> {}],
StyleBox["]", "NonInterpretableSummary"]}]},
"CopyTag",
DisplayFunction->(#& ),
InterpretationFunction->("InterpolatingFunction[{{0., 100.}}, <>]"& )],
False,
Editable->False,
SelectWithContents->True,
Selectable->False]}]}], "}"}], "}"}]], "Output",
CellChangeTimes->{3.6980865870613832`*^9}],
Cell[BoxData[
GraphicsBox[{{{}, {},
TagBox[
{RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
1.], LineBox[CompressedData["
1:eJwVWnc8l98XRxRRX3uFhOy9573H3vuziCIZWVnJ3kWyQoNESRkJyQiRZBYS
kqKioUKESBI/v7+e1/v1PPec93nfc8499/V6Dh0/ZedOQ0VF9W4XFdX/n0KZ
/24VCjxBj/zubGzjTGiklzjo51WEb8lEZD9FmfAkrKpE3Ksab1lZ+SXoZEK3
fKDnbq/HOM+pr3hJMxN604ZNVj1fYG7eYaZclUz4mcT0NaziJf7Bz9skp5wJ
0U+rJX7fGMFcX+NYuxUz4TKsVa2kv8bPz9BEf5TLhD7NxKZ5zw/4PXXNymuJ
TLirzRajUj2F9W5RBozEM2Fp4kRRSMU0dirWvFctmgkxvPQfl298wjJ5e8RP
CWfClStWnovpX7F7hm1yMV8mPLs4EfjNcxEbJ23Ode3LhJSrl/tyFH7iXYW5
B54xZoK0F3ZUqv6JCQHpfd0MmRC052JkUMUSXqfxVSiny4QtQ+W2nzdW8PyP
CN7//mUA99MI/R/pv7HzocM+hrMZ0JItPJrGuI4fKF9/+fRrBhw93n9C+vw6
Dmn+Maz+JQNKqAXO+ST9wWq63od2T2WAAjzp/R7xFzMyzrJJvsoAs9Y9VjOe
2/i21GrDmdad99vHVzl/buNX8tVn9zZnQJrgsUFhbSqw5le/mNmQASHHm08T
OqngmeCtlFPVGWAwE9BZP0oN9tUHuONvZIB6sgD3IDcNqCeKzRcVZIC0eL/P
jBMNaHX7m1RfzQA2bzE2ri80YPVFvaYwKwM+/njnGra6C16PaI01xWbAWMaF
+iwNWjj+5nI0KTIDnslpMJRH00LFMMvAdGgG1Abm1ryhowMuh1bPp34ZEL9q
tq3JuRtkF33T7zpmgOC/xoJ/qvQwoC4ZMCqfAeyF7j85Iulh7rZivLZUBtBj
NgPZx/RARXm2nHE4Axbj/OeOGjPAFw6moSWeDGijPazZTt4LBSb+dKXUGeDM
lP06IYwJno6Uq6f2p4PNPSx17RETtOZsNwZ3pYOB1Y+YB1T7wONd9yW9tnSQ
zjIR/ZyyD9Qz8rgyq9Nhk207xDBvP2iUxdx6k5UOA/tUGWidmUFlvT5lxjId
hs1Zrw5GMQN63ncz1zAdxs8viOYVMIMN+8B1UZ10+ExXaiA7yQwEkwrGVakd
e/+44ihHWOB04mMTJfp0kFz4s1bpwApRbbrVP5vTQF567OyZcFYwtE3xf1iT
Bqretex6eayQt/+KxLE7aaA34634epwV+BxKiSYX08Dh/YQfNYUN6m2dx+Pc
0yBlsPUzkcQOc7EBlYZ70iCDKT9YMJQdOk+sqYxvXIBcs1CauUvsYKrC/8Zg
4QIUdcsKxr1ih/atN1p3Ry9AQ1vRkQoCB7AeqnMrKboAX6rih7fsOKHEsugx
g+wFmJ13du0L4oRTAve2Dhy8AD8lNX/mZHNC4kj6713MF2DzzvI+iWFO8Gn1
FycspQJ7kZupvS0XdNYU7V68nwp6mYbtpdbcMFtg3vdBMhXizOrKh/25QezS
o79pPKnQRiec8y+dG0433kvcT58KmlFUHvb93HBKz2xXxufzoHSyhXHbhAd2
GZgERF47DyL6CmSSHi88CRC8+o/qPOz+zbdAq8IHXQcPWCZVJINB7YXXcgQ+
iNUnJN/OSYYEv412x2A+YKz66HctKhm2Pr3OqbnPB2EMrHN0Vsnweyhbw0mW
H3LvRCtMLpyDbxX0Zx+ICcBx4e71Z+LnwF1uY5/CPkHorJW8ZBKXBBYTybe5
BQTBc9dhL6PjSaCczKFDJScIu5UvqQsaJAHtB3nfFzaCUF4ZNOO9JwmKMzyf
+eUKwo8O370haYnw4cfouYoDh6Dn0WoD98UEcKis3haWFIJHKyEnwy7EAVB0
rjBqCUGPQ4DDY+84EKd9LrtiLgThq9dvU5vFwfqRGecOPyFQUuztiKKPg8uM
/I9c7gvB5MxNm/G9sTDinRp2XV0Ycil/JqyposFSwn2J01gEyl/YTqtSh8PG
hdMt70giIM4+zvO7OwzKFs6eLfEQAf2N24Hf08KAtu4Oj+I5EeAdrQn25gmD
ZvRN17JLBI67lvjRqpwBUYJPdqL+Ycgof2HXcuY0UMUGKv3EopAXaq62oRgE
VdPx/xqsRYGYnnTRhjoIjhhk90QfEwX60cT6jReB0MDw4AhjrCi8Lzd+ZeEf
CL65vxJF20QhKlVEaq0qAN6Unxl10haD5tXE5rh4f6gdjT7dpy4O9mpHliet
T0LZ2RPHP5qKgw6PjdnRTS8oVDW3/usoDnnP6A/XlXtB6lVuCZlocbiytY80
sdsL3JweTGY9EYcj+wyH6Z97APvnb/oUMwn4YkFstoo4AXsvvZAPPCIBHK8+
XR9QOgFURg38qb4SYGkT48Sx4AbzpYm/WzIkIP4Io8TfE27Q5SNw9+CIBDiO
BfyldjoOoSv2rF+PSIKdvFDykp8L+JZobm/7SgLdL/w8WMkFjhMPzXPHSMKr
E/4vSSvHwKpxocvshiQcYCW6eDMfA7HI8+FVnyXh1FdTH64jzjBO83jqtJ8U
TNsO927KOcLgg9sDGTFSILm/qpMw5LBT12nNpZlSEFgSNfwt0AHudzvkvrkv
BSd+e3M3PaRAauovY501Kfh6xSnmNoEMWqyS1bSx0rD3y51P9E8JMEjt2DeQ
KQ3b508YCZwhgMvS+U+Xb0iDUWDLVxVpAiQNfeeS6JCG9p9H+lzy7GEwozzO
kk4GbqlFS1Fi7cCFUcLu8gUZoN/Q5fc+YwPLGxTfYwUycC638O4RTRtImk05
J35PBhI271WZbFpDed+35uZBGajC+466MlvDcnKZ8AcWWag7tYfmg7klJNGK
r4pdlYWBRXe+x4xmwPmL/N9SmSyEbfyZUqowhbJPyRLNTbLwurbNgMHMFAY6
vjpbTMiCd64FdKabAGd8aXeAgByw8LKMUh02hvJ/olebSuTgcbNA6HiGAWj/
INUm1MvBuV+ZvgRdAxicPNdv3i0HLY0H2RZ+6cNyywzV+69ykML6McrnmD4I
6Iwox3PIQ6r7AK+orh54/s2+fpsiD9WsA7dDlAFyreafhrvIw/MLo09d1jC0
3zSctfSShzd6612eZzFwm6yrrp2RBy+Ff9TMZQh6c52HjK7IQxNv713NbW2Q
kBXf9e2VPPRG5iu1/NQAUny8xKN38sB9yWdmMlcDEkbfWmd9kYf/ohro9DQ1
4G1EeoH6qjxI9P0u7kpRh9SeZZXz7ArAm/H2kLayGsy5tHpJ2ivAy+Xjoe5P
lIG7jjNzy1EBavv+9r86owwGewLqh48rgGivabuZnDIU3BOiiQxSgHPtvsXK
N5XAYiP52vOLCjArXHmqNEsR7uXYDfoOKcAPC5bHUtfk4e3M3V8wrgAb9znM
SDtx79akO8AxpQD7LBi47ETl4ehUo2frogI0RSffOHFSDvbL8FPv/08RWDcc
3DPqZeBU91elaktFaFidyX+VIwntUOkjTVKEFMb2ojkZSWBuCbhVflQRbkr5
WVzok4Ca6j+sJacUwXbCOs+GTgKWrzIuX72oCBeCbgraZYiBHtuQBGe+Igi3
FOhyKItBTnqua06xInhLd1uqTIiCcgL/y/QHiuD76j/CGRlRCPWRq0l4pQgh
LPfqgqZFoPvLytftd4pQMlL+JP2iCHC5PDwYPaMIdRtn37zTE4GHBN3MM78V
QX1Lkri3XBj+6Nj7+/IoQUiGTg17ohBEMZ+RJjkrAa+o4lM1b0HoT9U6MeKu
BKGcntOXpQSBn46qwMZfCS7NP116MXEQ2v6kMJrHKgFjuE5+8awAUH3Kn8U7
ut+ki3rSxs8PNs7HhFrLlYBbaFhjrIoPbr4WdtSsVYL2KMWOPF0+0Ouv7FN+
qgQBKmic0fcAJNa3lUl8UYIfVukBCyM8MCqXOHXnhxJ0lhtvNAbygEiFMbfI
mhJ0/LvjusnMA52FQ8kCe5ThShFLygd7btid8tGDTUIZggJbQoJnOSHVcffh
f77K8NLFYVhxZ45A8kxr33cSvZR+wVpnPzss0bH2jEUpg4pTSbVZDxtQ7vOf
rElT3pl31pKPabPBYXqVe26VynD+3hYwKbJCe72byrM5ZbjrHoBiTZgh5MLJ
3Y0rymALxFye/5hBzPXU61t/lWFEXKp7173/duaVyPBoRhX4JJF3xvXFfjji
lt0qL6UCnD+cC31F98Ea8xPDK94q8Of85AwPz16omOnmTApSgdcRWhPhrQzg
/Kj/a0CECkRGVwcfcmOATs/x82apKnCNjeqiSi09XGxbHPhXrgJ+dO4xK0f3
gLSvAOnEdxXAK0unqRZpYUpXRMx2SQUGP7yqfnmbFnK4JNd1/qjszMFBKaXO
tLDRoZLPxaAK1Ckvue2GdkEvj+X7Z+Kq0Kr/PlDvCQ249UR6KHipwtu3ckp5
g1RwQlvUau2UKsz87efRuUQFHveHVFrOqIL94KXhMicq8Lp2mM4wWRU6aqTp
8ia2sf+pF7cod1Th2i0b7UWFLRzBJTwd+1kVnryJqSIe+ouj0gZ6DeZVITaM
xexw1waOpg6rYfilCrb/Bacre2/guLn+2JxdahDoT0NN/fAPPvc4VKBUSA1Y
Vudt3rqu42zPZ0cGXdRA0S6f7ffsKs6dDNHP8VKDBbdLHsM3VvEl24NSlAA1
6PnR1HyPvIqvaoVsTMeqgSP5bhWx9xcu/E8gb7VQDabtKn/sq1/BFY2BY3zv
1MDorik7oWkJV0rztU1/VoN7zFzjt6KX8L2b3bfvzKsBm0D2UrveEq65cOC0
/KYatFzcMqG8+IkbjnWxGRxQB409zruOfVrET/bw2Po4qIMF1gt0vT+PZ27J
mr1xVQd2kzEqDtN5zAQG+sbe6tD7JVnP4uMcpoSdUhGOVIf1rHxjWZ45vPit
i2eiQB3asxKyD179jvmfBX00m1IHx47cn+pjM1jPI2Xi4Td16Fa0GyxKnMFe
NIWjokvqMPqpvU1NcQbXafZ109BowCtj8effs79g87sH7zYLa0DWmKPmB9fP
ODytP1jSUwOKzt+Z8FT8iAvFP/pePaUBp/e+jZn/OY07O3+77w7TgMtT8uh4
zTRm/idE+ZiiAToh85msitO41C9cO79CAx68y0bCMIVfWYnS7V3QgNF+yaYo
+nd4Y1Zr68yaBhQrBymk3pjEgsm2v79sacBCadQkm+Yk9n0c9f3Jf5qQ6SfQ
eS5wAu+SHxkIV9CEwWunKshLb7ACS9zl7yGawCFcxC8k/xo3078W74nShOhs
0dg9n8awHpVsS0mSJlRdUHNYuTyGCYsTH47lasJrk/vfX1CP4dABNYmxB5pg
wESaUfs0ih+dX2zpWNaEy997OJZeDmPDeCOrGxua0Jpkl12eM4wHw65PRdNo
gd+J9+eUScP4g6c5nQarFny5brGr8d1LTGNUalWtoAXlDCb6ImtD2ITm2PS1
AC0oIt1buZU7iF/+qQ8OD9OCo0o/iuflBrHjEtNucpwWeBsYrqf0D2CfqWZJ
1iwtGO7KaqxgGMAZbZwhKdVaIH95PIwj5zkejRjcfXpBC6I25saevO/FzkGH
8+zWtCB18KEoTu/FMyejpOS3tOAU8pEa1O7FvymSNnP7tGE77p4ra1EP5lU7
l+cqow3+x8Otn5/qxi4rOtJWvtrwn/SBZBbcifexUg2NhmgDR75wZfLvp7hZ
viP4SLQ2xAtMhFLdf4rZ/Y1avDK0QTXjlPCs6FPc+83KLLFGG1IiZks0BTqw
/PtjXg9/acNbWZtDn9Xa8bvNQ0z4nzZw8xjaPP33GKce+FzdRacDTlZEV0r6
Y/yZ4vV7mFMHDvJJyWveb8NXRwLO/VDXgd9/h5YwYyum6osvEYrWAYekhFcJ
u5rxva/6JuVndUAZVTGmlzVhx9275+UydMDi8Tf9d1ZNuE4/VUmnSAfqkgX+
aBc9xF5t2R3kJzqgl/SDGGrViIce3JpKo0NQUXRQk/1THY4edk9i24/A7PGR
h/75dVhySUw8nxPBSl3A8HO7OpwkuzMXiCGo3G4zjOl+gNXL67aemCLgirpg
59ZUi28WdvH9TkcwYJVSbPqsBhc4Sb8mXkZQetffGl+swVd5cy7WFSJoDzis