resonans5.nb 45.3 KB
Newer Older
himyss's avatar
himyss committed
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
(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 13.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     46202,        920]
NotebookOptionsPosition[     44928,        893]
NotebookOutlinePosition[     45326,        909]
CellTagsIndexPosition[     45283,        906]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "Input", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"kin", " ", "=", " ", "20."}], ";"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"mass", "=", " ", "936."}], ";"}]}]], "Input",
 CellChangeTimes->{{3.8698997367340517`*^9, 3.869899843251645*^9}, {
   3.869900162194075*^9, 3.86990022514108*^9}, {3.869900263296324*^9, 
   3.869900278820773*^9}, {3.869900716799706*^9, 3.869900721797647*^9}, 
   3.869902565634445*^9, {3.869902624866561*^9, 3.869902626057645*^9}, {
   3.869904698619644*^9, 3.8699047019298077`*^9}, 3.870142822697193*^9, {
   3.870144764448131*^9, 3.87014476799769*^9}, {3.870144888379427*^9, 
   3.870144970243932*^9}, {3.870145792817855*^9, 3.8701457939369917`*^9}, {
   3.8713910972634153`*^9, 3.871391103619568*^9}, {3.8714518138981028`*^9, 
   3.871451832796793*^9}, 3.8714523038600388`*^9, {3.87145238207703*^9, 
   3.871452385670923*^9}, 3.8720846338697147`*^9, {3.8735416074490557`*^9, 
   3.873541607808901*^9}},
 CellLabel->"In[7]:=",ExpressionUUID->"4a74809f-4d6f-4bf5-84df-1c81e845ccce"],

Cell[BoxData[
 RowBox[{
  RowBox[{"sps", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Dashing", "[", "None", "]"}], ",", 
       RowBox[{"RGBColor", "[", 
        RowBox[{"0", ",", "0", ",", "0"}], "]"}]}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Dashing", "[", "None", "]"}], ",", 
       RowBox[{"RGBColor", "[", 
        RowBox[{"1", ",", "0", ",", "0"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Dashing", "[", "None", "]"}], ",", 
       RowBox[{"RGBColor", "[", 
        RowBox[{"0", ",", "0", ",", "1"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Dashing", "[", "None", "]"}], ",", 
       RowBox[{"RGBColor", "[", 
        RowBox[{"0", ",", "1", ",", "0"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Dashing", "[", "None", "]"}], ",", 
       RowBox[{"RGBColor", "[", 
        RowBox[{"0", ",", "1", ",", "1"}], "]"}]}], "}"}]}], "}"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{
  3.871701793434409*^9, {3.872065693342557*^9, 3.872065697733029*^9}},
 CellLabel->"In[8]:=",ExpressionUUID->"4abccf8a-d47b-424b-acb1-36e9461ec53e"],

Cell[BoxData[
 RowBox[{
  RowBox[{"kIn", "[", "e_", "]"}], " ", ":=", " ", 
  FractionBox[
   SqrtBox[
    RowBox[{"2", "*", "mass", "*", "e"}]], "197.327"]}]], "Input",
 CellChangeTimes->{
  3.870144988515087*^9, {3.870145025281687*^9, 3.870145047838258*^9}, {
   3.871451965683363*^9, 3.87145196928242*^9}, {3.872063208516612*^9, 
   3.8720632133894377`*^9}, {3.87354165429843*^9, 3.873541682028607*^9}, {
   3.8735428967851543`*^9, 3.873542921291761*^9}, {3.873729554724765*^9, 
   3.873729555224043*^9}},
 CellLabel->"In[9]:=",ExpressionUUID->"ddda9570-52bb-4f7f-8d3d-86bdfd8586a4"],

Cell[BoxData[
 RowBox[{
  RowBox[{"kOut", "[", 
   RowBox[{"e_", ",", "U0_"}], "]"}], " ", ":=", " ", 
  FractionBox[
   SqrtBox[
    RowBox[{"2", " ", "mass", " ", 
     RowBox[{"(", 
      RowBox[{"e", "-", " ", "U0"}], ")"}]}]], "197.327"]}]], "Input",
 CellChangeTimes->{{3.870145798453801*^9, 3.8701458265877657`*^9}, {
  3.871451836727564*^9, 3.871451840179946*^9}, {3.8714519734346457`*^9, 
  3.8714519739887877`*^9}, {3.872063215936698*^9, 3.8720632175080347`*^9}, {
  3.873541677338694*^9, 3.8735416793910713`*^9}, {3.873542899115011*^9, 
  3.8735429242501287`*^9}},
 CellLabel->"In[10]:=",ExpressionUUID->"86a69de3-cbb6-4ee5-8cee-9244eebdf060"],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Delta]", "[", 
   RowBox[{"e_", ",", "r_", ",", "u0_"}], "]"}], ":=", " ", 
  RowBox[{"ArcTan", "[", 
   FractionBox[
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       RowBox[{"kIn", "[", "e", "]"}]}], 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"kOut", "[", 
         RowBox[{"e", ",", "u0"}], "]"}], "r"}], "]"}], 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"kIn", "[", "e", "]"}], "r"}], "]"}]}], " ", "+", 
     RowBox[{
      RowBox[{"kOut", "[", 
       RowBox[{"e", ",", "u0"}], "]"}], 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"kOut", "[", 
         RowBox[{"e", ",", "u0"}], "]"}], "r"}], "]"}], 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"kIn", "[", "e", "]"}], "r"}], "]"}], " "}]}], 
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       RowBox[{"kOut", "[", 
        RowBox[{"e", ",", "u0"}], "]"}]}], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"kOut", "[", 
         RowBox[{"e", ",", "u0"}], "]"}], "r"}], "]"}], 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"kIn", "[", "e", "]"}], "r"}], "]"}]}], " ", "-", " ", 
     RowBox[{
      RowBox[{"kIn", "[", "e", "]"}], 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"kIn", "[", "e", "]"}], "r"}], "]"}], 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"kOut", "[", 
         RowBox[{"e", ",", "u0"}], "]"}], "r"}], "]"}]}]}]], "]"}]}]], "Input",\

 CellChangeTimes->{
  3.8735413196128597`*^9, {3.87354136431686*^9, 3.873541446430656*^9}, {
   3.8735415717184877`*^9, 3.873541576826765*^9}, {3.8735424454621572`*^9, 
   3.873542574461391*^9}, {3.873729533939823*^9, 3.8737297084349613`*^9}, {
   3.873729740827138*^9, 3.873729822160321*^9}, 
   3.873730350881184*^9},ExpressionUUID->"f9b93666-a5b0-48b1-94d6-\
6d1a85496d1b"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Delta]", "[", 
      RowBox[{"e", ",", "3", ",", 
       RowBox[{"-", "2"}]}], "]"}], ",", 
     RowBox[{"\[Delta]", "[", 
      RowBox[{"e", ",", "3", ",", 
       RowBox[{"-", "5"}]}], "]"}], ",", 
     RowBox[{"\[Delta]", "[", 
      RowBox[{"e", ",", "3", ",", 
       RowBox[{"-", "6.5"}]}], "]"}], ",", 
     RowBox[{"\[Delta]", "[", 
      RowBox[{"e", ",", "3", ",", 
       RowBox[{"-", "6.8"}]}], "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"e", ",", "0", ",", "5"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", "sps"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.873729834304262*^9, 3.873729962816473*^9}, {
  3.87372999361537*^9, 3.873730015311598*^9}, {3.8737300586286993`*^9, 
  3.8737300653434553`*^9}},
 CellLabel->"In[15]:=",ExpressionUUID->"aa0fae4a-3cd8-46d8-ac13-32fe6f27924f"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0, 0, 0], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwV13k4lF8bB3Bbxcw8Y5iQkiV+UiRZKpJzi8guSyJLRpYoRUhosyQp2dci
S0mL7HsyZI8wspR6lRZEsrSI8p75a67PNcvzPPe5z32+I8U4ZeHKxcHBsZ+T
g4P9evsCY/WOY8Fag7MyL+SOrEPPWZLGNpK2yC7c5v67g5FopaU+TFLyOFI/
Mtb3OTgWUeMrikhDvujRffpiumIiut1wQ2KV5Dk0Y3zo2nejFKQwcyxmxf0i
ijORPlHNl47eUjP/LQ+GorZfYrrxYbfQib3+J39JRKHGQ553n5y5g0r7R/MW
dK6jilqvnmqfbPTHy3hk1v0memPVJNTsnIN0u+L++zkYj17fKeBMkchDrCWP
xO8SaajFV44+ppqPAlrFlWaK0xFlo+ND5kQ+Wp/Q3zmtcwsZePq0q6XdR4yt
wP3VPQtFn395O2GiAAU+IlfPDeYgC+EPRsdsHyHD+M+r5+7kIsvTVhseTj1C
YoFMq9njeUhMQdN+LOQxatA5+31m6S4yG4iP6k8sRGvefNg8LVGA7sUZ/7W0
LEJDDU/9p8YLkA2/Iak5vgg9uJfa9LX4AbrqJGyv3FuETH1NnSZ1HqFvVI+i
bUbFKJmvJvmL+xM0yRyWLlEtQYK32rhmB0tQVJLJ2siPpWiJWlFwWrMUKZKl
awtEytDHy7nm3++Uoj2Cv39oGJahCvcLmTPHy9DOz6feVT8uQ3YqqhrflsrR
WCmzLsG7HGV3ZPl8lahGp6H04JXXFWghtXzr/65XI4mv8mcH+SqRnnvnWN9i
Nco5JJJYubsSTXH/sq5h1aCEnrNm25Mq0e69phpRkXVocreTq7ZBFep5ssy1
eeYZEkmotOzKqkYXOxQMdzg1oeWopqYNgnXofnGl7MOkJlQaHWbNvasO9aXu
4/7vRROim+mOPz5Sh2TcD9etU3+OzPctdprl1aEO7nBFDnoz+k66YvM/5adI
eO8bwZctLUhCL0j8w/56VPgk6s0JxQ5kOtcpvV6mAU0KTl43dO5Akn0bOnyV
G9DmAEMtucQO5ESKEayBBnRHk5w9ttiBqLrCekIODSi+/bqbXXMn4jKKU+RO
akABYzdn9e270A3VnyofOJkoi9wZJrC/B2UkcMeV9jCRO+pNm/JlIamwg0YK
Bk2II+dnzZtoFkoK03u6yqoJpfKIjXTksdDmmdnsV45NqK3dTbxggIXucT5f
Y+vXhOSslnNcNfrRyYSz675nNqFxj82P3nG+QuLBU3on55qQZ/z5+p7YAfTe
XCOGlPAceX/aMlZWOIwu3HTXjnnajEYdh3mZ7cMo6oaC1Y7WZmQ5dFWx6+Mw
agnI6X7R04x2d34592n9a7SO7Lp5eqwZcRffpQlHvkYRC0kZ/XwtKDVESuus
0xvEyXm8T9+qBTHpomnqtLdIo7UvjedjCyJJhW/+3+5RpOR3UbfueyuKoXVu
zT46hv6Yt5MqHrcjxhurZ9wvPqM31ZO648adyKc5nHFq8yTqEOvIjS98gYpP
UyzXGUyj/TY+11b1dqHsYK9JJZvviPBVWpn/0o1Cs5YPLkzNIvLBbdNBL18i
DtYo1AXNI6liqkekSQ+S9XMLtNr6A9nrm1AiinrQxx8ph8hPfyJTrdHHHty9
yPCT5/l4n9/ons7GDjGrXmT6MVZeXPwPsn8/+FskoxdNNjEutT1eQvpZF6yK
B3qRosDARh2nv0hFKFtomNaHSOiH88j8PySZUd2Uod2H9Kftfm2y5oBcbhXH
5dN9aKjxGKedAicU+1yoDEzvQ+VXnztK/uCEDb8UZgSe9iG3Vf6j3d1ckFJy
dFvwuz7k4htcTknnhvtCmjNSf/rQdnf5bmNnHuBnVEdtE2Shl4IKZH/1VfCu
W3tRZAsL1RaGZihwr4YMvx/k6j0sFMy51fzF29XQu2gksMmEhUgjQQe8y9bA
SuzG+dtHWGjn2QLiZuUaiHwfZvMS21lwqO9JzRo4fSEzncOehaoOqDnMNqyB
qEShpwxst/IZH7+uNaAKhN9/DizEjDmWEfR5DXxIkivJc2ShQG3Tb1fW8QL3
eFrhNWcW+nRPKjErhBfa27LSf7qx0FT4+0XuS7zAwbA/Je3OQnOMbCePMF5g
Pr8uYY7NISElv+MaL9hu9RMswF6fItnYmMILlj7ROrYeLGQeKfH9Uwkv+Ixo
jZQcZ6E6943GChO8oGHKo2x0goUS5UR5qq35wLbsS5mgDwvZO8rmjx3mg02r
xaxVsWUSVQyp9nygwFk1aI1dtmIS68Lgg1ZR/fBU7P6BUDHqKT5INs8JFsP7
kB4xpepylQ/OzOlPip9hofj3z1yJWj54XZv5j/BnITuRLt7d9XzgJbplSR57
k8nrhwwmH6RbiLwxwC6pWpitbOWDHCWt/RHYfTe3XGD084HI8vWmP9gCWgnJ
ldN8sCrmvfdoAAvFpru1OkuS4Pr5c+XJgSx0NiD0TIg0CVrN1t4vwna0yJRI
kSWBhaDopQ5sBdLA2RcKJKgzCRn9i91+br/cLnUSBLufcHE5x0JctjJRFAsS
OClzBW4JYiF/kQ+GlWHYuXyUnGBcn/m/P3uvkKBUWLO6Clv3pWjuVBQJPM+V
Gr/EFow8uCQVS4KSkWjZZezCnw0Pr98iwUmLzHtWISw0/iqLwignQVK/2DuO
87geiQ7dlC8kUBwuOKx3gYXESkkbvCZJoK2iFnYE+3+9le7t0yR4pHgy9jT2
MX5BzisLJJj6xW2Zju19rUVlhZMMHM9Kzk1hh17anjYrRobFo3usoi+yUMEJ
zmMDFmRQOpgRUHqJhU5EFxapHiLDLd+g1c3Yig+O/I23JQPZlzdoALvkS3my
2VEy2GSaUxax6xie7W0nybBbyvud5mUW6jnMUqy9SoYLw1XoGfZv3buLWfVk
GFN/q/AklIVCarPCnjPJsCpa6GAdNodyOjHxnAzbQi84tGOvkbgppdxJBubf
P4pj2MKLgQbPB8lQ+9egXCiMhVQeG6eNfydD/81UqwDsCml9aWKBDLsiq/9d
xtZI13684xcZnlaciLuBvS9yJzPoLxmc32efzsM2d5acoJAo0DI/I9yL7S20
sHuHNAV8gm15ZcPx/on+1mgtS4Hfut5aStj+XBPGQVsoYPTPz0EDO2Tm7dGm
7RRI/nbqkCl2dHvrVWtNCvyQZqr6Y98/nzF0zpoCGd9/5NdjKywkMTIPU8DD
Y+pUK3aRZ+xU4xEKDFUpS/VgV9pEcFAYFAiRaEHvsVt2nJLL9KaA+f4HJ7ki
WGjs077AxkgKzGr5xQH2vqzunK3RFJDlaY45gH3nsF1XfAwFprmGfcyxHTt9
NrkkUWBnpfuEE/brojsvuHMp8OnpZOZ5bHVPhV9e9ygwukmyMxw7VbpKqr+A
AgdSLd5GYx9KfhmQV0QBtTOdpWnYfcH/JPfXU8Dt9qstZdg71K4ZPWZSwNWw
7koNduw3oQChZgqkHhhub8A2dd7W+amTAjmP+ri6sDv07P2vDOP61EfBR+xG
gZr21gX8vO93zhJX8LncuX9h+28K3AlilQliXwrvFU9dosAxeREnEWytX+Nn
PLgIOFH+xF0Su2ZERJyPRsCTDeMGO7BFk3MP+NAJaFa6pKOGHWi2/cywMAEz
n39Kq2PvbNRrK9hIQFpPRhpgl+T7+xrKE7BkFX7YFFvAmfN2iSIBl/f2RR3E
Pr3+Rut6ZQLGmGo5VtiKN/LEvu4m4NPs4at22Df0lPQtNQmo0JY95IA9tVLr
U4sI+OMvQz6K/cCX1RKtR4BIXfdWV2w+BafZeQMCqo5eTXLH9vg0ucHehID9
/DFfjmPLHubyUbDCv5eSqu2NHSEQk5FoQwDn/AOD09gfO0Rblu0I4Doio+6L
nbN3x4YuZwJ+7XHv9sduNCWFLR0jwFmDZnsW+4PT2OQWDwLWNTh0B2Jz+dRZ
HPYiIEjbWjEIWzo0qeaKNwH0HcuBwdi6Cd6byn0IGFmxKQzBPpanf23Mj4AD
G727z2OHl0vOCQQSwP1LZ/gCdl7Loi0EE6C5+K77IvbzwT6m9wUCQmJ0n1xi
3+/4wy23LxNwVSIo6DI2z5/wuM5wAlj0CKVQbBmy4+JiJAHfejx72dYV2+Us
F00AUavkEMa+/jZa+6EYAkJVR1hsh2tNKEXEEbA301ctHPuuWWNqaSIBV2wW
L7PdcjSD40MKATz1/tVsf/bx86BlENC7bfoN26vDTHq0Mgk4yuk6wbZsouzu
k9kExMR9eM+23t2VrIw8AmrVvJrZdqsYWtORT0CqDk8i21dai0/9fkCADn+p
Cdv3hq4NyhYSYPg1ZI59fy0TLsi6mIBHcu6hbH/+o5kfVkbAqPDZv+znW00R
5i+pJOA2d4kL25s3zgSM1hCwbc/GcnZ9Dii2vaPWE6DLWz/Lrp8HytbbyySA
oyxNlO2r5kGFXs/xekeUK7Drfd/ZUji9Fa/fc2F59nq0+SpcaOsgYKikcS17
vcbDVn3+2UWAUFL9BHt95e5Vllv2EzC4q9b6HPt6lbEbQwcJMAhomWL3i0fb
8Yii1wTsydt5KgC7YHKDNfGegLPVdjvPYHcsLdTt+Yj7S1A02Ad7ktIt4/kF
94PWwYensLduv7TQMk1A9v90ur3Y/X/mY8KlPwS8qmgSd2F/P/zpUuFfAhjr
Pzxj74/JpGSXtxxUCOnKNXbElq86oKKxhgrMpv+22GI/Wn7UN0+nQsWqsm8m
2C+IKxqbRKhwZtX5aUP2fhN3yjFfTwWn04JD+tgK2gK+jyWpUPq74rQ29uMI
fwH3bVSoTvEvUmF/XtMhNEKJCgO+mZJK7M8v6M7nqlChPjrsvAL7fhlrX42q
U6HxRtJfGXZ9UWmqnR4V7GUcRYTY9fyZznfOgAr3Db5N09jrVxgalGJMhfSK
o8UUdj+IWdj3W1DBxWLlHzd7fyx+lzB1ooJVaX75LJ53H4uGYk8wqNAE3fNT
2NIeDZzRrlT42rYsPo6dM3BzrNWLCpxR4ZrvsLNKFfO1A6lw9ykzsR07/eQJ
RbU4KqjvW11/iz3PZayyLBOpML1zpDsZW3RkD803hQpj0Txdsex5bkieK7xN
hV9Mxzj2fE/a/KBc7gEVyDw1Gl7sef3+i6ZYExUkWilr1bAjrV2MuH/g5335
maMWnz85ZiP2rr+pMHQ3SKAEu87A2rt1Cdcv6BmlgH0+7tWPi+bih7GXvk3J
2Pay8oN0Gj8YJogr+2Ar/5pjyMjzg5t2oOgm7HepoUH7nflhoPvALj983i7G
L0XnH+MHKdXh+x7Ya2/43ebz4Ietr16scsA2uOzW0OXND/ez5y/vxy47brj6
UAg/vJscDGOf/9c0BOPdUvjBjL+mpBjnA7WR7ILILn7geR/x4TXOGzRFj99P
eviBEi5f14U9dVFRf4iF31/wCW/AzpOu/Sj3mh8mrOZG72Kv9eqX6PjMD8mx
uztOYf8J4TZazUED7qp4Zw7sYfJO5gFRGhzd969hLc47Lapp4yPraTA+NxOw
CrvMYZnfR4wGWTx8Yj9xnop90uiYJkED0x6myiD2AUuz5Yn/aBDkyUxLxa5K
99gVrUyD78s3ddZhp27JeNxlRANJyuRWMs5nC/N6d8NMaKCpQ7H9jfOfWf3c
LXUzGmxYt3LmE/ZqS8Prdy1oUDzhefoZtn/Ioud5WxqsLI+H+WJbvLSR2+ZO
g9WXx0Vf4TxJ8V+bez2UBrs+zctcx3nzYsONVKNqGig+ON7RiPPt5NuGxMUa
GvA8dWE+xD60NB+bX0eDtjEiP5Gdh3faRXE30KBFP1zLDXvooWxQbQsNAjYE
ZPJiK6c8OyLfT4Oqw2ayhmfx/Ds5K06eocFno2y/BpzHTTZY3+2QEQC50eW/
QTj/a3u3HhGIEoDkKDeHMi/c72IMmRc/BEBhm92Cpws+H5ryk02PCMKoYtfU
sB2u57bGrX31gnCO5pBUdRDP53me3vatdOBsthQn67MQo3hOMuomHeRd8166
auL8Gmq5bBhHh7IKoUFVbH/LskFKAh1u2Sj2cWNf/eEfE5tMh/FgRm42/j9X
qL64lHKbDk6dSiUjGrg/mSsD9x7S4b+MneHm6vj6fZQbz1vpeD330bbvxPsx
96THlXY6PJDtj19WY6Fcv26dA5102N7ly9eBXSEc+6ezmw4q8R/7jmGP2NE9
WK/oIBA2O5CmiufHmKjOhzE6TF7yvMepwkJKZUHieZ/o4L9G36tbGefXiDeL
rl/o0HeLIZOBbbz5dtHEJB2W9VzOqGL7ekmKz87SQa1oNNxtB86jmpcXS+bp
oNNBnlLGjiA+9Pv9oMPPsBj9FSXcX09yo3//poMtf/6bFOycSzzuNX/ocG1p
39pj2A8Puu4LWaYDo26fthJ22aaWjVr/6KDkUMBY3s5C9fOyiysrdDg8fi6w
Dfv/5P+Gdw==
       "]]},
     Annotation[#, "Charting`Private`Tag$11606#1"]& ], 
    TagBox[
     {RGBColor[1, 0, 0], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwV1mc4lm0YB3AjGmh4WiQrvSRpkKb+t0QZLVmpVDLKSkplNFWiSEJKIiGU
qFAkEdmK7IyMZ/MMUUjlvfp0H78P13Uc531c5/88VRyOWziJiYiIrBEVEfn3
jTvvILnS0X+T7dM6xSfMPShtUDa3Ud6LLeK10nK9HzBRVhiorHwM83zGLyw4
XYXpEblZ01q98cD2iJLRjDrEFYUqSSj7YvnYYFX14S/QEjiGTbhcwOJYb96S
Q43onP7w7++Wy4hdEWUqqtoMtZSper/Gr4Lbu87zak4L3PV9PEaUguEw+UXg
LLU2vGrsTho2vAm6wkyu3qmv+OVm3jHocguh32csqkpvR4JoR/CI0W3Iro85
41vdgS21txf/bIlAuGKdh3pbJ9gxW4uHj0ViqeP9xWrNXbjp+Gf/0HgUWF+t
9u4p/YaG8aORQqV72Fglc1v9XTdOlyuuELy4j7JXsSbGqj2Qv9NYzTN8gPXJ
6p/9LvfAQZMS73eJh+v+9pSS1b0QG/f7PByQgI+VzMSg4F40dlstHDZ6hJOG
9C+lbb04+0wq73tLIq7Ez/mb4tUH0wim5PeEx1jpuG4w73UfFM4WWw4eS8LB
FToVi373ocjwjFAwnowz5zvcEs/RMbm9V52nlIZDE66z57kw0Fr0zmeAnYaT
1Udiah4wkJ4SU9L/Ih3HR+3qmusY2OG94yDX8Bn6vyi/1NdlInpqfjTLJRMx
bSqZllwmss+Yep+cn4Udc8dfv53NwgPuGrmhgCzYP1k7eUKfBc86WZfvRi8g
WXbRwCyUBdkHFWKDLS/xRDFx0hJlNsan56Z5bXyFeXXGk55sZoN+6fEuYcIr
1EmdcbJwZCPX5fxDwbFsNNeWSe5KZsNOR3c9fzwHIva0jFRFDgxTVHo8DuXi
i0VdhfUGDrTkZlznlebi5G7WXWcbDv7+5TQNhL6GZIt/1/tQDh5VxZ/oV8rD
g9Jhm44hDoZjcjS/3cyDZtVIjcp0Loxdqvu+jOWh1PjwyRZ1LgbER6zyG/LR
XXTuZvJeLtbq71gfHFSA2I+SCQavubgh5TgUMFyA5RsP0c/XctHV5vvM6/A7
VC+1n2nex8WV08mKthsKsStdS/SxTD/qMn+LqQve43L6+wvRB/ohi6n2MtpF
8Je1W//Eqx8b/otcH+RYhEA1tUj7wH6EDacN+dUXwXjJeFjVk37o3m5ycnha
jHyt9q7GgX5cqNIyXXmwBPfWJl05dGwAqS9e//c0qgQr1J78rj0zgC8xm8UX
15Qg3+5a2MTVAai52BbMX1cK2Wl5btmPBlAlfkVbhPYRHWLBL683D2Cufrvs
57IyBHOl9b7q8UCpOQu2/SlD4OybtPmbeTgmNVj9Qacc2qKJMUu381DQJnk1
N6EcMp17RpsdeHA4vWokzq8Cl0PHj1jd5OF5ZnC7u3YVMkXbCiNaeeDKcm+a
Hq4CK8zZld/Ng/pp000akVWg69v0aXN4SNgo9ahvrApNP7ukD43xEFF509nu
YzUs9HPtbs3n43TfrcGt+2uhNeeplL8FHy+NhYmLb9XCvqasx3IvH4K0XZbi
H2pRIbshV+4QH0e9Zr0uVP+EzPqABCcPPuz+RASsHvqEisa102YE8REvVR04
y6gOl/K/5LBy+dik77//on0d/rY8PxtZwEeH59LVgjN1+P5rW/KKD3zM/3KD
WZNeh//EP2isreXjdoy5yfWZ9WjuuTw3pJePwMWfpk+01+PUrymLjk4TQNXm
PMvjRz0iWynVxTMEKLquXdQx/Qs+V0ZXNdIE+N0fdiLf4As8bKzLpRcK4PNy
Z6PPky+w0a9YQNcWwAX19wa8GxDqME93zy4BRBJ/5rffaMDG0AiF9ZYCxExS
6KhKasCsfYU35toKUFHprJjW3ADtr8rnsg4KoGH5O9FpfSOm9chMdfcUgH1U
/VmXaBPWJjUtXhAiwMVq89pa+SaoBbgfPBoqwHxtb36BThMsFOcyn4ULsG2o
YEWsUxO21mfpK9wVIO2cRY5NZRO2Scq/jn4sgGvEucK68GaIVmkUTH8rgNiP
xG/v05px7+rz9PZ3Aty3qRDJ/NCMIkl9s0dFAlQp0AxDh5uR1TkknFsmwNIn
qeUmti1gTCgyMuoEGHjbUPdBqRXP+/XM79MF8GQs6ct+3gb/LTv5AnEhuu3b
phRXtsHYplhsioQQe1qva9fS26DQbfVLQVKItdUsX4b8V2yJDH+lP0UI8RfJ
M+cGfcWNBW1a+6SFiAlQ2XTmYDviH1q2Ss8WYtrPuiOBfu24NHDo7U/ic8cv
Bt+Kasd1lb2jnXOEcHD41vikuh06JWfvJ80TQmtbnFvrmg4svPdGh7ZAiGKa
3L11MzvhNT9Hzl9FCN3QikKjpZ0w8r8crK8qRIrkWfpu407YplS9+kN8c6x5
uWtAJ5oCXLafURPC+ltk2X12J1RundbYpS4EN33m8HhRFw7rKbEeaAmxX61I
bkpHF9Isi5z0lwnxOe44Zo90Ya7XSn47cXb4pxCtZd/gvMpem7ZciAunb6oc
uPcNI04iTNeVpD6VK+rf1nbD0tz964fVQlRsvL38vFE3zludNtDXE+Ka7cM1
Cy26Ib5ZVyuXWCz8zdZ9bt3onMyJebxGiF9/B1xa4rpxVn6bqsc6Ifo7rFI/
i/UgfM2lU+82CpE+6pDlOYPMtU4rbxV9IY7O9nojo9CDOzSO8hViullIhdnq
HqQeuXvfaJMQHfmF7HKXHmj5vDfOgxA1MepLimp64L2Jeuy0WYiQbN2VB9t6
8PPO98Yc4m11Buv+Mnpwzaeyc5KhEKWT95tsnCDnt2efiycuOB1+7M3KXihJ
iuqUbxHCLyLuhM2mXtgKLtpNNyLv4Xm670/TXmD2RldL4mxGabCuYy/2iNwy
7STOsBxLy4ruhbT7R+NeYyHcvSRf7nzcixXrDrGVtwqheZOWz8/sRYbOzMv2
xCkly6q0Knsxo3dLcjPxw1UO3NTxXhwbUfyWv02IsJnVmo8O9WFPXvQrS1Mh
8tYszNnh2Yd9SeFRQcR0++P47d8HH3elI3nEGzJoljZ3+6A3bFEpZyYE22T/
OZlPfRh1sBuvJZ594vnkt+19eLH/6MlRYsSIRBzl9GHGtrUsFXMhoplJKSWT
6GjKbq08QWwYyPvsu4GOsy9WZ0lsF8IzHXbqJnTM2x22QYP4fv1teqM1HSzP
1HoTYqGS3thybzpeHty+OIRYYev1wM4LdBiUHmGmEm/1/CpzI5QOXq9FXhnx
w4LzqsxUOvxnJ8dOEFf21T27k0vHrxHv5/I7hBietmiNQSkds3YI2nSJzfaW
mz34RodmXNwlZ+LRHzN8LOYyMPAq0eA98aKFDhMTixjI/1C/tYF4x5bs4IyV
DBTkPDrJJE6OsH04ZTsDLzg2m6ftFKIuL10j146Bpct3suWJf3f/fnnkKAPv
fKbmaxJbLH9UXhjIQKl2CXMr8Xnrod3utxloNH2zxYo47ZxRh1w8A/M+7Ws8
TNyYdNe5/BkDbebKdz2IRWo4wlP5DAgVfcLOEi8d2uCvWsFA1tGyt5eJreXD
JOqaGOhadE75JvFlg+5b5/oY2K0x8S6SOOPoKvmlg6Qe3cyoB8Stt64ktf5l
QHVkIPUx8aTXzdrXpJlYO2fOWBrx8i6NPB15Jr7qXgrKJLaT8DfsUWeCNum6
VTbxNa3a2rDVTOw2OOf4hvjFHiXbjYZMuJa8yHlL3OF3opeziwnPSydNC4mn
JJa437Vn4sx9GbUiYp3KOSNb3JlwOlBhWExsL3S59N2XCb9dnPR/DpmXL5UQ
RPa4ta/s/jlnk3T09igmTqg4W/473+1krzyeyETIEdz7d790aFZ6ahYTcc98
NQuI12SLrbYuZMJEcvvUPGKHdsv34jVMSHTy9XKIp8X2eJxoY2LR7xXVWcQv
7TwVvjGZuFcszH36r1758SrzYSYml5uNJxOLfw3yzRdl4d3uA/fiiZ/dm62h
MYOF27abImKILfc+ao5SYGEgWIIVTpzUmq9zYg0L7K6o1AvE5jFbe7u2sFC2
5NK808TDNo3h5hYszB2LZroRG7XweOoeLOi4eab9ew/8aL8HUX4saLdQ8SbE
0daTzcSvs1DjeHBsIzGzSTm1K5EFt0M3ulSIw6IyrM2zWEhgxbjRiPWs1kvk
v2PBvfb7oUnEQY17HKJaWTDJ9blOJ+9Xo+Gagvl0NlKGMxlxxHURtOq8BWyI
3/6vLoT4rEWCr/oSNp5L3NA8Q1xZn9cstoWNlRuqVu0gdq0bCM/zZSN2v0/6
D9JvsuG+UA9io8L27bIu4vydkvzISDYWDMlp/OtPqc9KZl6ZbKTKRV2NJH5W
ayGhzmBj0/Xd45rElqHfXkV+Z0OGr7dzBvFvc3cHMREOlBsGlYdIfmyvufq+
U56DxBfqBnnE/Ko3vpG7ONjuPRaCf3kTYqQhZs9BwvvxPiXiTaZfmo+7cWB6
W7N4guTTrcp+HbNrHJxIVpAqJF5RocgXLeDAOvnYd13iEx+vOBwne3VoSPFL
aZJ/inRjyZ86XATppJmyTYSoFpuaHkBxMdOj+2IJ8WIqdDCY7N1WSd5VZ4nb
8qMuJt3ggnfcsaeL5K1BVsrDNgG577Vz1D2Sx7NiK74avulHU+Ok+BoyDwrz
Qs5Xlfbj7qUQ0Thit1Zz1d31/ch9t1fWg/jjnPpj9tx+DH6Po6SJfcPbRs8q
DCC2gNZhTOZLz1Xu3OcXB3ApZebxDAPSf17Se+ab8NCh4kffTOZZzX2ab4Y1
D625IzISxMxS+fjNjjzsZ2Uwy8j8k5db0u9+gQfB6/gnJsRXio0Ci3N4mND3
PGVG5qe17IVXrqp8TBn+GketJ/P21SDt3Tgf+R65msNkXj/yd+OusBHiRoB3
UxTZB1YUTR0NOiDEaosbm82JiyalSnw7IkTC2yvnxYi7Q+nKoV5ClHvJO3ou
Jf834YANJ1gIFu3jgKEmyf/SnaWPCkjdd4M1mWTfiJDRfSirOogLO06PTVlE
8iz+9+7hgUH0J8lyG+aTPGzopgr8hmCifOviFrJ//XfK+ayl5g9Iq/TIrGwU
gP7jrrXUu58wc9MQq04UwJThei7ixCjauktm7ncTYAc9fKmi4i/09U7PjNER
gFvicLEiYxz6Ti5dH//yoT2reaHhwT8oTzfdyirlYxp+HO4Y+osVL2UiIsL4
2MqzG1G1EqFG+kd3he7ho/WDo6idlii1fsE62QdyfORcL7VX/iFKBSoZhaf3
8eAs4dP96ZMYtdb3cpdjGg9HvP1zpO+LU14JN00MTvKw3GXpJ/PDk6gNsZVt
fpt4+CyrJeWzToIal2lyjZXg4e3zy7Fa4pJUoJn19oG6AfiLau6q6ZSkRP/g
lXHMAKZ1+G3zzJ5M+VrLFKsdHgAjRSUyPmAKxTi4nb9EewCRGnKT8qymUhXK
e9sOjvQj/L5z+WHladS+lsauV+Qd2kUe+CTNmkZ9KNVb5XenH6NbksfiC6Uo
34jWO7ft+9HH2Hz2Q5A09Xi6gdLzZf1I1F+5oPawDOXbIh+tM8ZFxlWfWS7L
plMP+AW2mRWkb6yOmIn/mE4xtBZ3R9zjYnXHo7Sg2hlU6D7T3b+cuYhZEptR
azaTcooN2rh+LRcXikJjzPJmUmYHXPVoklxsX2CVXKU2i7rjEdVm28SBgWf5
vlnBs6hU9a5mm1QO4hUc1Gp+zKKsVz1eOeUsB2olT6J37JOl6hdpXjhsxkH4
sg+aXwplqRkdElmSyhw4D02qr9SkUee3Tkk8OciGw4vvysG3aNQnMacGtzI2
3C/v+W16m0axpGIePvvIhs+e7BbpOzSKkz90glfKxvUfPmHh0TTq85pFK4+X
kFxcNzZ+N45GFXCsZp0sYmOseKI55SmNqpqRZBuQz0b4F+nQ0nIa9dro5UTk
czbuP/Y4eq2SRh1gzvvcnMHG41OfDLdV06iTUXdT5hPnzg3/Vf2JRr0z6/N4
8JSNDjva0YYmGhW4J/5AUiob6n1yhr19NMo6ybA4J5GNFdl+ikkMGpV8nCY7
+oiN9Vfbx5xYNErxzw/X9cTm6nFZHC6NOndjUPd9PBvebsqKg4M0qjcrZ0nl
AzYCNl4aezlEo3ZWTH8rRXxVprfx1A8a1T5+2nJHLBsxmY9vjI7SqH107/iG
e2wkXpzkkv+LRrVkSuyeS/x0t9PmgN806lRp8rS9MWxkq5Yt3PSXRl3IsKiJ
vctG4dB/YxMTNMrgp1R0VzQb/wOdwtfF
       "]]},
     Annotation[#, "Charting`Private`Tag$11606#2"]& ], 
    TagBox[
     {RGBColor[0, 0, 1], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwVUHc81f0ftSIj6/ajYfVkhaR6SlG9v6EIqYyECjcjiRbJSIOMn3AJScqI
IqPrScN95GNUZO+RkXVd8w73Xr7XyK/fH+d1/jnjdc424jVrdwE+Pr4ff/B/
fhFKFN7tFnzkbZe2BY18Fn3tULG0V3GA49UFYYqJ39Da94owFRUvCDBMknYz
a0aSiR/JYr03IVa9cZ+vVzt6URmrvE4lEGacrQ6EKXQhHYZb3JrnPXDZ1rhp
5XoPGpR8+Xul5yFsWgoLjn7Qh1Rfi+5fWn4EIhajt3Uu9qOrh/19FpWj4cve
Yx1XZAbR+87hHI7xY7jq3iCakz2ElrwtB1ie8XCE1TeVKT6MTJoS1BZ6EiGa
Qq1RvjCCJlNNqzheSRCXuVi7f3EEPXZbPc9eTgbruZ5VwehR1LF8OYmp/Awc
P9llbn06hm7XKukxStLgn/BbJVUbx9GWJ50Nc8bpQDts5zL6eBwRtTDBGc8M
sGnyl3bxpSKB5aAWTkgmNHLc99zvpaLOYTtFzrEscPNpqco9MoHuFIqXzfdk
g/Fw19+b+WjIPHFCeD7zFZwYjD8w4kRDCneqbFleOfA0VWQwp5SGKo0DmIzl
XFjIyGjouTCJRPpHNeaU86Gz2cBF9ckU6q384j87mQ8zOy9PqPZPobevU2tm
St5CS9RGR8Nt08jqppXztHEhiN6/LbM+bxqliFJSaJ7v4IJfwBun7BlUGmB+
89YmMlywunTcdnAGpU/rb2aHkGFJMP7oktws8m2V9Zw/VgKa8jbL4RGzSDa9
ToDV8w9YxHsFT9nPoWXJj/nXD72HcINDD95Ez6HxB69OMzPfQ9rBMwK+lDn0
0TP0JcOrFII1WpT1N9OR496/DejLH+DTIeJwUiMdZdVn3JhRLoM57lrf1vVM
xEn9oPXrcRn4kvxeBygw0XHPhrF2XhkEHD/RHqLHRLOCi3aUDgo4NduXxtoz
0YHDVgbRkeWg+OVNuEQ2E7W+WxHQYCDYPKU2Y6bFQrIgenGDbiWEBuS8xQxY
yFA9ySDSrRIi6kqU+c1ZKI6Tzw5qqwTLCng1cpmF/k7ocicWVEFo17V/TXJY
6F69jvlu5xoISIm2jZOdR3kln9QLkmtAZbhGJUN5HrWnGgmqNdZAdqb3aLTO
PFL1PFe+6eBXuLRzLJB5bB7VC4br8hG+waOGJwT5gHkkd7hftuX7d/igJmfh
2jGPit9F91/VrYdnI2TK1wA2mpadfmzuWg8iU6+ycu+zkcZt8yOaSfUgpVNo
7xfNRpmHxLPGePXQUpxhMZDGRok/Hns4fmuA929EygbL2ej2WDzL9HwT7LT9
emZwlY0yxBvCZI61gpHf8j5Tfw46cjj4/P2LrZBzzf6LdQgHDfhq72MEtILs
yEimdRgHbWqPmWh82wo3fzrI6CVwUEKq5Yko6Tb4KTwW7l7IQWFqzZJr/W0Q
Mmuk2DvMQZ7Q9mz2ZgcUbGk7WWDERXzZC5T+mA74xTW6RjTjolQhhYH6nA4o
Y6GfMlZcVPfDQym/+4/eeWbEwYGLNG1Xst0NOiE6/vYQ6RoXTV7WKBzi7wKp
JrZ4QRoXXUm8W9FK6oZJvfQy4VkuEuBm/0L53ZBEl6+eYnJRmn0d37vqbmjj
u6tRy+WiegWCcSynGw7W6R2+scZF2m/yak+c64ELt8TwMNkFNPtvR2u1ci84
ikYbF+svIF/qjrHS4j7oG3RzuhS6gIYv9q2v+tEHvLF9fUJhC8imN0q3abwP
Jq8ePJkVuYAONNACqVt+gjb2wKSBtIAES3Kl5SJ/AjVrLnM8ewGlhmw7EuDc
D6sNx3bs/76Aqgibnx2UHoTZPGlWiugiEtsWrvHrwDBMPd9dIPVoEdUdStgV
emwYdCcFPG5HLaKIcy/1Fa2HwXzgif/PmEUkQPps6uQ9DCTnxzlpiYto6fes
Z8+LYcjcTqYsZyyimQG7vBaBEShxsfYxoSyixlSNHZWNIzCYoq8tPreI4qQb
tLJcxmBjHc14pyWOyvQVP1j5joGHQ7+a8SkcjV+8BivBY/Bs9dvoWWscGRYR
bO2fjoExv/VSwDkcTZ44f3dD8xjMBwXgmW44Mg6bawk0HIfqbcQBUgiOcK6U
v7UcFdLv51qvvMXRdkXi2tp2Kmi6pZ3sKMKRlUlpdNFuKhBPL93NI+MoN/Hc
y/UnqUDWbok6+RFH1ruyaivC/viNwv77sApHRZf3bNFmUUGUK2uZ1oOj3vjw
nN7ff/IkSBS7nzgS+tStGyExAfZDGWFSgzhyXBdsPKIxAT7s7bqhozhan11z
9enFCThLjFIxm8MRsd8WCTZOwPSRcguSAA+JPR/xudE3AVrUu/Lq63joH0df
hV8TE2DjTb5KEeEhwZ+RgRR+Ggyfe0cZkuChwmcbNTWlaHDqt9wnXykesnXI
6k5WoMGslLrJqgwP5fRS9t7Qp4FguaW8nDwPWaaajg6Z0KDYneqSuZmHOPad
JEtrGljlaxjtUOChdHlXoDjTQE6psZysxEPHeubmNHxo0G8Y1b5/Gw/RU4LS
k4NokLIYuN1QhQcpZ0UsBKNo8OROUPsnRR4ckUviXU+mwS4nu5F9W3kw0aWS
N5RNA5KTq9WHTTyISy46a0n+k6eHy++X48F+O4N1lC80OBpx+9BnAg8iO22I
yb00aMLfxldK8kCzI0LBUnISarsU9owI8eBK6yypLHAStGOei7qycJAlBYJG
5CR8jc7+ukDHgXJKmJ6UNAnlur1VsbM4iLcoW1x/Nwn6zashiIZDYZP1Og3q
JJi4CDTvHsKBXv85MOn0FITZumSdbsDhxrdw4jWNafiPm2mZz2sclMaPCy/s
nYbMyTuyp3NwaBAQfRuCTcPDYRfnvdk4qGGxrGiHaaiICG1deYFDHyX5fk7M
NGjf9eKmJeNwlPz6ZR9jGtLznYIEH+Eg87zup/HnGZjPJB+854ZDyXUJm00n
5uB4k3UXWw2HrGDvaT17JiwIRtoY5S3Cw4yVM5xZFhBjXwfqHl4Evo5hrDyI
DcwXWTNbBhZA3c/jjq0WFwgT52p0ghZgnPv0rPiXBeD5pcm8V10Ac+qVu4k3
cOh0+WRIauKC1ThJW0lpCcZ99lwfDOHCdA3xfl3RMgQG3Ch138cFXZluRWPn
VZjh6L+MmeGAGHBdB9i/oVe8eSM5nwOmc46Lf9nxYSvHhcXt3DnQW+3G76jD
jx0dVxpcVOXAh6ivF1W4/Jhj1Plp9Rk2eKzzH25uFsD4/2MqTy5mw6WbwR8k
0gQxoaGuavU7bNjlqd1s6SqESU1G/HXOhA0tsjri/gfXYVHrLTMGNrDh3+KH
z3UEhTFOh5PLo4F5CObXOt04KIxdXQkvSMibB7GBIDPfUhGsfrtNBH/APFBf
b0vKCFmP3dqT6798Yh6SNDcLldmJYhem/xJ22zIPpDSPWlcVMcxdUbhIksEC
x6QLzRI0Mexw2mpgzDcW4Ca5vIwKcazQ7MzjiGcsGKMa3amOlMCeE7faX7zB
guzDu7c2uW7Afie9E6o7zoKiR/4ynjslMa/8OP4BZRZE2l2yEORKYm20X8l+
S0zYN5CVH9kkhXl8eUkmtzMhdcfzoiYLaayYVrH9ZDET7lXGplqUSWMHNsj4
Z0Qy4eRWu9x6VRnsVGpfqCuRCUd9a51komWwJpkuYxOMCRkKRNVGrgy2ls+u
oW9lgmrNmxQrJ1ns1Xm22T4uA0g7q7XaK2Qxib1cb3opAzzYQm0/tAhYYJBV
iMN1BhBL5lWi4wlYy3mpcI/dDLj60GbFPIGAnSrfVHRKjwH+NqU9Ek8IWNqV
xKkDuxgQxfWPI6UQsHTZhnTxnQwoPshbfvqCgKnZ8M2WaDKAV7XW/bqAgMmp
sbUEVP70t0vEfq0lYDHPqXIlkgxIe+VzOeIHAfPiuOWlb2DAK79mY7MGAuaw
hltGSTDgoxxpqaGZgE1N2rQ4izFgwJFwuaOLgD0jRnyWEmaAxthm49ExAsYQ
8d5ya5UOeqVBSjlUAha3n9rvvEIHg0f9PHcaAaO6mJItl+lgqfGCPDVNwHhP
Pieo8+hw01tFicUiYHzi/Pr9HDqEHHrA+4dNwCzyfeLr2HR4tGG0049LwE7r
FC59mKdD6rtXMThOwIJGn4gmMOmQfV/Ik7L0559U/aJQBh0KzrgbhawQMMMD
2Zeu0ulQ+td3xSO//+xZ16bjOEeHCrY6b22NgHkbV603m6XD/wD6wp1J
       "]]},
     Annotation[#, "Charting`Private`Tag$11606#3"]& ], 
    TagBox[
     {RGBColor[0, 1, 0], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwV0mk4lH0bBnAhWSKMhDQmRFEq7cR1l5RsLUhpkT1pkRJCpOxkiUKWUApZ
k0rlT6lkKYzd9Bh71rln4x7E6/1wHefx+3KeX6719tdPOvHz8fHVLt3/M+2u
vdB2R199yumiNKbuCVRDpZhaU85AvZn1pmr3KrT4vfI+heIKv93qKyM2/kAS
8eXFop0eEHAsyOcl1oDSqqKVllN8YMjYRO8NNKHNDMeHiy4BEFirI7Ig14L+
SKQvzHcEQfBZ8WfyVVSkmiOye3YuGFpOxg4MY23oip7n1RmlcHDn/hPakt2O
3rTSn3MMoiAq76yQzUQHmnUzpTFdYmDmngjDmdKFDjXGbZjuiAfRV2/8gy17
0N+kI9Uc1wSoz+tGLadpKMrx3zn2XCKcmY9S/nH8D6LOXUrAlZLh9TJO4H1K
L7r9g7yNUZICVt78Nwxme5HCo9b6SYNUcLyt2tfuREf2GpjAuEsGRHrrb5fY
3Ie8X4t9YHVkAYwmONTF9yPj+GEh1rNsGAw2iagb6keK3tWWTNfnYNp+9bv7
7gFUZeCFM+ZewKjNrRO7mwfQip5+9UmlXCDtOFT6iTGIOqs+e078zYW5Bf6e
dTuGUF5O0tfxkjzwXJVwgHNrCJl7mNuOGbyGMDez5MOcIfRYpOLxiEsRTH4P
99McGEZlXsYeN+WKgSu2wW9s7QhKHdsjz/YrhnUFBSKxFiPoWpO0C8uwBA5H
eZu/RyNIOrWWn9lRClYlUfOBUX/RnER5rvv+NxBimJjh9/kvGryXfRx/9gbK
/erink/+ReUud9MZrmUw7rjiZILJKLLZsVNnau4tyKZWeKUvjKLMuowb40of
wLNMxb7fcBxxkt5q9EZ9ANaY1hP36+PosEv9QAvvA2QraJMdksbRhMCMVQW1
Ar6tvqIQ8ncc7dUz1wkP/QS7Nw5abg+aQE1F8/zqDARXqUNlpPRJJA0iF8S1
qgBTigl6jCaRrlqCTqhjFbjjf1NL6JPoISeXfae5CpaR80Uq1k+hnXFtTvb5
1XB/b1br8YwpFFC32Xi77Vf4k/xxjVEEA70qeaeWn/gVZkISgg1yGKgl6aDA
hoavkFJUS+ZWM5Cqy+lPcvtq4DcyCAgkGKhO4IEWH+kbKF/qrdLajiNZvR7p
39+/w92jNbhSIo4Ki8J7rmjVgbfy1vutRkw0Jj0WZWxXB8L8fuJJJ5lI/bax
/saEOlA5ILji6DkmerZfLHOAVwdONPzEretMFP8zytnmWz00m2oemklgotsD
Mcwj5xrh1BvBoH4aE2WI1d+XMmyCLP9CarkdC+nr+Z4LvNAEO4R7ah+5shDt
muYuhlcTPODqeTneYCG5lsjhhrwmOD8zE9cXyEJxSaZHwySb4YV8QJp5Bgvd
3/BLYrGnGXTW1KtIdbOQCzQnT3hQYVeXdbSQERvxZU1X9ERSYWOqefJtczZK
ElSk1T2ngkxB9fH/LNmo9qczObedCiW5vqsT7dhoo+V8lpNOKzQT0Vuy7rDR
30vqr/9b1gaqmluumOSz0eV4/8qm2HbwrqZ0W6zgIH5uVi/KbYfYrxaNtJUc
lGJdy1f0pR2M2K+ibaU5qE6RZBDNaYckM/nXJ9dxkObLVz+Onu4Avoenqbg2
B018pDZ9UeoEA58g+przHHRtaNNAWWEX6PbclrYv4CD6hS7h6p9dQPWp3PGi
lIMsOsO0Gge7ILxN913/Ow7aWz/iM6TQDZQkaxmzLxwkUPJCUja0G8KGS9zo
7RyU5Lde38u2B2yLvAzEFjmomiSfvE/yDyh5p7W+N+Ei0fUP1Hv30sH2DznZ
upOLavfHbb1rSIfFheGgPTQuCjmdvmfdSTpYSUQ4yNC5iD/2/ZGzbnSYln2j
+G2Ei2YXJlw60uigTbZn889w0TjN6tVv/j6IGszdxZGZRg1J6puqGvrA+XoT
ud58Gj2UrNfIvDgAt5rpkeSP08i+xxIJNAzD24CYCL+IGXTj2wP76+pjoM+9
oL1Jh0DkwcNC0zvG4FSMh/LD/QSq5xfJ88PGwGFt6TRTn0AbsGhm+Jkx2BNI
Pl5uQKCuisTA55FjUOO8c2i7GYEOFOekdzHGIFbYNZRrSyCpp7XdBu/HYZ1k
WtvbBwSq/BBxt65mHN4Xap/ghRDIrdNU+UTzOHzOsW3WDSfQt9XNrhfGxsH+
ydaBz9EE8ontIrwVJ4D+eatz+WMC9QWPyRYGTsBzdk+JzysClbivtJA7OgmO
H6Ltj9UTqCGF5FNwahIKA+x9shsJNFyjkHHQcRJGdxaGTf8mkIL8pvErAZNg
uDkzIbmVQA+qDe9Xv52ES6+c/Kl/CHRKOuDNZeUpyHB3cV2DE8h9f0jX4tYp
KOW+8z3LIlCEc/Rigt4UrL3+MTGdQyBU8dQEnZ6CTWv9BtfzCLTR4f0AKWYK
EhUEhdbx89DBaCSSmzoF/+Zlrp8R5KFz775v1c+bgjXbZ4cThHgoTqzN99K3
KcgW/SgkIsZDs2+YpM9zU0APUSMPk3hI5j9i30kRBjy6yVq9VpaHtIT5Lo7I
MoBY46JyTI6HjLRXhPipMmCw8cyRIAUecjgn8VpKmwHHfqUHlSny0N2Q1S05
wIBdO0Tah8g8lFSsSOiaMaCCd/mg7HoeKu1WITfbMGDDKXZvMYUHDYKah5wv
MaBXfzsEkXkwrKV9ec6TASGXhO9bKvKA78y+2Nj7DFi1Oz9PXYEHuwqO0CrS
GSD/UjmraTUP0q0cb61qY8Bw3qe1W8R58D7ALeV5PwPiNRpdhMR4QM31qNqH
MyDS3z6OLswDwQBMqHQ5DtfaxO48EeTBXqr5SRoFBxVJlahV8wRk+rqNbbPG
QU1nfOXGcQK2VYkQoedxqNk1E6Q9SkCV4KvlvQ44uJ2jFumNEECPHqREu+PQ
pFiseWqAAPKz89aj4ThkeoinxvQQkFJzrCbzEw5iDYr9Kg0ExIvvTJdWZkJf
5CfdXwUEBGXMn+BMMEHXWoIufpMAPiod+3SHDRWxv70jhAhQu+XsbanBhcyY
C5wbyTMwyH1ySuzzNGQ5/XczXGcGjIcu+8ffIOBW4funst3TYD4Yq0kmz4Ka
zMdNovemYeyrfWBtwRyU/rwY4qQ9DVpS7esMbP+B0cJNi9N9XBAFrh2NvQDd
FwMGIZkLRyZtZpSt+DBXlztmyeZc6PziuMxm8zLMzusC8hbmwtuwmgsU7jLM
2DG4OPUnB5yXe9J//eLHNjz0i18WxgEHD9+3K1MEsAi1F9tYphzY6qL5y9RO
ELvsn5WMy3Dgt/RmMc99y7F2s5v+qt1s+FgY9HSzgBC2bSLffO45G3yXaRxv
+COEaZI93eWuskGUdsfoWtkK7Mg1xQe0PWwYylmfkOEnjFGWMbbFr2BDwkZ5
wQ9WItjVVazbSa0siE1x/mFHEcX68qt/6uWwwCbh/K+VI6IY+wf98IAPC4hD
L3gZlWJYYKXRoVFjFgwMHfT+EroSe3zxSIcChQVZetvXNtqJY1D1lVXNYkJB
sKeUyxYJTO+v7vmzP5kQauVgIsCVwCw2Ss49zWTCLlpmbmjjKqwlcTHT0YsJ
SZueFjSaSGKvH4TJNRxnQkBVdJLJB0nM+syrY7/VmWC21upFnaoUZsi22mHH
x4QD136clQqXwnxD/W8b03DIULRXbeBKYbuV1jj1v8VB9evLx+ZnpbFDrrtj
CuNxiN3yRaOlUhqzP1pc+WXp75zZgs0/NUjY9d02Y09McbAvYVHCY0hY8IE7
zm6aOFwJspg3jiNhVLt7gZ0aOHhalHWsfETCcry1XxguOYzr+TD2MQlTmYyU
oGzCoXAfb+5JGgmbqZq82KaGA696sT0nn4R1/pvnP6i8tN+yMrrmBwkb/NK7
TUUeh5Tsq5dCfpIwyzleapwcDtm3fhkY1ZMwkzURJL4ll8vGztb/ImFdj5ZT
/sjiQLMhXaK2kbDayf3lT2RwUB+QN+gfIGFpHKN2CUkctpXdIT8fImHZJ/hH
/VfhoBPcw3MaIWG9GnJCkxI4mKqnFY+OkbBCIQvbenEcPNwoZCZzqf9yV3SY
GA5+++/xStkkTKnQmD0jikOweH/rLS4Ju9cn5uCy5KSi7EiCIGF63V52h0Vw
yAoUdKmYJWFu95Twt8I45J9wOug3T8KcaZLhG5Zcpvx9nf4CCVM8b7w5cQUO
lWw13uIiCaN51nUJLvl/vLiFRw==
       "]]},
     Annotation[#, "Charting`Private`Tag$11606#4"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{All, All},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.8737300176193027`*^9, {3.8737300596459723`*^9, 3.8737300665241117`*^9}},
 CellLabel->"Out[15]=",ExpressionUUID->"3524f72a-730c-44fb-9dc6-30e13e9dc263"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"CS", "[", 
   RowBox[{"e_", ",", "u0_"}], "]"}], ":=", 
  RowBox[{
   FractionBox[
    RowBox[{"4", " ", "\[Pi]", " "}], 
    RowBox[{"2", "mass", " ", "e", " ", "100"}]], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"Sin", "[", 
      RowBox[{"\[Delta]", "[", 
       RowBox[{"e", ",", "3", ",", "u0"}], "]"}], "]"}], ")"}], 
    "2"]}]}]], "Input",
 CellChangeTimes->{{3.8737301794533577`*^9, 3.873730202629624*^9}},
 CellLabel->"In[29]:=",ExpressionUUID->"1d64fb32-c3b8-4ea4-8f06-4b01e68febc7"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"CS", "[", 
      RowBox[{"ee", ",", 
       RowBox[{"-", "2"}]}], "]"}], ",", 
     RowBox[{"CS", "[", 
      RowBox[{"ee", ",", 
       RowBox[{"-", "5"}]}], "]"}], ",", 
     RowBox[{"CS", "[", 
      RowBox[{"ee", ",", 
       RowBox[{"-", "6.5"}]}], "]"}], ",", 
     RowBox[{"CS", "[", 
      RowBox[{"ee", ",", 
       RowBox[{"-", "6.8"}]}], "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"ee", ",", "0", ",", "1"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
   RowBox[{"PlotStyle", "->", "sps"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8737302193640547`*^9, 3.873730285086452*^9}},
 CellLabel->"In[38]:=",ExpressionUUID->"5bc7a9f5-c5f5-4792-852c-3f1f56196559"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0, 0, 0], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwVjmtUjAkYgLtRLjuJsCNlsxVFrok2vW9klcas7aZtNXQqoSVKy3K2Zag2
umGTUqhcYkeNdFmrccuYlPTNN9/caJrJNKX7l3NkpDPNtj+e85zn3+MYfTB4
t5mJiQl7gv8duLuHfNobBhLfnPP3SAVY+M07pGF6oAMrY/MbsQJWXysormZu
Qqf61ZZmE816rbh/hRmKpGMRL6RVAeGx+kcXmLHI33vhpL5ZAe3L//TLZibj
S7vtP28QKqBk9rb0s8xU9J1FitrrFMC1ql50hpmHhjDVj17FCvgpopQ//OEG
Si/XGWxiFDDNyI30d+Cj2UiE1eFlCnAOGup5nlmDIR7PyjPH5HC01QsLPv6L
hwynfWxfyGFujY3a1PQJMtidZ4cvyKFO69dCrG/Aqfyk3HCOHBrir9hFDAix
ZY6gPNlVDiMuPtYVBxqRyTE33B+Vgdkob4adczPa7lHqC0Uy8DRG1YVUtuBi
97DAjEsyaEn1aBW5ERi8JYD6vE8GC871j94ZEiO54Ga7u5cM+vUatWyFBIcY
Dymf6TK40VffzMih8PPxF1xXlRQY9uW0TifFf1ZzmvIrpJC1rzauxEmO9oVt
PjWnpPDJl6fbe0KBa+FIVmOoFJJXCWvmPFfijry4r391k4KmUphJOr7FIHZK
m7uBgrF4aZHr8TZccuBdNJukYMkGy/G4ehXOi7R27bhNgX+Z456quWrMlKZt
zE+h4Kj9roXvIjW49bfKi7wwCpYOhE1TsDvwi7NJvnEZBbOsh62EXzoweGFZ
fPYkCiyeFPG/LX6HnkkB2nCNBH7xqH8zhaVF0fCazoQ6CfC2RJc8/qRF5fwG
gTBHAiP982Z053filFy7pjXxEqgN8TfP26TDMpFlgZWfBE4EpF2r6dbhzDmz
E7ocJODYHqi8ntaFEeyRh+oxEr5bm5rHXNWNd/nvS9lyEiIdiFFTqhs59duF
t6tJGGnO1h7ivsfSKdNHz+SQUGCxMznYqQd7lIpC//0kzFLXlpkJe3BsXaPP
DBYJk1K4oV6JvdimmRuXuIgEkzs7+dqZfehTlP7HX5NJSLrMz57Z0Id+U4v2
cTrF8GhJhpVNQj+uDPIeWvVMDDaBSUtPMQbweIpDXWWpGFTuYuqYYADDztu9
6jwhhntX0Um3cxAz7rocrIoSA4fvaWM/PoiXJlevP7hRDCfTSrmqv4ewyoX1
4JOjGCp67IsYW2m0kr9UjpuLIeNKWpT6BxpVXVueHpvo2OBB58ogGm0DbSo/
molhvkBwjx1O4+7YxOI+UzGczd0hyoqmcTxGnlFhJGCvZ+GHqcdoTKsnuFF6
Aham2gZY3qJxba6XhOgmYHzd718pbtNYxvBY2dhFwJvBTsktHo2SspLixzoC
zoVXczZX0ZgjaCmo0BJgdAs+nCagkV8Y4ZKlJkBF5l41p2hM5mRfZMkJeJCu
j6FkNPZ7ba73kxGQ573L9bqSxkBDzKC3lADWzWU1G9U0ZtLNR5ZKCHj42+sm
bi+NbhnLHRitBOS7r8kNGqAxoeTzrcmvCUjUFoc60jQmWX+z3viKgMVb92ue
fqSRyUvPHG4iwMJUeuOcnsaTHsLve18S0FHrHR/1ZeK3NXK6tpEAQfz15SsM
NJ4u36Z5KyLg0oJpI0Yjjb5DRY+pFwT8B4Difmg=
       "]]},
     Annotation[#, "Charting`Private`Tag$17084#1"]& ], 
    TagBox[
     {RGBColor[1, 0, 0], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwVlGk8lAsDR5HiITtZKlK4CpPEkO0/4RpZcm0pISERkcstJes1bVLShFCI
VCRljRnRZKmunVL2LGXt0SJR8fZ+OL/z9Xw6ip5B9od4uLi40n7zf1seGm+v
nXAylm/iZCw024LXVO7YoKw2rmTrX/fNscX2zNSMElkzyMw9p1qctIVVc3fx
DVlHTAr5+dur2MLZe746SdYbGRF6dG/GbvRvPWuaIBsK3q8zhPIeGzAzZWqV
9/+DC9wy3EPbbGApkm/4JOM4POVn078I2aD8Y5MOqXAS5ZEGZQqN1rhUKK5q
pxIJ7rHaRD8ja2RJ2Z65IBsH1rq4zw16VghvTqnbUxWHkhPshWhZK+xhDPFs
2s+AkMA1nthFS6yeC45iZZwB5/p+F6dqS4R1McOmFc7j88+Mjj/MLWF3tefI
bpVL+O7tLf3Eaxdi+Ev+OC/LxHdDu8aQF3RsoV85V3yMCcXcC23RRXR0MIIm
ehuZiJrYWFt/jQ6lFeoFlOPXUE35zvjqSUfDr1xKV0cyoqaYlR3L5lj9JVlH
4eJ11IXbaq8wNUdq/ynTsqWbKP3FP4opM+xcty93wCkTa84OlSp0m2HSRXcl
f2Em2mZrWMrPzCA/fdKrUCMLwTcYrFvpZgjyoypz9WVhr1K4gNNuM4h7F929
rXcLzv/WLZpWmWLvvuyi2U+54JX293l7zwRkr6PALbnbCJqJe/MmzQRn3PkP
OZjeBo98qcRsvAlKvILkypm38WjmH95DgSYQDjRihOvm4dSNhlZPbRPUx751
WRVxB34+Z7Si63ZC677YqnX8+XgdrBH5+QsNLyj1Hs2a+ZA50Ttx/gMNHo/C
WJH78jEfOmmwvZeGhPLBY+/y83HQWX6hikPDeG1hb55NAfjl/SrTkmi4+WrX
o21X74OP2lsarkOD4HKMK12+CJ41li2lQkCqn7O7unkROoOzHbBgDKUudQ+x
wCLQjtl5j44aw/het1dvdRGM76q6RLCNUXkvuayC8hC6fgekxgKMMeMwzbZq
eYg5jaNOBzuM4HQ35b8QoWJEpJ3d0/7AEBlusdY+isXYE0gtsL9hiGHxo817
dYoxULGwNB5viKAIk1Yjt2Kw/M2EHP0Mcd5+poOvsBjd5yRYCsqGqP5p0pNu
XYLj8/4qgVkGULb7OP4svhQLrsELv3L1IRdF1z+cXQpfB6VhJlMfooVZ8QIV
pcg+x7qLOH0s8jtQ7IZLUWa2iDYvfbTVVoQM6JVhTfg29ttN+jilGbP0fawM
7xg6Ksfu7ECLiKQkhVaBmpeirk0sPTwzOnqo3akCifLyDpmFenjs31Ae6l+B
CXqKcWymHnIaw/axkitQFR9/9sS/egiL6r9pMVMBrguPl9dZ62EjmbfZK+0x
nFe06xoN6eJEyw6kfq0EX0wtXU1CF8KUsrEX/FXoDnHqyeTTRW6C5sUf66pw
xUKSo/CTinbrP966/1mFtZNXd+8co0LtP8kQ5eQq/JlBv0I8pmKwkcwr1mXB
4HxSe9YBKuhPbws3n2Ljpv/z6ehyHUiXig1wc9dAPGdzjkS0Nlwv/nAcXVOD
aJs7aarHtXHLe+y/BvUaCHwOnTcP0Ia6VGVl/N4afNRuLsjcqw2T4x7XpB7W
YJNjdGLNNm1QM5zzmFtrEbe3WfHB2HYcdM9zuztQi4klq65PDttRPmza1GrI
QfHBR87FdC28toj5lEDn4Gek+J1JaOHbgydrrO05uG8wGbhZTwu64TsOPj/M
gYZBG71KVQuPJbZ+q03iQK1TnqYqoIVKM7kNxeMc6P3jG0tp2QZ23mwIk/kM
ggkSiQ9dtoFz5MbafdN1SLJrPCsWp4mSCX9/9291iJ7lerQ+UhO5vvosL656
iMS4E5QwTTB8ul0CJeshFTmjsz9AExaeYmn/Gtaj/u8LvB8cNdG6jyHz4GI9
5man2LUqmuizOCq1QqMBjeYxIzVpWzGnYiRSeLQRF4tN1n6xpIBnoUB0rfJL
qOwaP538XA3UZY9yhwdNWM0zpvK1cTOa4rRbGra04uvp9Ub3xVWhkDi1cO9j
G+Kf2jwtDFHB1PzgwCvNDjx8mmV9Ul4ZuZOsl8KXOlHSGwva600QXn+HHB3t
QgOtc4dXxUZc9CvzyVJ6DaJmdii9QRHfaAWjvlHdGFwVNR/QswGhWnWla569
wdirncFG4hsw+KAuvl2xB55smuN0qjx+HOlK33yqFyMVmwTcqOuhtpNvyYfV
B+KDnYbMqnWg31I8/Eh6ACbKFtmVr+RwYv2Bje9cB1Fs7yPZ3yIL9WknwW6b
IaxN5D/NOygDCZFZ/rrFIcSuHGJxpqXBW5NetCnjHbh+BrT1yUrDX5v1lrAa
xnZHwcPbzdagYJdn1pNvw1jZ7NUvHS2FuSk50ffJI3ixQSIFXZIoc6CvYJqN
oj2ekjuiJIkoC0Zm6ftRNF36kp17QQKK/ZZvchhj2N9paXb5hzj0deOYslrv
MSJ16mWBrzhc5VsXuDvfg3t0ern/vRjmXiYMH4v5gFivxLSgg2JI5XUPtVca
B/gMOGpTopAYKLvFUzeOGnWLvsrTolgZEeO4I3gCTrBbmycgCq577kXD4pMo
/i6a6mopgr/TihLEOZNwv+rsmBImjGq1c/xigVNgn9Vypf4eqpjl3+qxwtOQ
ExET5F4phD6Nts6T7GkUZsR/2WK7Gg9vQmnUfQas50P9V+8Lwq2IKrZ+aQYj
rTP+URKCiGZkx/Tlf4Rq57Uc0VABFI6vTxe2JjFvOX5TfYzAuRsMj4HdJGiv
lL5XDRPwtp9RfmBHwmv7Z9auIQLr2OyHNs4kmtQ6Y317CVy4vL/hoicJLdsq
gYJ2Ar7U658ETpLo+evGlOkTAhvjJC348kiosS7p1qUQWNI7LdR9l8TTacWi
A9cIvJ0Z6cgrIPF+41zCjyQCic4lbuaPSFh2fYnWvkRgeYt9CINNorYsYkdB
HIG+9ss3V3SSoMgZDBQEE3h8Zt6r8xUJo8XyapsgAkyDA5tz3pCwv93NMxtA
wOo2pdRkgIR5ZXa1ji+BqrDmFzETJKIMr36udyeQrKFz2W76d29G0M8jrgSC
hzMcFUkSAW4q/qIuBFStAwZrv5K47/fkuKsTAV7urtzEeRKOUUclVjgQGCoz
OOKxSML1fJJk/l8E2Edytmr+InE+RvKU3W4CKQqCc8vLJASP9JguWBH4H0gX
d+s=
       "]]},
     Annotation[#, "Charting`Private`Tag$17084#2"]& ], 
    TagBox[
     {RGBColor[0, 0, 1], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwVVGc81QsDVtybVPib2ckoo5Pr5pD12MSJDkqEXKsiVEooSaEUik5GWVmJ
7D2SZJRQyLVnchH+dMOVyuv98Pyez88Ud/Qyd9nMxMT0aAP/Z2OXqY666aOa
NzctRvs9twSLruC5EYEDWH/w7oxXtCX+TI5LKBbQw9NaY4FMX0uYtPUUJQpY
Ymv6U51OfUtYOa+8iBZwxkz5eOD2TxYY2n9LN0LgImSuUX1pey3ASN5ZJ3Xi
Esyt0gymCAsYc2Sr1yb4YPWUVkf5mjnK5luVSDE/nLU1sCTfmyMyl2svXfoa
PPqY6Ln+5kjhNQu9IxCM9itKeb29dFxpi204VhUMEmb8To10HAsZ3SxxIgRH
9ph68RfRsX3pfGB1QihqVaw71+7S4fuR4TsrFoZX063yGdp00B/0u5lKR6JC
/6w6R8kRBLEW7wkTYGCKnvgzytkMsoZRt4vOMeCxrf/ngK4ZOkO8pgeaGVhk
i6WdkTCDJLN8DsXnIQo+GiQafzJF0890ysfOGHzqiqZyuZpi+78xSmLh8Wiw
lY7nvHwYcUP+uqW/kvAjtOykYr0JtIWt04ePJuOWfal4RJYJZmyUf2PNTUZ8
nE0G3z0TiM76OeXuS8FEf/1Ssa0JvM5QpZgGUyAgbaPG9t0YXM75WRkqqUil
xfbXqhvjuPWT/IXFdIzK8jKcRoxADliypQpmwHOreYtQqxFC7VldLHQzQO8d
8V+sMEKxk5dgGSMDgosSmt+ijcDuqRFyRTkTa4pzYkeMjNB4o8/m94CnqDni
aSdbYQjF58TvwqzZKEx/aNKdbYC3lEaHNoVs9DIiDrx5ZACHQt/qa9bZuJK0
NbjjjgEiykbOjWVnQ0yBJUHa3QBTdbkDmYdzwM+mU9Qtb4Ck7kOFfzx4jqFA
R2HxEn1sWw+yNRTNR+HZFXV6hx7izljZyxvk45vDtCFzgx4kP8o7EJ75OMt8
k2gq04Pmsx6ngRf5yO8dcAxI0EPls5jSckoB2i77J+Wf1sOcxWyNSXsBnOku
F1hZ9HA0K/ad944ilPY5P3mrr4sEuxs0V/EiXJLLnjyhqotxLo+240pFqBKs
Clqj6MIrQOe9hl0Rkj/TbB136iLMfK5zS24RFM/5SLZ+0cGLHzr9j2nF2CM0
VvyDoQMp+vzU67sl4C2MN9Je0IZgoKHqqSclMDlwQ/zaZ21w5qbcZSsvgYh1
aE5zvza+s1pQ6OMloFQG+l1v1MaHunLvYZXSjR5Ox4Q+1oa/QtCv/z6Xgv12
8qHOQ9po5+DhoWiVwysrt2wwTwuvNTxcOo6WI0rT+/utDC1UuDeVXXQvR9Lz
qQCNBC2kNftaV8eUozoxnak5TAu+gUNJRnPl6JfPkDZ30cJuMlPG6VEFXtVH
fFUV1cLl9oOI+1aJElubgP59ADul9PNb1iq0cj89+VgUSI9QCF8TrgKaV8bd
OIAO2p4+e/0qOHTW+SotakLuHY+3VEwV3ItPJfmUamKkmcwsUq7Gkfir6rFa
mjB8lcHe5l8D/xdrD1odNMBfQgxv2vQScb5LNw81qsE2fM1ygu8lolZ4k8Qq
1JDq/Pldk/xLhPBbLWzKUYM8b2Xl3eMvYcnntTpzXw06Pg4PeQteImJhfwvN
Tg3UBKtMxv66jZ2uy2r/p4q/7DPtsobrIHZLZ/E6VRVl47qt79XrURp5C/QO
FfxtFLQYYViPTtG3Xl+bVLCcV8tHM6+H7rzszoQaFShfOfjXm1P14A349ylr
lgoquPcv10XXY8W4NkM/UAWVeoK7iqbq8VNo1sFgvwpqMhe8GYzX0KeQqrse
KqPeLVHIerYBxhG7dpj5UFE87e5uv9wAau6VqeOeVKSfVq12YmqEqcsW/dOu
VIS49th48jRi6ecXt9hjVBg5Eo9uqjci7J5mK1WZivfWITvzwhvhNlLHeWlV
CYNGHrzM+5rgyb4Y2BKkhCVpDY5cj2bY6ETuako5gM2rOZxCUi0ITQ74dWtR
EdR1hzKLvFa8KVz2S477A63BB9qbZN+jZS3y6th5BYjd/7L6bP4D5KKuu5wp
puDLyshwt0InFqZY57ZMyiN9prqFPbILMRJiP1Jl5MAu8pScmPgIUUfrB2/S
ZRB+ptQ1RfJv9E+zjhH2e7GslTNxOrAHMtPjvjSRPbio2FDC97oXchShJIUd
0hjJa7jbId6PuK8d7VJMUlhz+/hYxn8AM6B8/cAuCTntLb9cqwchnDi58bIS
MEwVP1XIP4zg0pqTMud247LIyd1jtiMoEilOnqsRh/zs0W09h0fx36gms5mA
OLg5Flgbvo/iH5mdz1e8d4Hl5eN8iYQxqOx74jzfLAb3A9V9W03GQWWrjGRN
F0XOIceU2uVxqFQc+3tznAiWvghyTsZ8wp9P+ZM5c4VRamHIzNCbwNQ731S1
F0IINApJLpmcwI+SRA+lGUGIDxn3poV8hh1PjIqrtCBUlYMZAoqT8EmzP2Fh
JwBb0ferm7omQXmWqRGbsxNLLRHj54L+wcC7prs93/gRx2J/0VxyCon9zvJd
5vzgHi5N3dwwhTb2HelXq/nwW0CQ5cHz00i9byggIcYHpmf2+eNcMxBqEuVk
Y/DiwqP8CK76GVit+2iqsvDihdxtVsLzCwaPs7RGBfCAML4gf4N9Ftw3Li2e
+o0Hg/s+dPnVzEK2e7Xw+W1uFCRBcsJ+Dqs/bFR3CXDDLp9KiPyaQ2/E8R7e
Ai5cD3kSNJg9j5kSGl+ODhdyp0Qes9NIcA/1N8WMEridGOIwbEqizBSHp4cJ
OJvPSeXRSSwvBwurDxEQrqkpOGxFgmOP28PxPgJ37p1oCnckEcYcP6/YReA0
NX6RzY/EhYlgylgjgd3BPEZbMsmN/CpcvHMI/FK5uqMni0T4tR8T7c8I9M19
6szMIXEoeiRbNovAfatiO4NCEjxqamyf0gmsy5p7h9SQEK6w4bVKIjDYcS+J
uYvEsQ/Dh2lRBCpCV5y6uknIK3dQnt8jwFA7KZPWS6KJx/jo9kgCJhmUEp1h
EqFzYRfa7hCo8m17GzRNQuuU/ax5MIGYfUr36LMk1nJTVkpuEDg/nmApTm7o
vTBiyB9EYC/t7EjdNxL9/+woHwwgwLLpY/r9FRLs9gPduEpgtFTNzeH7hh8j
DgfT/AnUuKXtV/hJQto2Yux3PwKxYtuW1tdJ6LeatbtdJvA/+v6BYA==
       "]]},
     Annotation[#, "Charting`Private`Tag$17084#3"]& ], 
    TagBox[
     {RGBColor[0, 1, 0], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwVxXk8lAkDB3BJmyg8KM241jHkmqgwUX7OJiG5thSTig4dypGjYi2jFKt3
m1CInEmMY6bCFJVYFIlCNIQswsOWRMX7vn98P1+1gwGufqIiIiIh//P/d/iN
tNWOeliMFM1eWfmOCTEb6uk+yiYEePZfjmlkYmNmanoFxRYVPwzUiEomHF52
lmdQ3HHR9fm+gFQmdvvOPvqL4guTSynUjXuYeL/+ok0iJRhmSstel/VsQ9Zq
57jLlFjU3A+2CJq2w7mXKXW/VcWCTX3gfWLQDr+x+0U19rFxfnQdLfSNHVbO
nImqTo9D5L0R7SeVdgjr4ISNq8ZjMn/+pEaMHVyuvfPfqfUnJJ9uiV9GsUO0
eIV2PIUDPdIoLGiXLXSZ/7lUfpoDjXNvdXxtbPGaHTDa08DByjj+fj8TW2gu
1S+in72OY5kBPalKtqj/mUvveJ0M/9CxqtpRG6z8nGysmnAD9owB5/k4G6S+
j7DhL9zC7Ys3BO4vrGGl5Jkr9MiE8ahTn8oTa4ztNV0mXpyJ+oVS2le+NVTG
ww8VG2QhvP/yk+ZMawQcM6GJ9GZB0HjokUywNWR9uXfyGNlouZt2b1LFGns8
b3OnpnMRbzR05vgFK5A97hLZ1DywxHq6PIKtEMcS93OzycPa7q8x9setUHEo
gHqfk4e9thqnHD2tIHVqK/ucaT68pYxW1Rhb4fkf3Xt/uVAAe8niAocpS2y4
R/yiJH4Xmm6i9MLjlpBcjPZiqnBRsIZfOjVtgdRju1n627jImnCUnhq0gGaH
vg9xiguua1L39zcWsCjsPNTziAuepUkAo8oClYXJ/Af0Uvxs3tVlE2OBCbdx
gUNLKdq91f++rmABjzspzUGryuFCrWgSbN8KmsvkyLMrPLhCtkOxyRzUKKbZ
kds8KGnJki2PzCFTnHVF4gEP6XbzpQll5pgXd6O7DPAw99XrX+0b5nhV+yBI
yOCjSCdhsueoOSIMoxe+feRD0VfWiC5hjhZpeXm65QP0rI6N1HA3Q2jLZqR+
qYSQdV2rVXQzpOj8j43iVTjCqzwaOsdAbqJhwnelKqw4cfGc1hQDbY7a3Sy7
KjwvOziT+p4BvWb5IFpyFR4HOp5oechAXwOZX25aDTO2uGLIaQaYT/KkXkYI
YNZ22iB6yBQKPEK4ZEkN/h5krTAUmuD+gM2L1i1Pkcbuymoe24Sn/hmKnuN1
8FE9udfeYSNmtLZKF59sQODya684P4wgOlcko0hrAr2ezFavNYTJos99t5IX
UPDLClIOX48XsZta6nVbMW3OLRflGUD16qe5wslXuLv67PDsDz18mu0TvjF8
DecPIyrGibrIHatukvqzHfu4D7dH03QgpVxADg114MzHgEQMayPhGP9wluZb
sOIZ7NK3WvhqWTR0NKoT0dWbdhGNNARvqOOtedaFap+0pPkPmugrqbvSpvYO
6ZXx1a7fNfDdvyNNJ6IHpWmyy67qakDPavnC4epeHLj0zKQ/RB3MbLUjZQpC
uFPHF0pr1BCqvF/9g1cflDLOH6CrqUF/3EOy06kfbZ6PH6ck/Qo56Snxuvl+
aOkVFBvPqkKsJo2rkf4Bvdeb5xY/q+D4puruFQ4DCD284dboJ2UU2R/Mevx1
ABrFmXIbvylh5hNVZjh5EA2Tv4qnqCmB78ZcyrEdgtPNiJG/7BQRtZ2dyRse
wrQwUcUohgq19zu6ctgfQXDK5de1UGBmGsuhbBiGyBZtAUeOAi+V1rkl7cOw
VinQVzi1FjNNiQOno//B+bkIHuuZAlLFWMGumiNw+/LySbyeAuSE/GzRuhHs
c/dy3ZmxBssuRLtvPjOKe/qtOf+uWAORQhZ3QHYMZ79F6kXGrUbgTW6i7NMx
iJwMOWu2II9HepfEiVOfsPOfOV3nSHkQOwL1/5Aax5rIkqlBSXn0GrxqDxeM
43bmT+d3HDmU3oLmEGsCnhz7Wh1dOXhzTQjlhQmsPXpAPOq5LH5n347uvTsJ
S4Yp7aSXLIpHlNOkHEmsMnAekl4gcCmD7SPcSSJ5l45T4A8Cvq4TtBIXElHq
w9865gkoCQSlTrtJUN0b1qbNErictK8+4SCJnmy/MN1pAkdNbkxLhJN4GyJT
6DFEQD1WfvvyfBJh71QOvG8isMA4v6rzDom+zsDl2xsJdE8Mvs4vIrHalTtZ
3kDg6u4K721lJIT1OVHxdQQWdV2D2AISrRJNhZsfE+htS7q1tJ2ExNtkIr+M
wMO42UPtb0jETG2ukS8lwDHfr5PTReLHcOu9mBICDnl0nrWQRG5hCvNgEYGq
sJeN0aMklrImw9XzCCQbGCe5jJMwkP6ddi2HwJmBdHc1kgTSr8iIZRNY53ii
r/YLiWrNjtzhWwTElnTkXp0l4Ub57LQng0A/39zfZ56Efq+zSVMaAYF/znrD
nyR4IT3eW24SSFGVnFlcJCH1Oa6uJJXAfwHuv2rW
       "]]},
     Annotation[#, "Charting`Private`Tag$17084#4"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{All, All},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.8737302421940727`*^9, 3.8737302865099287`*^9}},
 CellLabel->"Out[38]=",ExpressionUUID->"faad53b7-2ecc-458c-bd43-c2e2528af9c2"]
}, Open  ]]
},
WindowSize->{1389.75, 768.75},
WindowMargins->{{0, Automatic}, {0, Automatic}},
FrontEndVersion->"13.0 for Linux x86 (64-bit) (December 2, 2021)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"a31d9556-d2b9-4444-bce5-a086ed60159c"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 1110, 20, 71, "Input",ExpressionUUID->"4a74809f-4d6f-4bf5-84df-1c81e845ccce"],
Cell[1671, 42, 1195, 33, 51, "Input",ExpressionUUID->"4abccf8a-d47b-424b-acb1-36e9461ec53e"],
Cell[2869, 77, 586, 12, 50, "Input",ExpressionUUID->"ddda9570-52bb-4f7f-8d3d-86bdfd8586a4"],
Cell[3458, 91, 654, 14, 50, "Input",ExpressionUUID->"86a69de3-cbb6-4ee5-8cee-9244eebdf060"],
Cell[4115, 107, 1856, 55, 49, "Input",ExpressionUUID->"f9b93666-a5b0-48b1-94d6-6d1a85496d1b"],
Cell[CellGroupData[{
Cell[5996, 166, 950, 24, 29, "Input",ExpressionUUID->"aa0fae4a-3cd8-46d8-ac13-32fe6f27924f"],
Cell[6949, 192, 23640, 414, 230, "Output",ExpressionUUID->"3524f72a-730c-44fb-9dc6-30e13e9dc263"]
}, Open  ]],
Cell[30604, 609, 543, 15, 45, "Input",ExpressionUUID->"1d64fb32-c3b8-4ea4-8f06-4b01e68febc7"],
Cell[CellGroupData[{
Cell[31172, 628, 787, 22, 29, "Input",ExpressionUUID->"5bc7a9f5-c5f5-4792-852c-3f1f56196559"],
Cell[31962, 652, 12950, 238, 227, "Output",ExpressionUUID->"faad53b7-2ecc-458c-bd43-c2e2528af9c2"]
}, Open  ]]
}
]
*)