-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathAreaofEllipse.nb
More file actions
1204 lines (1165 loc) · 54.1 KB
/
AreaofEllipse.nb
File metadata and controls
1204 lines (1165 loc) · 54.1 KB
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 12.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 55236, 1196]
NotebookOptionsPosition[ 52514, 1139]
NotebookOutlinePosition[ 52910, 1155]
CellTagsIndexPosition[ 52867, 1152]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ",
RowBox[{
RowBox[{
"Area", " ", "of", " ", "an", " ", "ellipse", " ", "using", " ",
"\"\<polar coordinates\>\"", " ",
RowBox[{"integration", ".", "\[IndentingNewLine]",
RowBox[{"https", ":"}]}]}], "//",
RowBox[{
RowBox[{
RowBox[{"www", ".", "quora", ".", "com"}], "/", "How"}], "-", "can",
"-", "I", "-", "find", "-", "the", "-", "area", "-", "enclosed", "-",
"by", "-",
RowBox[{"11", "x"}], "-", "2", "-", "4", "-", "sqrt", "-", "3", "-",
"xy", "-",
RowBox[{"7", "y"}], "-", "2", "-", "1", "-", "0", "-", "via", "-",
"polar", "-", "integration"}]}], "\[IndentingNewLine]", " ", "*)"}],
"\[IndentingNewLine]",
RowBox[{"eq", "=",
RowBox[{
RowBox[{
RowBox[{"11",
SuperscriptBox["x", "2"]}], "+",
RowBox[{"4",
SqrtBox["3"], "x", " ", "y"}], "+",
RowBox[{"7",
SuperscriptBox["y", "2"]}], "-", "1"}], "\[Equal]", "0"}]}]}]], "Input",\
CellChangeTimes->{{3.8748063022322397`*^9, 3.8748063295996513`*^9}, {
3.874809168226102*^9, 3.874809212839238*^9}},
CellLabel->"In[1]:=",ExpressionUUID->"92135a08-6b08-4843-895f-17c9841b36bf"],
Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"-", "1"}], "+",
RowBox[{"11", " ",
SuperscriptBox["x", "2"]}], "+",
RowBox[{"4", " ",
SqrtBox["3"], " ", "x", " ", "y"}], "+",
RowBox[{"7", " ",
SuperscriptBox["y", "2"]}]}], "\[Equal]", "0"}]], "Output",
CellChangeTimes->{3.8748063307260113`*^9, 3.874808620618083*^9,
3.874809217255077*^9, 3.874849217739643*^9, 3.9769453343847113`*^9,
3.9769466789148417`*^9, 3.9769474908863773`*^9},
CellLabel->"Out[1]=",ExpressionUUID->"69eef7c4-d0b1-4488-91ac-e2725a826c66"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"fs", "=",
RowBox[{"Map", "[",
RowBox[{
RowBox[{
RowBox[{"#", "[",
RowBox[{"[", "2", "]"}], "]"}], "&"}], ",",
RowBox[{
RowBox[{"Solve", "[",
RowBox[{"eq", ",", " ", "y"}], "]"}], "//", "Flatten"}]}],
"]"}]}]], "Input",
CellChangeTimes->{{3.874806350986723*^9, 3.8748063704091063`*^9}, {
3.8748064132781477`*^9, 3.8748064740282917`*^9}},
CellLabel->"In[2]:=",ExpressionUUID->"f05dcec5-cae2-423c-8918-6a76b17b348a"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{
FractionBox["1", "7"], " ",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"-", "2"}], " ",
SqrtBox["3"], " ", "x"}], "-",
SqrtBox[
RowBox[{"7", "-",
RowBox[{"65", " ",
SuperscriptBox["x", "2"]}]}]]}], ")"}]}], ",",
RowBox[{
FractionBox["1", "7"], " ",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"-", "2"}], " ",
SqrtBox["3"], " ", "x"}], "+",
SqrtBox[
RowBox[{"7", "-",
RowBox[{"65", " ",
SuperscriptBox["x", "2"]}]}]]}], ")"}]}]}], "}"}]], "Output",
CellChangeTimes->{{3.874806356388144*^9, 3.874806376468141*^9}, {
3.874806419149021*^9, 3.874806474603272*^9}, 3.8748086219052153`*^9,
3.874809218754685*^9, 3.874849220349002*^9, 3.9769453361769533`*^9,
3.97694668009085*^9, 3.976947491859004*^9},
CellLabel->"Out[2]=",ExpressionUUID->"6ceb8085-aa7a-4e0c-b101-a582e1d155dc"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Plot", "[",
RowBox[{"fs", ",", " ",
RowBox[{"{",
RowBox[{"x", ",", " ",
RowBox[{"-", "1"}], ",", " ", "1"}], "}"}]}], "]"}]], "Input",
CellChangeTimes->{{3.874806332124896*^9, 3.874806343242053*^9}, {
3.874806379531989*^9, 3.8748063939931517`*^9}, {3.874806478236711*^9,
3.8748064823951693`*^9}},
CellLabel->"In[3]:=",ExpressionUUID->"5984b411-9c8a-445a-9602-288db7cc966b"],
Cell[BoxData[
GraphicsBox[{{{}, {},
TagBox[
{RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
1.], LineBox[CompressedData["
1:eJwV13k0lOsfAHAlWhBJSsmWlC5CpRJ9iyKVJVu2qwhtyK1IkiUlFMmWa5dd
WZI9fGfGvpuZ15gZoSgpkqUs3dLv/f31ns853+ec5zzvd3keWcdrps7LeXh4
5Jfx8Pz/2zRu+8Cbpxe3jbOMXR0bwPF2Rp7PHIF2V/7zqW+kwhL/eKfvBIH3
DS6UdzygQEKM+ozfMIEL9zn55lsQNOR8xAPZBPb//bl+0bMWGEVUzaAuAhVa
td8W8dSAu/bqcw8aCBR7dW8mo6Aa1rSbBD2sJtB2cuZUws0q0B0dan+UTWCm
+kXJV1oV8O7GjqnwJAJFHZ932uwvB99l18QiowjUXT03IHO4DEoll+xi/Ah8
e8iZ9vXya5A3k5xMsiTQd0JQaZtUMVDeXRBNPU3gwU91n6ilhWDn/kIjXYdA
gUTtiXP7CyAmVNM/S4XA1X4v23is8mEFxUqkkJ/AzjKmsOlcJnxUilXHciZu
TJuJE0qMg11+WZNjL5j4ISvxdsyuGHDvLssXTWeiYrpsjHf9U1j4hyXn8oiJ
Gs9rV7tJh8Gaio3rhc4z0SA3dlWejAeoHk2YtVnNRB6KdN3d0FD0jMovur/E
QNGBdu/2kHCsHqm+WjjLQNM3ju6DdZGoG9w/smyQgbZn2c8HM2PQsl2SyC1h
YOLXDLXRA0l4xyK1dM6WgVuT53uNDbOQkl3kIXOGgQWWuTYif2XjigVUOqnH
wMW7/oubVuZgRMK7zGRVBq7oerp7Z3Uupg3Jxh7jY+BcsO+GPUIvselypmdU
IR0fX+IXY2m8wnV+eRoqy+g4PnrXW1q2AoPLd7aIzPegkaC0r4ZnBf6czLGe
nejBK+q3ZoZaKnD4XLZvJbsHN3i3vP3mXoklRzNoOsU9qBE4YqtYXoWmfMlG
Z8/14NNM5uJG9RrcISHbyXOnG8/Qk36k36CgksY3pazZDqyJovvFVzWgi4q/
f9VAB8r5JlGc3jZgqoIwo7O5A9eHdWWu/dOAoht335pL7MB/9h2//Eu3Eefn
3an6xzpQRyp2v35HI9KqJi2+xLSj4gYDbTqrCa20J/12a7Rh2Ga3ZJ3eFgw6
9rWnyrsZa2KVTHMCO5B2YsJ474VmfB4kf3guqQN5DMe7Cg2bce/PnGW7KzvQ
1+JzR4ZcM/qln/K1n+xAL5ePLeEdTUhduKZYad2Jl0IGKA4yTTjmcqTMVakL
T3d0Fq9uacDnVS4jfrRu3GBeGGm9kYZRtvyJhA0dSxQeVR/npWF7eMUOQ2c6
mixe/KD2jYrxTK99ddfoGJYie2BNMxU91j35EnqfjktjMUPVXlT0dDO3aCmg
4+jdOyqSLAqu2Lai9+MSHd2Nn8bXKiKK3cg2aEpkYMlco3CLfDUGbjqVJUJl
YkT+c9rhxSrcrDf+LKyNiVft/T3LOqvwaNjCr19MJso3Heh/7lWF+jtOqrSP
MjEu9kX23ZZK/CAmlygkSKDP3qfae1wrMNNqw5iyBYE61+2uprwuxdjhfeXf
3xOY1f8x/PP9QmRYV+npjPYi11c8NXlPIS63/qng8bUX10rpF58ZLsBtr/XE
4r/3otf5XEYVFGDAuJ0kazkL9T5e2Rj28wUGSo4JCUizcGzyW9oujzz8wJdM
tbZk4a7lP19ftc3E+guurdEUFtpn7GqUXpOJISFNvPzNLIw6ZstiVmbg8eKR
m9c7WfgzuGbhkHgGfrS58ucgl4XtAv7agvR0rL3e1hM0y0L2rlUR1uEp6PQr
cTfv9j78+4xkaf+FGPzIcTIrD+rDidkfBz0bolFJfEbPPLQPfeO6ce32aMw2
G6RORPRhUv+99qOjTzFGLtVYKKEPS9sZjcb0CIxq9mMJFvdh0XfXbWPc+5g2
aT7C19+HAcROYku9NmqXHq9/r8RGJeGBK3zK8sBf9UaTrsZG1vS9VJGFE+Ak
TRx5o8HGO9F6bVf3nodG04EFvyNsrJcJfZgsch0O2UhubDRj48WLu1zKl+5B
o2rM3KnbbKxac6eX5hkNoQq/6DoUNqYpb2AWycVA9eW++aEGNgabFPUkdcdA
x9aFl16tbDSPG2n3UowD1odYWiSDjVOyp2m7BuJB89p0yK0PbFQ4uLUo6lgK
JK0+Mbh2JQdbX5879io0Exypwjfs9Tmo1ahUv+NdJvBEcyb+OsXBYtbi0RSN
LMgzEWHPGHEwbjEaHo1kwQu7Z79dLTnoBC2azto5sOJS+DlRZw7+aVNVk5jO
g7tpQkdkAjm4b3iZVKBVEWyVtDqqXcpBleqz06p5RXDqIdWWt4KDClEFDe8W
i4As1qj6Kg4KpL/cl2ZYDBIRFkbqyMFfJlZv/0wWw/KSv7tbWjk4UFykiGol
cOhqmgt7gIPcJJv1G21KIH6Vyg14x0FWCN9v93slUKb60y9tmINd5217pJgl
cEXL6LHJJw7Wiay85XfzNRT/KT52aYqDqR72DdoVpVBXdtLk5jIuJtmtLowd
KoX7fkoWPrxcjD9R+uzryjKw3Van4cPHxUiZNVeTrcpg/drA4IuruRjQU7bu
92IZWPjLbFgQ4aJvzfn/zOXKYZWwXk2HKBe9cwU+vjxZDktfEpUSxbjoEeBQ
aZdYDho/0u23b+Kig6rQuRpyLs6oquGkFBftJStPiDlXgM166YtPZbhos+qC
umt4BdyM2FGvLMdF03eVfJKDFeA16f/YbDsXdSKdXt7xr4Sa1LgNZn9xcdvU
m0VNWhW8zKXelNYg99Ope//71yoQ4grqeO3nIiu/fW2hRDUM/i1wvfEAF8Oc
+7fJ/lMN9yJaP58+xMUp7k/DlTJvYN3OsPo/wMWTlffYlFNvIG1sXk3iKBcz
Y9dc8Ln1Bk66m7zfpcPFsyabvb92vYEOuuC0xjEu1jUefM68WwPrpXUsp/W5
uCmDqhSeVwMFtGHDjhNcvB5gUKHXWwP/sNkOqQZc3H7IuqNKqRamPRPalU9x
0W/Te8sbVrXgJFyy8TNp9o9L75Xu10Iob4Vr0mkuPi72nkvtr4V03/it44Zc
HA3nCbReWQe6En9fCzbi4pGrIQLr99SBz7KhMgljLiacEInrsK8DEYGZ8UzS
s9vjZYLD6iDgm53gThMunuaVeQHldbDhL866DNLZ73L2Lb6vg5B+0fkNZ8jz
T6w46XoQYZnbX6dGSJd6Q+92Z4QRX6kibVMurrVsPjcUiTC7af7tE9KX9hh/
ia9BUPKX7OaQpor03TQdQ6h888l/ixkXhW4V/eFbQwF2UMZnc9LVzrMNw9IU
0BJ5zB9M+qL5/jDcR4H6ua6eQtJiuneMk05RIIv3kU43aaoait12oAC/aI/F
J9LuMrxci1sUmBlJWTFPeouwfqp6OAVeLo6dWCLd8jvMSTiDAm6XIiV/k/ac
6FKcqKTAjeKLt2dJy/WLfmvpooB7lJr1e9JdrZalWR8ocPhDWlkT6TuVCbfv
/aTAGmerkAzSO3MGD58ToUJK5Je2W6SJWLkVWgpU2H6d46VLOvC+S+smLSqI
esU94CetciM/4scZKjgfujdFIc+j32HSjHGRCipOE8XXSYeYqEsU3aXCMef1
lC2k94HX4KNoKmhZyG+sIc97WLk641IeFdKqNYvMSD+RXLp0HKlwwvp84DD5
v7QEdFTkeqkws/Q4+BLpz4sPZpe+UOG3UXbFR/J/x421Vvbz0MAw3F/IhrRu
n5BfpTgNzlziBDSS+TLVeEY3VokGijSHNTtJG2RwOo2saDCubKnDIPNt6kGh
S6MbDQL4pc23ko6/FPTnUBANcvOYr86R+TmmrKymWEiDqm+te9rIfI4UWd6W
Wk8DLfujBjNkvh+YZTmKc2jwoFY4fB3pkMqAGN4V9WCpYRC4T4+LqokWyj4S
9SD320hO8ziZ73d3NU2p1MOjuHeMfWQ97dQl5ges68E5d0lHlKy3ls4dtpVF
9WBX/5+vvTbZT4p/zao01sNPXlrlZi2y3qLpj7O49TCTmUzr1CTz0epOXRRf
A/ya1ZTaRNb7qpEuWTfbBggXi/yjpc7F4saM6hGPBrB2kKooVOWiVa63mU1w
A5i6F/0Q383FXDe5B/qvGiCAeDDcRvYX/QXPMdmVjbDNRcED5bk4yT0ZGC/Z
CA5X4i0+kP0prlZ6s7B6I8TKNN78Q/av0Xutp37ZNUJcRsH+TVu5GCy0tYhV
0gi8L+QKGGQ/VJma1jdsaYSvw08O55L9ksVoelc/0Ahq86FlN8h+qhDvIfpq
VRN4vGHtmBDgYpNcg+cj+ya4MND0zmE5mf98/65dfrMJ5D86On78w0HxMbcc
79AmsHV46Wb/m4MuBRs5zqVN0J40MKW6wEH+A1e1jqxphhdt64wtv3LwuKEo
74+yZjhkFr3lWS8HIxxWTy+2NoOPtKN4JoODbE+eoaWBZgj/9upITjcHr6RM
Vq/ib4FtCpnfn5DzKOJb23VJyxZQNxMI+VxDxj8NGtada4FiZyvrknQOymb7
9JxY0wqBAaZ3JFPI9dX/1BlKtUKKwcs/dxLI+TZ8LuGsXiskxHHnNkeT8Xu1
TK/GtQJNXXp1XxAZz/pOjdJoA993dC/NCxx8/WW86NnJNpiextpZe3L90nBy
kn0bvLjhcy3dhtzPDvrt7OA2uNIc1MU5Q8Z7F6hVs9pg3dq2xzVAxm92SX/v
1Q7qTW7EQQkOPjnXF6Ba0QHZrw9vvk8j7xMTWjDa0QGddebC22vY6Hf7+e/E
4Q4w7FrULStjo3u0m8/KtZ1Q92xo5HkuG42aeW8OOnUCJ2GrQlE4G9fuVr34
WLQLtAKy5q+Ys5G/Jna7zs4uUIihVkyeZuPSif9G5rW7IIbINj53nI2Tjk3n
L1zuAoM35835yftRV5ydjSalC06WPHnuIs7GiKUQwzG3bqC+PTl+jNmHwY8m
BVKCuuG7sb67f1sfOS/M28z+7Qbx05mxGdQ+dFeX1seGbgiPGhrKKerDHa9t
ZNqFeoD/k13Av2F9uPLTlMRzpx4wun7f4uOhPmw22ipkIkqHkjmp8c+hLBRT
eS9+XZIOBtm3jsT5sfC8UJZMjAIdOOM+pdtvkPfPdqW9bE06oH+47YQtC5UM
tGzPX6BDZsoQVVCJhZG6tnkepXR4tP/aC9WmXnwrJ/06qo4OhjKfDNiVvai4
fKSmtIUOl+VNZG1f9CKNcqVn4S0d8m+x5r896cUf2j7zAXwMcHS2Cdc924s2
B/49/tSSAZ+e9Yl9HiAwe+Pfxq/PM0DZyOfTRCeBs3My1r1XGKB/aLc41hL4
uCzXVcKfASJxMfvbyHc2qldGp+cwIEp9cdTGikBBUd/k+lcMmLrhvP+VHoFW
05Dz8Q0Dvva1C3bsJXCmqLlasYcB5wMETA1ECJRX7nv/aoEB/RK6pXsbmOgh
mDjOXM6EYo+834HFTKwdt//xQ5AJYl1bpx4kMdEyf3SVpiwTVuf5qOTdYGLI
jrndNAMmKMqILd8sxcSvshv8GAlMkKldJynowMADd9L/oWYy4cKYVeUPfQbe
I5SdiwuZsNN7g5qzCgPFH+qdjqAxISE/3YqySMfDE7c2n/xCrr9d8LMylI4h
x1esPfidCcaBUl/lXenISIlctnOJCcpPDhYMGtLR5UzeGJ8oActEqdtEhen4
pJxbTj1IQEZUfJrYwx7kCl/ML9YlQFOWqI1w7EH5y7PJqYYEVPsZuA5q9WDl
FsEHdx0IWHSKjN7d343LPeO9Xa8SIFTqrh7ysBtPd8m72noS4BZwtThKtRvf
BWibHQwl40VT4uBWF+7iturtjCaAZ1+i5fT6Lry5x1JzYzIBAlmF3S4vOrHu
8bAyfw4BLsIzLKZWJ64adZf9UUyARhlrVI18F5vCf2IfqgkwVHfN327QgUnx
D1cxGwgYdLpLY1DacXR6/S9qFwFpXiVhj/9qR7VTad+K2QQcLYiN0HzQhr6Z
SiOpwwQQGYpnfje1YtPvSlbEBAH3tN0/+I20oMjZ42135wjYKsdvlV/ZjDbF
9FpXnl5YWTi/cEirEf8H5ZLiWg==
"]]},
Annotation[#, "Charting`Private`Tag$5181#1"]& ],
TagBox[
{RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], Opacity[
1.], LineBox[CompressedData["
1:eJwV13c01f//AHCpRCIRpWGkElF2KJ6VpIx2IjJSyiyKyv4oUqJCfW3ZRdl7
Pa/r4nKvcYdxrzSIhhGZSfq9f3+9z+Oc9/uc13nO11v28o0zV3l5eHjUlvHw
/P+zacQq5C5PF0o7+ifrqjTC5Xvpb3xm2fh20kRjS2UzLPGNtPmNstFka0+l
yRAV4mPUfgUMsHFzPqkvq6UFtLb5SPzXy8abbbtFIsJagZlfr3u/nY2RadMU
5700cNcTsA2hsNFGp9FHoZEGq2mn7j+sYmPwDp61X83oYDD8kRaexcaxjW93
mR1ug0+35CciEtn4N3zZOtXCNvBbdmP9syg27m5obeHb1A4lW5asYwLY6Ab3
OM/ft8OpHMP/Xt5mo11CBMl9XweM7ovIjHVmY9Fotm5KZAdsP7tlPNGcOG+L
9vp3kp1A+uQgmmLKxvK3tbxVBzvB2j1XK/UwG1OsstOCHTsh5pFuYOYeNq58
TDfKzuuEFSQLkTw+NsZVCKWeUGfAK7NXGgWLLPTfpZzicIYBB/q+WhT9YqGL
SX9O9U0G3J71Ti37wEJl6Uux1rkMGFJ6oYZlLHRTuSpkIskExYDM8W+5LGR1
lXbEajLBvaM0RzSVhQO5YXWnTjNh3qN7m2M4C52O3Doo+ZAJq8s3iAnZsfBL
YXhTzigTTvLv6tA6z0KNHfI9IytYEGOpHW5nzMJ344njpC0s2LposbxEg4Vq
ny0PaJmwQOVQ/NRFARamlNx9n5POAq+onPwHS0w0MWnxGalgQdVglUveFBP7
zK5ocdpYYBDaN7jsAxNlxdLiG2dZYE7bwn5dxMSodN8uKUM2+J5PKZm1YuIt
voUJly42kLLyb8qcZuL1s4JqzsNsWDGPSsZHmbje3d5bdY4NkfGfMpJUmJgg
t6lOemMXvPoo++LISiY2nqa8j7vQBU1OGV5ReQzcJjARzG3rAsHqEtWadAYu
rheTG+jrglNrGseGYhkY4qvqV/+9Czh5Q1d17jNQY6+qyOqV3TA6tdP80wUG
9iScTmPrdMO6gDdae5YxcORI/WBHYjeElu2iisx14rXWQ43M7G5YGM+2nBrt
RBGRmwZVRd0wYJvlV9HbibUcjUBtajcUHUonHy7oxIgn1vZ/JrvhzMqkExds
O/GWd7zQp4M90KS35ZPO+U68J0EWO2vcA7reCR5bTDpxXdy5kPyzPbDta1zM
Z61O5PkUlr3DsQemqC+5LsKdOJrc/ljocQ/IS8q28fh24MaKO2J/23pA5LH0
RKNjB+b8OvBHtrsH5he2ioWf6cAIxvca9Q89QH2/yVJcsQPF66dDpcZ74Pqr
9UMKnHb0OvY/CQmhXsiR5/97Zl87xkWzcguNekFJ66dS5hQdhZIkXRILesFx
T2BgZT8dDy3F+imU90LKzrXMtmY6Vk+a2GbW9oLohr13ZhPo6CZTvc25tRfm
5tzrjY7Qsb/J2HT3YC+QK8fP/4ihoc3OaVb2Og4sFga8XgqgoVp3prTqBg5o
5gj/EXWi4cI4r1LuFg68jt/zar8eDWsjvv/ykOdAhK/7jydDrTjxpn2j2QEO
WOiNB+zVasXdkxkRJx04EKUZwDCQacUSih994hoHaMrC2y1Wt6JDesOZB64c
0Jfa0xL0oQWH7IfDg704IPfPTYwZ2oJZjzwW+0M5ME4ay77dQ8V5d7Wz27I5
IF/pvxBWT0VVn00xkzkcsCsUMkvKpaLLvhDzkjwOsFKVfzUGUtF36GPm+jIO
VAS7Hdiwi4pmdXdND1M4cP/IWGfl3Wb8o3/4j+0HDpCPjZ7UcGhGuTWjk3af
OcBjNtKeZ9aM6UlWbPMvHPA7/52evq0ZS7YKRMv+4IC34xA1gt6EKmlnXDVm
OFDi/MWIv7wJJXp4932Y48Av98Gm4NQmPPrvrrn/Agdu3PlMuePdhCKPXRsT
/nHgelg/yV6mCacfvym/KMCFrCfv4f3qJjweTvEvEOTCl2d9deYzjShY+nD+
jxAX7OM4NSatjWjz3FfntigXLuZ0VWjebsS/f67KjGzigim9rUCASsHSDfGc
PwpceNxJ3/ugiIKDZTVLhbu5QGXT8pYSKeipRA2yUeaCYX/L2ykPCv44wq2P
V+GC/njj6/7NFBzdK6OfqsUFv18UeQs+CjaY3jktqM2FqtmGLOZEA8YZvRN0
1eHCvn/1Gc2NDSj+78v/1h3ggopI3avCGw2oKUrlGT/IBff1tdJKFxuQfXGZ
nvBhLrzdWJOcdaQBfZx2X9lhwAUF2arEBMkGfM+uitY25ML1HZWbJVY0oJiZ
zFuNo0Q8FCrin42TkXP4Toe8ERe2qZXFhjSQ0Ym00XzqGBfstEo38OSR8dEw
iUo7zoUU3ZKXPrFk1ASJqwnGXNhiUBRzw42M2ze/05YxJeJlVCj24wIZXxgW
BnQRjjUpiLpymIyZ9t8Eg824IH4u75nlBjIe7M/eWH+CC0U7w6sMl5NRaFRO
wvgkF079vvZF9Wc9egjnxdAIj9GOCEv11aOU51Cc4Ski3smy2qub69FzpfXu
CsLyHkv2s0X1SIt/eULuNBcoBn3hA8nE94JCAg8J20tUlLY/rschS0W7QcJL
32I+VnnXI51laql9hgsJ1R4C2ZfrUU+96XcIYe3IE+rRJ+rxNfYdpRPustt9
KVC3HjXcmwwFz3LBU53/ocvOeqxaVvznMOG1fEMFF0TrUUWR7HaL8Nveeq7B
EgmL0/lTEwkfz01eofKDhBLnn8bXER72992zpZuEK9zP2nII3z9lYcFPJmHl
g6PTo4Rl5DSDp9+RcH/xCZvfhGtn1r39FEfCIyVH05YIX6SOd9FDSJj37y9l
kfBcPO1fhQcJHwxr06YIx7i9Vsi8RMJXeQWVg4RVD4acfX6chL/fib2gE24X
vezvr0lC/6srL+cRdhnSz3aSJWFMgtSOR4T5KzYzzguRUN97bMCasLH45Mz9
X4j9MTMpCoRDqnKOW7ciJioU2P0k4kOyvZKkkYZ4W/SJQh7hhRVSk2t8EMcm
x5auEtbI6TkydBrRy0hlSIKw+8nnsbUKiE4P0z+RiHy8njYefbEMsUEtcsaB
8EDcioPunDoMGlfaxUPY/Mudr1KP6lDI6dicHJHvZ49U98/Z1WHRlFBiLlEf
rXtGIju063DsjOQNJcJwz0Yr6Fstkn8oFkgS9eUjtfGxBakWgyP37Awl6q+k
gdGvEluLlqtucUeI+twlbBjy2agW71adE0o34YJI+m724ewatMpIW1VP1L/x
sSH5zYE1GPVyfd8c0R8PxpJ9p8xrkDQoJLiL8Pw+UbkMvhq8/lL86W2ivz7T
5z1WOlajmp10eN0hoj88ixr79apRcFVeFonoV/MNrpJl4tXYrz+iWwPE+e0/
khwbq/Ct+wg3RY/oh9nGtdTtVeiH6eclif6PzEkj6/+uxNU+Io/a9hH5sgn0
Km2rxB0mP/b4EPNje5N2X5p3JUrwdn2uUefCyxe5Wf7UCgyV+3j6GzF/PI+H
WU4nVmBSwPlQOyUunPx7ZY2LRwVqSK03ZSgS+b4q5WmxqQI9XmyNT5Qn4qXx
XE/dtRzH2rUyuTLEeb+5Tbw5WI6b4/hbNkhzQT3ROF1GvBxvXn8Rb7aVC6PL
VwgI15VhTIfzuXRJLtiwvNnf1pahy6n+wLl1XDjsae2SXFyCe3bJBBku44LU
Th0p8bASvD3UsnwtMZ//cMQZ4dYlGBfgdo2xSMz3Qx1a91aWYMV9Ye7BeQ7s
XHeY99yFYnTfM1nVOs6B1QXycQJ/CnFhJKThCZcDmX1DEd8f5KGAvqZi31sO
cP0kUpLU8zCi03b/whsOCEsZFZweeIfL615YihH7zNvuNbMS3iHf1sAG9VQO
HB1y3vB4IRctEgruKr/gwLfxn68Ub75BK+mi4pV+HFDkXSh2scrAu0WbWs8Y
csAmXbFRenUGPne20Sw7ROzbI1bdrIp0TN1xtUdMnwMLoTXz+yXS8cdsXUOt
FrF/BQP11jBSUdS60YJO7O9eRf5Iy4hkNA8a1O/i58Cl01tK+hxicGO0oUdC
Uy9x/5rR8aJE41XTzh2XyL3g97IDhXdE43e+KnvJul5I7AumHRp+jsMdXh7/
lfZCCY3ZeJIRiXxOJotrM3ohf9pV7hv3AVreFVonH9QLQexd7M0Neugv3M0D
GsR9Zm2/80rl7WD27LjOmb290D0ZnCIyfwxKbtny2Cn2gm/00VYXDTvYleqg
6yLTCw0yjx4miXhCbWKs9cE1vXDtmqJj2VIwPLFaNBga6IHK1b5dZK9o6JY9
sHfjkx54pSzOyt8WAx5REsO8oT0Qeiq/M7EjBn7c0lL6GtgD514O0rwVXgLl
BdMh8VYPTMiakhX7Y6G1S9G98mIP7NTZmh91JBlc+3xMauV7oKXY9kjhowwo
d9KcGKrpBs2BZVL/WeTDHQ+HD96sLthTdWFS5U0+uIh8v3eI3gU7o95RPv3O
B50rjap8jcT9N/Wt5iuzAnhWxPjpX9YFi6cs3v8bL4CYTGOvnXFd0F+Qr4Cq
RbDqr7v4xUtdkHLThqJXXgKL99+L6H9mg9xE9W9dciV0/ORY3OtmgV+bwYPp
sUq4sNLgzhKdBd05NOE8ySo46WY779vAgsdX++RkPargoKOlqF0BCya4C2ar
ZKrB+dlzFfojFtQ16qSx/Gsg/bW2qcp+FlxMKDd21UH4p81zdlUsE46nc9pO
WJCBEvDVjaTNgImQPMdGNzJo7Biwc1VmQOz1+//23yeD9ktqytptDPimrKyq
kEcGnrvbvuwXZEBYRVDM8hUNIDh7vUyxvxOobfJWFfkNMDr5pvpiYCcYzXt9
k13VCOwvA/ej73WAoZno8pnSZvigqSLvq06HSHuByd8tzdDp7YIUSTr0evF8
XOpvBmd8XbbwjwbOyeNV/HxUCPanOuym0yDyZ6vnFnMqeN4TEk9yoEHv8/sD
BrNUEB2IoedGtoJz93R9lFYr/La5rFrGpkLxj5H8/xm3wjnF9/YWFVRYXBpI
SrRphYPd49cGE6gQKc+4lxXaClntB1yKLxPv332nWtXdCk1/wXPfz2ZY3OSY
+tmbBnmFtyXylzfDU9ueIJVyOgyKTdL0JBshdPQADNPpcFW9bp53jgIB99L+
JgzQId7qTUYWmwLu0W4+q4Tb4AtPuVPYUwqcaF5++8OVNhjdXT7/bjkFhPeq
XHsi2g7nDnYUwBAZIpfCzL65dUDyE0WyfDgJQsPHBZPvd4COH2W1zBUSBGw8
13o2rgOK/aUTRw6QwF1N2ggpHSCiW+Ay04cgX3xRhibUCXFd73i7E+pg1dcJ
ybQrnTD/wvo/HYEaaD6xVeiUKAME27/uLnleBuv3fJbw3MIA0wqeZqMDZWAn
lCkTs5P4D15oeH9tuBQWaEoavboMaEv1WWGkUwpKxw9Y2TkwYI6eUHKaUwzP
DKze3CxhQMnT8LApnkJ4v026OKqOAZk2ssdupxeAAu9gTQmVAb/jfUfsDQuA
THLunH/PgEeBvGlTwvkwo+czF7SSCU6mLN80/rdwUTvO8Lk5EzqTpcWmqjMh
a8Olk8V2TKC9TlXvHsmAqVkZyy5nJlxyKtVQ2pQBT0pfu0oGMmHfzwpzYc80
QLWK6NRsJpCuSP0z2ZwM25V7PhfOM2Esd4W4Q85TuLkmYYTFy4KooaAbUlPh
UDtiMzOzhgXHJ2tH8tTCwDxnmF9XlgWiwjWi+18EQZj87F7ycRZs5wi9bJ7a
D2Oy4gHMeBaU+/GHN115htq+qR71GYT9A/J8oqMwmK18tSCPBY+0klczWmJQ
4uFR00gyC6rvvbWvOhGH+qN3Nhn/YMG02vEJ9cupGGa4QlhnmgVCpqWnfJPS
kJn8bNmuJRbk7mMYnPuQjo6n33xbKcoGjW7/O3yeWfi0jFtWr8MG7/3RdpGj
uchdey2nwIANlrK5zlds3+F2p6mkFDM2jDyd/hHRl4cVm9eE+Nuz4T/WrVth
vwuQ1yv2rqsLG1IuF+YrUgvRtH27q5UXGxTnCpWdkorwU5DeWZ1HbMhUOnS+
zKYEFbktR3dFs+Gj+KGTWaaleFvdXHdDEhvEHNb3CBuWYd2TAWW+bDakVz88
HHa8HPmH3WVnCthwIqkiUdq6As/An/Vfqtjg3PnwaXhgJSbGPuRnUdjwLNaW
86C8CocnxRbr29lgnvaTfEWgBlVNXv0s6GXDpctymtb+teiXoTSYMsCGB2kz
CoqyiE1/K7ojR9mgc6JYv/YhCUUuGLb6z7LBIGhj2P6merxYwKh15emCv6u3
Z5w1p+D/AQLiD1I=
"]]},
Annotation[#, "Charting`Private`Tag$5181#1"]& ]}, {}},
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->{{-1, 1}, {-0.4113765441906263, 0.41137662741931913`}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.05],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{3.874806482951189*^9, 3.874808623108183*^9,
3.8748092196232777`*^9, 3.874849221186489*^9, 3.976945337963664*^9,
3.976946680995391*^9, 3.9769475007288857`*^9},
CellLabel->"Out[3]=",ExpressionUUID->"bbf8b6f3-790f-4919-8056-f9384a4be5ac"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"p", "=",
RowBox[{
RowBox[{"TransformedField", "[",
RowBox[{
RowBox[{"\"\<Cartesian\>\"", "\[Rule]", "\"\<Polar\>\""}], ",",
RowBox[{"eq", "[",
RowBox[{"[", "1", "]"}], "]"}], ",", " ",
RowBox[{
RowBox[{"{",
RowBox[{"x", ",", "y"}], "}"}], "\[Rule]",
RowBox[{"{",
RowBox[{"r", ",", "\[Theta]"}], "}"}]}]}], "]"}], "//",
"FullSimplify"}]}]], "Input",
CellChangeTimes->{{3.8748076439096813`*^9, 3.874807707176753*^9}},
CellLabel->"In[4]:=",ExpressionUUID->"e93da1ab-6cc7-423a-b31c-eef2b211ca7f"],
Cell[BoxData[
RowBox[{
RowBox[{"-", "1"}], "+",
RowBox[{"9", " ",
SuperscriptBox["r", "2"]}], "+",
RowBox[{"2", " ",
SuperscriptBox["r", "2"], " ",
RowBox[{"(",
RowBox[{
RowBox[{"Cos", "[",
RowBox[{"2", " ", "\[Theta]"}], "]"}], "+",
RowBox[{
SqrtBox["3"], " ",
RowBox[{"Sin", "[",
RowBox[{"2", " ", "\[Theta]"}], "]"}]}]}], ")"}]}]}]], "Output",
CellChangeTimes->{
3.874807645215337*^9, {3.874807687183955*^9, 3.874807714146865*^9},
3.874808624379738*^9, 3.874809221121003*^9, 3.874849222825783*^9,
3.976945339467812*^9, 3.976946682280624*^9, 3.976947503205727*^9},
CellLabel->"Out[4]=",ExpressionUUID->"fb2a88b0-0e56-40ce-b234-c567e7abf710"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"ps", "=",
RowBox[{"Map", "[",
RowBox[{
RowBox[{
RowBox[{"#", "[",
RowBox[{"[", "2", "]"}], "]"}], "&"}], ",",
RowBox[{
RowBox[{"Solve", "[",
RowBox[{
RowBox[{"p", "\[Equal]", "0"}], ",", " ", "r"}], "]"}], "//",
"Flatten"}]}], "]"}]}]], "Input",
CellChangeTimes->{{3.874807718073971*^9, 3.8748077539936743`*^9}},
CellLabel->"In[5]:=",ExpressionUUID->"3fd0327d-7eb8-40e4-abf1-4f0bc152b75c"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"-",
FractionBox["1",
SqrtBox[
RowBox[{"9", "+",
RowBox[{"2", " ",
RowBox[{"Cos", "[",
RowBox[{"2", " ", "\[Theta]"}], "]"}]}], "+",
RowBox[{"2", " ",
SqrtBox["3"], " ",
RowBox[{"Sin", "[",
RowBox[{"2", " ", "\[Theta]"}], "]"}]}]}]]]}], ",",
FractionBox["1",
SqrtBox[
RowBox[{"9", "+",
RowBox[{"2", " ",
RowBox[{"Cos", "[",
RowBox[{"2", " ", "\[Theta]"}], "]"}]}], "+",
RowBox[{"2", " ",
SqrtBox["3"], " ",
RowBox[{"Sin", "[",
RowBox[{"2", " ", "\[Theta]"}], "]"}]}]}]]]}], "}"}]], "Output",
CellChangeTimes->{3.874807723753869*^9, 3.8748077543791*^9,
3.8748086252075033`*^9, 3.87480922375686*^9, 3.976945340536889*^9,
3.976947504029359*^9},
CellLabel->"Out[5]=",ExpressionUUID->"8fedd494-3885-4f5c-81f8-d9048cc3204f"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"PolarPlot", "[",
RowBox[{"ps", ",", " ",
RowBox[{"{",
RowBox[{"\[Theta]", ",", " ", "0", ",", " ", "\[Pi]"}], "}"}]}],
"]"}]], "Input",
CellChangeTimes->{{3.8748077627936153`*^9, 3.874807826304165*^9}},
CellLabel->"In[6]:=",ExpressionUUID->"3da898aa-730a-49c2-8505-bc1036bac955"],
Cell[BoxData[
GraphicsBox[{{{{}, {}}, {}, {{}, {},
TagBox[
{RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6],
LineBox[CompressedData["
1:eJw91XtMl1Ucx/HfclwGtkYGudKAxtC5dVHMK/alkCBAIcxAMEEuJWqK4oWc
ggYCIvtJM0ItYILBVgQEBObtUyLBuAm/85ibogsxQo0oQ5QA7TA++cez197b
s3OenfOcHdfozSFxT5hMpo36Gdd9neOSoQAFKjvnegXv8Feg4n2tceyenwKV
aoezMwd9FKh4xl+tvyMKVIbGBhY1eShQ8WnxTZn/vAKVUxlrZgz3WUDFP3B9
1es7LaCy+/bpowdsLKBSmem3yzW3E1SmeT8X/+2UTlAJasicfjG3A1Q8XJp+
r7HpAJWTG9pXz/2jHVROWW/p6upoA5VvZnWmfn+oFVRi+/NnmVe2gEq/u5/v
maebQSXnJcy2KWoCFbPX+Zln5jWCytbjReE7cxtAZfvSjfDprgeVmoje7Y9e
OQ8qN6MHCroX/ggqAzfK53yccQZUFmYVHTs8VAcqP4U5Td7tWg0qLj+3x335
aRmo3E5M87Z+6zionOhp9G+4lAEqbyyblr1462ORsT39bkVSsVDsGbX6JXxD
pVA4n7PLfWawTihaFnXesMuEULhdyjgWFlwvFIEjN31rnRqFot7trsW8rVko
XosLcO5LbxOKC65Hp/p7dwjFqdT7eyN3dQrF5ej+wOmfWYQiNerW4ZQPlFAs
OF2Y8L6VIRTzg2ceHNxrCEXSiZhJFS2GUETdemF93pAhFLO3dCUO66Z4tjIy
edV9Qyjs9xd8MeWBIRR+S192SRg2hGJb2tSjLSOGUKT3IdbpkSEUI31RoYt1
U9hGtCSt0U1xxDFnWaFuilTLD3H2pktCsc4j+56rboqzOZ6Yp5ui5rdqqwjd
FIPnPfebdVO82+Vtl6ebovzNp4x83RRjCQ96inRTuLuFLS7RTbEvPP1qqW6K
5SuuNIw3xULz1yNf6aa41hyYPD4ehY+tU9D4fBQpaeGbP9dNUR5fVJw6/j0T
IqDG2eSvm+LX1dbnFuimyLpSPNldN8W+HUuiRvX6UXjErH2Yp5viRa9XTdkP
DaF4Mr/fb5PeL4p/rwRtWvG3IRT9Q6rziP5fKAY/MZ9szzOEoqtnTm9AmCEU
rcV20wtsDKGIMxkm52olFBeT3qs//J0Sipxh9ZdNpRKK5LJQ859lSiiiL2Rd
rylRQuFb7j3L/ZgSCtvWeLuCvUooquJj+0wpSiiS8mcsidmjhOJR7dpA911K
KFxrTfYliUoo0tTykGx9niiiHT660Bur558QB6IyqrxilFC8ndM/6Z9IJRTx
yY5p76xSQnGwsOx6aagSioSQirqxlUoo2uxm9JaGKKG4c3GO52iwEorWkssO
wUFKKNq6Iz8sXqbfn/DxPfX/vfUfsKuZJw==
"]]},
Annotation[#, "Charting`Private`Tag$8706#1"]& ],
TagBox[
{RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6],
LineBox[CompressedData["
1:eJxdz31M1AUcx3EgikqZiCkuHkSXYj5ASoSm8glBjQLuPDuf0qU8DELBE9St
JgRMDXSp4FmaW5iACJSmU2SCfkeSJ3IQ3MHB/dBLcSqOaZOpEU/2a33+6rt9
9tr7z+/kuC26RBcnJ6dkdf/6/1vw9tdR2a8GC8Wk+3vCX3GLEIodaRr9ONcV
QqH9sWD03JE4oTDcsl3O6UsXisZZrZEpXblC0ZWyXrPhdKFQ3H5nlXewV7FQ
WN6KjcadSqHQea60a/zOCcXV0eWn5u+pEor7ZR9knMm7JBTnpppNYxuuCMWa
zlDfRkOdUIRY+3NO6a4KxV8llaumr/pNKMqSHniFhJiEwtl46FikW4NQPH1U
kr/8xQ2huNRvXenbbhaKdTAFlFmahcI4YebFrvgWobjisfDAnzGtQmE3bpjv
GW0Rirb4lvkj86xCcX6g997qeW1C4X5zZsvzoHahqNLWvDs01iYUd5eNXRc6
YhOK5F+dHmiaO4QCO83HivI7haKgsK5mzgK7UHQXKFdn1dmFwuH8wjhjmiIU
mzK8j8cEKEIx8VFLeOp0RSgWfTMQdmKGIhQ9B72V7kBFKJpquxJq31OEIur4
1jyJUIRirTn1qCFSEYqP81P7fJcoQpE79Wl66lJFKIq32jb98aEiFC+HVQe5
xyhCEaed9PgrnSIUE5ZluZ1Tm8Ivc3XpbbUpgoqzk4JXKEKxVzeq/bTaFN8l
ldypV5sivaInoUNtiqozldefqE2xcVNj6d9qU4zRdvYNq00Rsjng9IjaFM4V
Y2yDalNcsz5MfqY2RfQRj8RetSnWLKouNKtN8dzx+rfxalO8dmtKZ/dyRSj8
0jNb9y9W//1PuKZEXDzfaxeKF3OWLO7Q24Wit8TZ0ZrbKRQ7t5/V3FrbIRRF
sWZPFw+bULiH3c28/EWbUAR1xT50c1iEYtBkHDSEtgrF5Ia+oI3XmoXihOdP
X678tFEoXH0vmlwumIQi4XrF2o3z6oXijai8zJpxdUJRX+AoP/l9jVBkVX7u
2CcXhMLL5J10Er8IxfamwIzE8aVC8cP+nHD3hsNCUe2/I+vZlsNC8XPa+16P
JxwWCvvdlgiPRKNQZH3k7/ByOSQUg8FF1YE4KBSTApMNRXt3C0XGUM9w2s1d
QrHTZ3/V+sBdQjHU27/jiCVXKCx+tdYcn2yhKJ9o7Z6zebtQaGOPHqgYt00o
et+MG/msJl0opvgn12aPMgjF2dk+92avTxSKZd0Fv987EycU7nr/svEvbRCK
2U+iiqv1eqGoH79t4MawRiiy9Ka5zrolQvGJa93SadoQoSgbKm7yitpdSxf+
A6mzWA0=
"]]},
Annotation[#,
"Charting`Private`Tag$8706#1"]& ]}, {}, {{{}, {}, {}, {}}, {}}}, {}},
Axes->True,
AxesOrigin->{0, 0},
CoordinatesToolOptions:>{"DisplayFunction" -> ({
Sqrt[Part[#, 1]^2 + Part[#, 2]^2],
Mod[
ArcTan[
Part[#, 1],
Part[#, 2]], 2 Pi]}& ), "CopiedValueFunction" -> ({
Sqrt[Part[#, 1]^2 + Part[#, 2]^2],
Mod[
ArcTan[
Part[#, 1],
Part[#, 2]], 2 Pi]}& )},
DisplayFunction->Identity,
FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
ImagePadding->All,
Method->{
"AxisPadding" -> Scaled[0.02], "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], "DefaultPlotStyle" -> {
Directive[
RGBColor[0.368417, 0.506779, 0.709798],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.880722, 0.611041, 0.142051],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.560181, 0.691569, 0.194885],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.922526, 0.385626, 0.209179],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.528488, 0.470624, 0.701351],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.772079, 0.431554, 0.102387],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.363898, 0.618501, 0.782349],
AbsoluteThickness[1.6]],
Directive[
RGBColor[1, 0.75, 0],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.647624, 0.37816, 0.614037],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.571589, 0.586483, 0.],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.915, 0.3325, 0.2125],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.40082222609352647`, 0.5220066643438841, 0.85],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.9728288904374106, 0.621644452187053, 0.07336199581899142],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.736782672705901, 0.358, 0.5030266573755369],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.28026441037696703`, 0.715, 0.4292089322474965],
AbsoluteThickness[1.6]]}, "DomainPadding" -> Scaled[0.02],
"RangePadding" -> Scaled[0.05]},
PlotRange->{{Automatic, Automatic}, {Automatic, Automatic}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.02],
Scaled[0.02]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{{3.874807778789596*^9, 3.874807826725697*^9},
3.8748086261123962`*^9, 3.8748092270356083`*^9, 3.874856585260891*^9,
3.976945341685665*^9, 3.976946731037601*^9, 3.9769475080732803`*^9},
CellLabel->"Out[6]=",ExpressionUUID->"2932cff6-7698-447d-9c18-0a9c0223139c"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"PolarPlot", "[",
RowBox[{
RowBox[{"ps", "[",
RowBox[{"[", "2", "]"}], "]"}], ",", " ",
RowBox[{"{",
RowBox[{"\[Theta]", ",", " ", "0", ",", " ",
RowBox[{"2", "\[Pi]"}]}], "}"}]}], "]"}]], "Input",
CellChangeTimes->{{3.8748092368038*^9, 3.8748092433749237`*^9}},
CellLabel->"In[7]:=",ExpressionUUID->"550ce462-3157-48b5-a1c9-c1ef4441c8f3"],
Cell[BoxData[
GraphicsBox[{{{{}, {}}, {}, {{}, {},
TagBox[
{RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6],
LineBox[CompressedData["
1:eJwt2nc4le8bAHCyKiWpJCIrIyFUFLntrCLSLiuVSqUhFRkRUkLKKpnJ3nvc
OI49z3lxzqkoJFJIVlH93u91/fzB9bne847nucfzvNch4XDVymkZGxsblfz1
31+ZCxv2zZnRge3/P/ofGv/MGtNhy8hDXW4efShYWyU3Y0gHtysWNus4rUHL
+R1lHOhgmRC2SvWvA8z9mdzbpEaHax96q3ymr4Nh634vdRE6tG7vNrj4zhfK
A87I/hqlwbuLpy3sssNhccvylWbeNDAI+Hup/1giBD28OVMf3A2XdMKk46/n
g9GalviR8C6wrtCtynhdDi8+n2wY2NAJzl7VKzJka8E+ylH+5Zs2OOtq62Sm
RoW1fdGrz8u1QLvr+bDK0SYwe3JcmyHXCFvE/ulEjLRBeL3ZmTt9FNh5rsDD
wrkLGi6+iXI5XAN5JzffXQc0+HjglbLaeAkc4fkUWL6XgAx73oQVTSVQ+1jp
xTo9ApYpO5q2JZcA9+GVM5dMCJCPvcsveKYE9rDveSZ8jICrTIE/ZZ3FMJdV
8evKLQIUNvmwa+QUwbiJqeCPHAIUN40HFpwpgFMP3Er2SfdAU+ndbK7fmeDA
f/bz2PJeyOE89i0pPRPCm0rPdq3pheS1G+8KnsyE07sKuIsEe2F15Y597hUZ
UNbNF+0u3Qvm8Wo3t99Nh5wDN3VHoRdqLv/6aDyZCp3fi2od3XqBMX1395eS
RFjOtuF69Kde0LOwfDzlkAhH79zq3D3aCz+edJ5uW50IA+JR57onegEmikwX
HRLgh/p38b+LvbBpkbpvF288KN+sYe4U7AO7F7T90m9i4IP+Qv0a0z6gH5vT
38z+FBT1Pa01c0hvurbJ+fYT8A7ytZkt6gOKeljm+2/BwNNMLUiv7IOmoQPK
Gj2BsCScasHT0gdCHLTtdbG+EPBlOOjhcB8selYffsd2Aeq/B38MFWLAHZev
nyUZD1DPb4F+xYMBg/mfs8o8/FGmdK7I3ZcBdkFpOlKSAbjyXJ6wRyAD1tEL
pV2uPsIhPua2qxEMKJL9wvwSFIrRLzVecGYywHxd34483UhkzLh2fmAyQHTH
/gX7+4loqjY/0qHMhNaX36odRhJR8uKhrNCdTCAcuhcfH0zCl5VrPpnuYcJ1
6SKR01uSceZJr9VbPSZMa5+OHalLQRdBhwwuGyZENv+av7Y2DV98vHBWy50J
42WWe3/cS8MtUftTPD2YMCKdFhM5koZP9ikeKfZmQr5ZHdpUpqMQx+bG1UFM
MEpXyA11zsSWmtLN1tFMMHjwqaihJQebuSVV5UuZICNC4aOK5CK7WwyFWcGE
I6F+OepWuSjPubfKG8n75V4Vs8BcZFpNpuQ1MCGQP/isXkweurSvuP2KYMLF
opFdst15eFkeb/zrY4LEzyeFIsvz8RJT2OroOyZo297lV3PLx3Y29Wujn5jw
YRfNX9C6AOPCHh8LmWBC3CH+5TcfFSClUyjo1Q8mdCT96xuoLcAju+bWJc4w
QUjvaeqnHYW4onNo5tFvJjg22fUFrilC6ZV2K1hcLKCd32i7t7MYJydc8s4I
s+Dt9+0M0VUlKKVMf/d9Mwv2vx914TMpQTn/1wZXt7CgfVUAcwOlBLn2aXgY
S7MgIejDj6KSUtz+VUfFSpEF2wyDn/+eKcVC53tL15RZ0PZdV8hUtQy5mhmn
fFRY8CB2JGR5VhlenLIJctvFAtCv+fI7sRx3lNqZ5Wqx4MgVDeWQj+V47O2R
O8e1WRDDfeybvFgF1h0efPETWHA5UerhnegKPHPAweK3HguStW7VbO+rQF+7
T7XnDFiwU79C4cv6SrxgtZpOMWTBtLWhybXQSkxt8fOxMGaB2Zc/1todlShM
iMx7mrDgvsgLHYFVVShsu1Dx0pQF7/bllY2bVCF/r29lhhkLqJ8OObcEVCGF
Ap/TzVnQu+HdQha1Cu8md0nGHGCBTtLatS84qpGrdMb+zkHy+Y6KnvDVrcZD
RNBTYwsWnD726N51r2qs0HgSxWNJzs8D083nqqrxs/3mi8WklXuKB08tVqN4
/60fhw+R4x0IvXlkD6Kl5G+RIdI3DDlvW91GBNvtfXZWLDhp8DbEsog8frSA
v4N01n4fW8ufiE92nijeZs0CLn2fhrdCNfjjzNpiN9JjBM+HNPUanLWa/JNH
ekDXyyDjSA3qu6t4vCctf7atPfNWDQYsaG2ZI901HmKSHUEef/G65y/p3v2S
gTkFNVj4gD/kv+M5EXHHc2k1+HyPq/oH0l+fMp/k/qjBjgMdNfmk+cLOz+Tw
1+J3E7tNt0l3Brw8la1ci04+j3f893zm/x4/zzhYi3e/fJ5tJZ8/QWbNlVSX
WuwbVjx+irRPf0VhwuNazNxuY/qOHP+/1kn1mIxaPBrmUWZMOuBaGzW0pRZv
MmZfJJPz12O0Z7v/WC3ahxjQJ8j5Hh8e1ri9vA5nGWK2sqT1pYzKzsnWYZzJ
FYWDZHyoxq/uWhvVoW+2pbwjGT+a5qzJPqc6/DooZHqOjK9iqfSMlF8dRsZZ
eB0l498Q1GnIk1SHKgaF5bvJ/DjVdp73S20d+mTvm+Ag8+eTsvJaysc6fLsk
vAr3k/lXfn71VTEKjvZmdP0l89F35PXFbC0KWq2Tuf1InwWyPkEOYycoKM56
78FJ5m+clerJQ5EUPLRF241C5rvsjdl5zyIKGnY8t1q+jwWCjjfPp9ApGNXs
b6upyYKPiiqzQ2vq8YzBioFL6iyoZXfK/q1YjzlcW7ScyXpi6Olf4zWvx/df
aXGH1VjgwZ6mJxxYj7b6JeazSizY9UbyyZd/9SimELbv9FYWlIj7zTeJUlFU
ma+sWZIFzro3uhM1qain9G5SSpzM3xPfuOA2Fb2611xNIPtDWk1n67oJKo6h
2UfkY4Fm3JvVBbwN+Khg2C6alwV+9zMlTOQb0MGJ19FuOQv60uanjp1twLAY
Y6KYnQUhhmw+7KwGDLETIvx+MkHYYPrtofkGFN7b0tQ8yYSaxZuckesbcbp5
THBpnAlrxxUC2CwaUY6Tnis3zISq/cXrTSiN6M9l659AZ4Iz/0nVfR8b0Uf9
0XfLTiYcr+8bkf7TiNdln0Z/byH78XnY37G7CTm+8ZRO1DKBVrHf4kFaE5oW
DSyuySH74XzqY76GJhTYNuYin84Ew4X9zcGDTZj38J7EjhQmvOApdjsp0ozb
no0e4oplgmifrazp42bMj+Fc2enPhOR/V8L3XGrBjCCbK4NHmTCjsl9z2q8F
p4TLL++1YkJsbfmql3EteMSa44CvORNCLgU8retuwW0aCo2fdJjgk+U5dkS9
FSPVnPZ+lGNCW93YHo9/reixmyv14ywDOox4m3cJtWGcS/PzPZMMiNOW3NG/
ow3rYo52+44yoNqp4tdv+zaca8s6NMJiQFKHap4dpQ29Gjndp6sZICIXqo2+
7fh2RFBMzo8Bm5bunTkS047fXyaee0Su38vGs/RYee1YSZxK+HCTAcek89+n
fmxHDa3DYkedGDAMCSuO7uvAuO66qEhDBqwQV/owPt2BHV/tDp7gYEDQh4eP
hWW7cH5jVs/Bm31QQ9u+/a5qF7rR5lhiF/vAStTrWdu+LuTR2nphwLYPVHUn
/E0Od+GZ54/cVM36oO7c+YRE7y6U3aY9KynRBzFzbnWNvV04a+v9Vq+lFyou
Jdgf/NSFwQ/5nWSwF/Jb+4Uax7sw1LJux1IBuZ+69TwzhK0b5zR0T99/1Qtz
L4743pLvRl2R7NGBq72QSTcJjHXvxsRVdrLKAr3wUENQUMe3G29+cTLQ4e4F
kTK/A73B3egkaiVm+LsHzDvvv+mLI71JyH3bYA+ILDMKVqzvRlbnY9q5vB7Y
UpitI7WKhjd47/VVH+iBH8UhLIUwGt7LNdEqvE+AYGXE13tRNPyd6at+5DoB
p04PnMHXNCTes8dNORHgIdjtIZ5NQ3/UWLnmAAEJyTk0sxYa/sCV32aEyf1r
XkSgCRsdP161Fr9TRIeLNrGK3dx0ZG5hL8h8S4crnRcvmq+mI4f0jCQrlg78
G/ZWiQrTUSIob1bSlw5Zw3Sz7Wp0/Hcv+eOVg3TwCbkd/ceBjmzn/jlzjNCg
xYBbsPYCHa9cMzo1wKDBs4TD3u5X6Hj+fXVocSsNNu3sVKfcoWO6O8fGo3k0
UF4TM970lI5WvkE/HTxoYHQwPyC2jI5jXAN3SteS19v9aWNCNR33tjY8tuKk
wZO/gstjKXTklyNUx+a6gZl7Re5yOx29Lx0PWvW+G7IenzD3+0hHrS+aYYpv
uoGqb+L1jJPAuUce2l/3dENxIy1TewWBHF95tXYqdAOvllzP+9UEShQuJN3b
3A38uaUG04IE1gbIrOX41wW/a5+O+coRaJSXNj5H6QKRhG0sQxMCFZTKXR8e
6IIywYOWMgcI9DnY2pqn3QVxhT0PflsSeONixNw75S7Y1vvN5v4xAkNeywso
CnTBojrX6tlzBEa8/GsYlNUJI4ekJNN9CDxrsepIEXbAZg1LjzY/AuMF7wv+
SuiAKSOTs0MB5PUTlKe0/DqgyHFEYPIJgSrZV32oxh0gXjVv5xVNoGv8Pp7+
7nZY+Ob9ev9LAh+6f/GSKWyHD+VfezlfE3iJoaN47UU70E33nDuTTKB9t2UW
z8l2SP+m6/00m8ArGrpDpsNtICO6+x9bHoE/Xbckxze0QZz47OD5AgLPhT9/
vfC2DZT8B/gESwlczBiry3FpgzUvM6/P1hCY/mnTWu35VlBdZqLR1kXgCwu7
CzMrW8n3qFSOcBqBdtzsG7gmWoDqd8jpIEHgzJRW0qbuFtBVkhLL6CNwW//P
hoORLbBH8F/D7AcC5dil28alWiC8aMMb1wECy4Wuy4nxtIB4tor38EcCjymV
Stl8bQaD71eP5gwReGbG1bI7txncz9lPsEYJbAoYHOLUbgb6eZll/F8JPPxR
NfaMRDN4j6rd1BoncP+nzpdVnM0g7mSQfO87gS3x5zuDW5tgTL7ubsgEgQ3B
07KL2U2gkOr5MWqSwPUwvulaeBMESps9ivhBoIvOScFLx5vgxfcXI/7TBCpe
MDT6odUEuTnPll39SeA048EJT/Em0HMoTJObJTC3ZoV76pdGEF2zlUNzgcCV
Qiauhm6N8OhlY8Y4ab7TxTlLxxth8qv1jYhfBHqbL50s39cIFC7+kI7fBE7x
i2ibcDXCYVr0gu0igWt/n3omPNYAQ5K1uWOkeUJNn0+3NQBFV0FjZIlANe9j
74ojGmClwS2x438I9Ez4pZXk3gBCixe960kfYxk/fHGqAfbo3nKW+0ug1XlX
aqhOA8TzTQ/7kRZOclweLt0ALnY7/zJJJ97qORWzvAEcxP1qZP8RaErR60r7
RoXjq8XUr5A+8v771ZouKjglq7lkkbb7NmDcX0iF+JXjLp9JzzuPneeIpkLF
/BkDQbYetHyYyFT2pIJYQtwikP798lK6oz0VWlOy4hxIFyW5jMQZUuHD7wca
90lLxsREf5KnQj5jNRFOOvpXVqcCHxVK3ip4xJP+J60Y4TldD60/8zRTSTsw
zk/39taD/SFn0f9cHOU+q15RDwJH+bb/9/njUoZv41/Xw0Lgoev/XS/uw/0N
a/3q4Z5p329P0uIGSfuDLtQDs1O967/nWXViyXjFgXoQ4/+3pENaSShBKkyl
HrpKMoOFSH9+d+uTuGA9rN9c+HCMHO/nnu2PS39TICjmcX0J6WDubHmlTxQ4
cfHpCW/SFfcbKtMbKbD2ssghfdKJVz6ZK2VT4BnXw1w20qrGhp9KIygwKCzi
U0rO/4tV631M71HAjEe09SLpR9FmakP2FOBk0/1WTcbTKWHPV3llCmzuL6Hb
kaa+27PUt4EC1xQMbf6Q+TA+t1vzyVIdpJcNGsmQ3qWkdXx1Sx1YvvrdkE/m
j8tXdwNGbh1Qtmuw7yUtbrv5VnpkHXC3NFE0yPxzfhVSddapDsore07mkvmZ
afvwo7l5HdwQ82iXJG1or26rpVYHZc9fnVqYJ+8/kOenwl4Hdksm7OfJ/A/Y
Fr5481Ut0C577iueIZ/fqDYwyq8WNpw9FcpGek7dOopyqRYUHo5HPiTrqa4w
22LH3lrY1dyxdIysvys0rdCzfTVw9rxUlQ9Zn52GFon11TXQaVfbkULWr/Ot
7zbb39TA3xqXZR/J+uYTMcoRuFUDFYSrqyjZDyqta30119VAZCb4SnwhsKf7
1UzeLJL7M+/NEiPk/ddoHl2qRwg2K1u+bpjsT+VmweiIYD9l61lH9h/RoWkV
tsRq+FF9OTKpn8DBxmLdEtdqUNdNveZN9itX46ejnrrVcDWt5a7iOwKHFnNW
7P1UBRPNRYk2vQTGzX+ddpGogv5F4eahVrK/JG+LakiugK+Gy5edayFQUDPy
qKdrBQyZe+YNNRF4fDH3uZF2BVQvGYR3UwnkrJcN52eUwx3Kxe3uSK5fQrxx
vqvKYQXluuQ82a8X4ru+9LuXwlOJtl07c8l6lzicUmVUClIreNVcyP4+d83w
Xfb6UkjNPc1NTyew5PbRt025JbDpUGbrlSQC+7f88GOOFYPgSSO/b88IFIoT
zXSxLYLd7FZm42EEZnwsTLmuVAQu7RccR54SqLzz9nDQn0LYEvIlrjuYXI+e
2m+ZjC2E1azfYx4PyP6afJIzglkAWoMagwrXCUzOFxnTO5YPrTFiWWvJ9fJ0
tKb/uoZsaHm6orTDlECt+MnBlivZUJa7VyrImMACZmVp/MZsGPLskp/TJ1BJ
bkVUqnMWsMnS7fP2kvXCvTJbfE0mbBaMLMqRJXB2ylWw0DYNOq0ztYb/0fGI
0AuuX4LJEBQp/oKZRkc7P9fVPwufwf1KbkIwlY7O7jpLtI3PQHLcx8gqmY5C
xtp1n++GQ5BjTTcljo7vrnGaJemGgeOz1ZVhz+hYNRAkOpv4BDQs2lJYHnTc
QXiYzs88AHaOYl6lg3ScfLVfZmynFkQ8m1nc/Y2GnMyHG4aEwlG+wpvaMkpD
KZUh9tLicKw5qT5+6jMNY5NDvU4efoY/qK2WHv003DBVN8z3NAJfrpfISO2i
oaNw5cHtyyKxPGukNq+Qhrujzlk6D8RitbBDeew90psTuCIqElHCpfrUIicN
KQNm7yoEMvFI0Jb0e0LdOHNS6kDV6WLkOFoaw7asC/l4Xzw48bMaNUPiGxVz
2/Gudc/NmEAKyvUUK2871Irx2/LkijUb8fcE/UDT3SZ84xJ80UKtFQMkL98s
8adiFvHjpunnDjQQXOhO1K/DpwTzSIJ+Nx51fG5/ya0S33nWvR7QomNpRKeT
wJFK3LcpOVLDgJynqKnS1F2VKPrd3CDMjI76OYJCr2YqUCF3Jl/3BBmXS9wW
ea4V+DdkODnuNp18D95idPNCOdY9GC83yqejNetG/ArjUpxQ4ZrUIuOcN33X
NbmtAGPk1SZXrenBukwT4+PJBahs8uxkr2AP9sq6xszeK8AUL1/XOLEenOZ6
c3N4WwFi1pDPNqUe9Le+bRsdkI/jXm0BKgd6MPODpr61Vh7aTnwNmH3Ug2Vq
fITp1iy87fgtTpqrFy90yoie6sjElpA7GsWretGk1HVC43YmJnlVRRqu78XI
LpuWw00Z+NtgwP+0VC+qS6W/E7iQjkrCe1RcdXtR6lrG74hXqWgt9YJT3LMX
z3hLsWQ+JqDmGTGF49O92MWBWsWBCbiBVnl0469eTDaw696gkoBt0cX+tH+9
+C1h+QVDn3jcdcypFlb3YW9K4YMW7Tg86X2kaEa2D0XglpRuShTKugTp1Jzq
QxBa4A2wfoTbPMKiIil9WG42lJLxIRB5szl+DTb34aXtrDm+8wGYMRq1fltX
H8odqZwJu+2H6fkXLDLe9+H0hZVvBbd4YuI411P32T4UV2FoCrw7jAJNHd1q
Mgxcu+9ajMp7f8i9KWPS+4CBzWnHWqQcAmDBUWCtUxADh/rpE82jgeC6eE9z
IoSBE4sJ1TsWgkHsxKsHk9EM/Gj+UD9nTRgcuqZ0LTOHgfvVUxVc90dC/wh1
kYPFQMvHbHw/10SB38/rGaoDDLzhfia+lREF8zcPfT4xzMDfMdbdXpdiQO2s
FUZOMDDUxH7uW8Qr+FhYe792GROXeNR+ti0mwNv7Kx83yTNRI7o2Nt8hEaQC
fbYGKTHRyfqyVm1zIvQlcdjrqTFxYeXD/GNRSXCTd6tvnBYTTdxbTlxST4HB
NTTq4kEmFgTNS814vCXXuRtTWteZKLLnbJDVyFvoyrMK8HZjYuDKHSW9Fmng
pHylpeIuE2PtC3wspNKBl3FQcNMDJhIs9j1LrRlQWrp3hcszJq4uDLu8UjoH
bvNa21/OZ2LMSclljVE5IH21xv5GMRPdK8xe9XDnwjbvnnNXy5m4fFVbyMy9
XNg+SigY1TFxb6isgND5PAjODjh+uZuJD5L3pa1JzIPJkW71TT1MPIlrHm36
kAfVdZQb5QwmFj6+fszROh/yXbJymQNM/LLinkWzTgH0l2+bDvzOxDUnHe0O
eBTAIZ2Lrr1TTOwbS1s+WFIA8fEibzbOkOeXH96jpVwIlx9Evrv5m4l3F4wz
ssSKYMv0ZfkobhZOzO/YZMFWAmmNt79d3cxCQ6e1hjVaJVCR9fPLhBgLt14r
k9C+UwJD4tomdhIs9Ipe+f7MdAmYvFb+ukmGhbtVqoq2fy6FtkfXHUaUWEgt
S2hflCiD4j/rTOZ2sPAgQ92h50wZvMkJnZpXZeH+p82sNwzyuHL8jq7dLHyh
NdlT3lYOjX+MtWn7WCh27kZyz4oK4GgMn3EBFvKYmIYtGFVALHi1/NZhoeQJ
90c2tRWQtpX/4Rd9FsZ99JgP+VsBuzZbG+w3ZOGvlzsMOjUrwfLEbXqkEQtf
ZwuqOxZXwuzFBEFeExauN53RKP5ZCcmnOZMVTVmYWsPs5lOpghdOnEU6Zix0
NVmffflKFVDpYbsNzFn4U+7Inc6MKmjLLt2gcYCFjwZK+nePVYGVRwJsPsjC
YoeywESZatjCG54wRfqMrb6EwNlqaN83IVdkwcJZ/7eH/BOqgb9Xg+JsyUIt
0XTTQVY1nGoPHuI4RI6PPbBbdz3Cxv63XKGk015fjUs6gBDvSQjzWbEwK1fc
mScAIf3z5fX3SU/qTYxerkGoTvk6NEBaoxsJ+i8E1s5Pd9WsWWiqFPTj1aYa
MHrwr8WdtETj0lmr3TUQ8pi/IYe00LXQLTzWNRA0tniaQbrLaWSx4moN+Do0
+E2Thn8uxLXHNXDKJ3z7X9JMs7MPtqbVQLlUwv7F/84vIgaZ1BqQWW3VNEaa
Oym67clgDWy8I53QTFpc+rSw7r8aGM8sqY8hXR4VnfJTpBZqglq3niYdbx+z
P0WjFlaI6Zbw/3f9MA+mjU0t8KqLXC4ix/MoSXkn1/Va4J6+r25G+u4Otx0F
IbUgd6BwFY2cn9SskTe2GbUgUCozYEKaY9u3qyubakE+2Skx/7/55W90LBiu
hTPrt5quIs377cmJE+Q+Nm67dssRMh6vVZW2/xOtg8qaHMFwMl4Lq11yE/bW
wcRgnFQ1Gd/TCqPVukfr4LXr0c8sMv7bnggqD9yog65kc8tRMj+KbE0Yd0Lr
IC9m0/ERMn/YuzPD1mbVQdWQzwxB5tfUjWVh6iN1kKjMVvJgPwsftsT6UZdR
oOfa5SZ9Mh+PdhQNWWyhQFqEhO6MAQuT4ID28WMU2NY4pi2jR16f32Wk7yYF
rGUzc1LJfH+vctjaKowCeUIbbTeR9RBId6fvaaHA9A4163ZNFlrkuIa/GaGA
x80/fmv2kvX7IN2Fj6MeUv2TOnU1WPhgVKOzQ7MelB+XBbvtZGEIh9mL1Vn1
cDdyY86W7Sx0s/V4erC5HoJUHxR/kmdhmd01i6DP9WDOvT3imSwL1x13vz0u
SgURiVeFDZKk+6pzTodQIW+jSoOzEAsTVffIXkinQtKno+VpG1jIiGjlvdRA
BX29gYNMARae0Iu7cvwvFZ6XiujxrGah+YqHLLpLA/xyKq58xcbC6k838ieD
G+BfjvvNdX+Y+FLox3L29AZ4c2d33a1fZP++YXKYY6QBLquF8LJNM1FUx/WY
9+lG0FaNZ20bZOKVSoaQ3r1GqHdcTuXsZ6Kk3PSr+ahGCO27sbWVSfbHP7uE
9xGNUE130RYl+23bmV3vFEybYLX711e91WR/fLJY43++CeRnZDjWk/16OC2p
qNuvCaTfS9D3FTFxCyt7WhuboFv1T8GRDCaqmrcQF9SaQSK/LvTVCyZ+0Cj3
bxVqgQf9sSE+zkw8uttZW3FXCywzKTv30JGJi+7cRZ6HWkBVx59++zQTPevq
+b4GtcDP2Agf8UNMjN4o8/nrrxZIh7SBIHUmHkul+Nn3tkLJvnI9Chu5vhUs
RUX+aIVzwmIVEb8YaO3lY1O2qg02ex4zsZ5m4PHEslWo1wavBM7SI4cYaMZ+
7XJHVht0nukbBioD4zfwj/V4tQNz/pbRRn8GqldqT2lEt8PzjZM8XzwYmLC0
rNg/vx00OrhPpdxk4B/BkaF3w+3keyPfr++ODAzpoWa5GXdAYsGa2lxdBnKk
WN4z4u2EFKspXd7ffSh8sthsp1QnvJwVHG6d6sMmZYeUFZqdIKepst3jSx8e
/3DL7eLFTlBL2riUQfThntKA3exNncCSmWgyze7DtEX7uZ3QBennQ0wiz/Rh
2LOtOvamXeBcNzp88nAfmitfcPG26YIMnZ+315r2oenae3qRl7pghe6Xjxa7
+vDXzrmIrZFdcFHdw2wzbx82yLvFVIyT558v9Vle2IvLedMEQx91Q8URyeD+
xR58+cz8n+Mz0sFp4jjVg+emlWIUXnbDQ5Cxf/65B6lOCsznWd2wpPFcU6qz
B4s57sv2dnaD5dm71IaEHjQ5DPMcAjRQnKV06Bj2INcSf/taYRr4jZo9fbin
B8fO5aqtk6TBJcNrDhTFHswu4M+bUqFBcmVMrCS5v+1coDgJHqLBriPZfXfI
99aXrq6rVj+hgXTVwauv/Qn06tVcdjiCBtS3livt75DvpTG54SGxNPBlfJYW
cSHw1+jz4P40Gigs61G7c5hAWcsY0a9UGqhdHNs5IEkgG1MF8hdpEKw79eRc
NR03hf15JnGGDlxLh3r/TtFQ/ksm95mzdHi9RoEYHKShZpOAz+OLdFjfbW5Q
S5DvNcuvZ6IbHRicrl9cSmk4eyW91iuEDm3JZ2WDvGh46M3We1cq6MBe8OR1
/yoarnI3ceisocPvOSVlt7/dqNpZ3C/ZQIfJxQ/0FVPdmDuhjAlddNgQX16z
ld6Nb0v5NNk+08F3q+KAanQ3Gqj7HTPhJeB4gdfVEulunIt1uMlYQ0DO09tN
IoLdGLl0ouvEegJaKQElnjzd+Fz/ttseUQIMfNYW7P3aheqx+TW2SgSUPlj/
PDanC3XY4vJPHCTA/bvW9lV7uzBjz7yNlBUB1YmaeocVutDl+LuFARsCdE2L
jWM3d2GbpOLVXacJWJdmPGX+txM773/xtLpEAPKzRM8HduLWcVlDO38CjuS0
cN942YEnD75X1Qwk4NQOfo4x7w6sf957clUwAUJrR9RtnTqwMoXjUVgoARs7
1KrNlTrwrcSYt1wsAZTNS9bKNe2oL2pzvfgVARJpcorpye1ID1GZ2RtPQBxa
iMkEteOoKw8hm0LAgcngAgmrdpzs2LwjPJuAqYSBwG3Dbeh9YqqrJ5eArWo1
y/Ka2vB60/RR/gICjH/tOLsnqw2Pdd3Qv1hCgJ/KiI6lWxsWsX+vzEcCtAoW
V0Qvb0PFWMPY1x0ElGnets5QbEWBSGrKkS4Cfji/PHtSoBX9pBp0uWgENO78
184334Jlqr3HTHoIeHLW8Y9XTQtOfbPefuUdAft2cKwOsW5B8cHV3wbfEyCX
sl78pEYLztG4wy36CQju2K24XbQFA7qayvk/EXBb08iXNdKMj1RXp6iNEJCR
oS918W4zhj08VXnhCwGjgpLrj9s2Y8n78NUvRgmonb0fbm7QjLfTpQOYXwkY
/KbrqcvXjNX18vWT4wTYGH/m1f/ZhNqfLpz9+40AT2XfyyaMJsxa+2aec5IA
07DWz2cTm3D00+Z/S6R1y5ZtuRvQhKxv1eHfpwiYP21yKeJyE77VUwoonSYg
vLbzHmN3E+41/oBzswQoV7fJxbc2op3hu9jcOfL8wkGtc3mN6DAYEes4TwB3
R3OySmQjJv4oUChaIEBQWjqF4dCIcvvqO6x/ERAUn3O00LgR5bM/l4+TvmK1
L/y5UiMWgMMznkUCVLZEdl/43YCBnErhwaTf9a+ePPGxAd35vv1dvkTAjbz3
VKuGBtzIrB30It0/N2pvmdmADWrtRlOkRdvHO6zDG3BESXrniT/k/J+8LnTa
vQF9V/QUVJH2Euk3uXymAZczbEJc/5L55SUQ9nIb+fkI4cVa0l6nGiqr+Buw
JEFSbNU/Ar7c2L00PEdFHunDbJakXfdPHBf4QMU43ifpT0hLXrX4aEihoi1H
0lYqaepwWLRXGhU9HRzvzpGusqh/ik+pqCD2LE+CrQccLNvauNyomFsxQTci
7dt67oTVKSp+1tv21Yn0nt1yeil6VJRT+bF0n/SdTx5P/spRMf7yyPpw0q6G
aYZn1lCRXSN572vS9w+03Kifrcc4mUq3N6RDk3rEVd7X4wx/SXcq6SeL0TbJ
dfXk+DccTiItaovCYmn1KGh1gDua9Lh7o3vc03q8eqhjKvA/D6x23+pWj3uD
/ojcIF3q0CRecKoe/R6rhx4l/dmK595+/Xp0mOu02U16urMoZFC+HkXs195a
Q1pGIMbRj78ew0/vmxskx19bZMOuOE/BDZcef8olHcx17ur7DxS8P6kId0gf
NJQuCq+nYL6y8UYt0i+XqEyLDApeNl5w/UXGY0nn9PC6cAqelrU6mkd69B/j
Xb87BTP97DodSYuacDXm2lKw8L30ID9pg6jgzCAjCopdjIwsJeN9VexnuLMi
BTe92bNygcwPDibjnvZiHVZbKvQ9JR35tPKR2mAdwhaLc5Kkg9YdKFBursPQ
E6Zf1cn8811BXNgXWYd/JqMmyn8TECggz25xvw4VfvIQGqQVUkebzzvVoY10
h91WMp8lgrmW5ajV4TGVQxsjyPyvC6Dc+yBch3Nnzdv/kPVxYyFr77pldWi+
PtamjqwfuVK32IiuWswwEEzlnyHgfL3xxUmXWjQYTDMw/kmA/Ofk7XrWtbh+
76LBXbIePaud3r7SqMVvqwLD28l6VbCW/HiFoxbPDerYqn4nIEHkCnE3ugat
eDas1Sbr3+b195F1XjXIq1zgYkD2h7PHePuKztZgvdkJU90xAjZredqu31GD
mqpdD5aGyf44Pj3mHoLI3dx2v3+IAHYv6/w9JxA5Vh13rhgk4OcBncurZBCd
u4KVzn0kQESe0fKlqhqbD9yw8/6v/729jW4TVcjv03Ktg+yfxiIym39ZVWI/
26udjp0EOAvIrH23pRL9w/rO/2wnwFqhV63zWwW+iq/ay91KwIC16YqJhxXo
9lzUca6egAs11E3cleX4fmxBhKeUgCHNA5bx8mV41GvX/iPF5PoS/tU/Zb4U
qZ/EbRIKyXhdYX2vpJai/jWFnQp5BBjt0LKQdCjFgdTAcr50sr5zN8s1xpbg
m8gnr4TJ9edA8n7Z2wLFmCX8wsw8moDlzmyKg8NF+Oo0X/XdSALyKY5fTpUU
oWgpN2fHMwK+Jfq7+JwqwsIb+18efExAZRvaPUwtxKYIPulKDwIWRkOiPXQK
0C+Mc/2FkwTISI8ILnrmomPELl3N42S/vvS35LBVLt68wFXLe5QAuxznu1Uy
uch5X6spmVyPl3O8eeH+OgeXBALtS0wIaJBLGX5vkI2luyp5HdTJ9Wmp5ZH2
8wxsn6IrjvAToBhu1sFt/gbzH82Xi/OR9VL5XDVvIgUlOYLuHyP3D2d92Ld6
haVgvMTKcuQioKdQ9Yx/XzKWP5y57fKbTvaDohoNpySkc0bGmgzRIfKGBEs2
IB6d2Ai2LQX/fQ/uo7u6+TmG/qJP8eTSofjK3o0Tgs/RM/NoyEQmHZ4Pdenz
O0Wgff2j/sI3dNhjKt6/cdkzNMzW3yYTQ4citdelShCKy9ucV8Z502Gn0oVr
rx/54+1XsvscPemgvzmk+LSSH7KV2JvL3KWDwLcFtyiaL0oVs/G+uUEH/y2V
dJ/N3uhHP2j1+Bwdrm6iD6pcvoXBdgH5Oo50sNjs8Ne24jqahn7n+GlLhwPi
Fyq9ea+h8/0NfoeO04FHafNnxdNO6GqVU/LHhg5RR8RTN3DYYcdK2ZFUKzow
f5gkldrYYOubvrWWFnT4e6hRld3KEGX////DaUtJ7RtNnlf+Dzphn4Q=
"]]},
Annotation[#,
"Charting`Private`Tag$9149#1"]& ]}, {}, {{{}, {}, {}, {}}, {}}}, {}},
Axes->True,
AxesOrigin->{0, 0},
CoordinatesToolOptions:>{"DisplayFunction" -> ({
Sqrt[Part[#, 1]^2 + Part[#, 2]^2],
Mod[
ArcTan[
Part[#, 1],
Part[#, 2]], 2 Pi]}& ), "CopiedValueFunction" -> ({
Sqrt[Part[#, 1]^2 + Part[#, 2]^2],
Mod[
ArcTan[
Part[#, 1],
Part[#, 2]], 2 Pi]}& )},
DisplayFunction->Identity,
FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
ImagePadding->All,
Method->{
"AxisPadding" -> Scaled[0.02], "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], "DefaultPlotStyle" -> {
Directive[
RGBColor[0.368417, 0.506779, 0.709798],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.880722, 0.611041, 0.142051],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.560181, 0.691569, 0.194885],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.922526, 0.385626, 0.209179],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.528488, 0.470624, 0.701351],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.772079, 0.431554, 0.102387],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.363898, 0.618501, 0.782349],
AbsoluteThickness[1.6]],
Directive[
RGBColor[1, 0.75, 0],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.647624, 0.37816, 0.614037],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.571589, 0.586483, 0.],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.915, 0.3325, 0.2125],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.40082222609352647`, 0.5220066643438841, 0.85],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.9728288904374106, 0.621644452187053, 0.07336199581899142],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.736782672705901, 0.358, 0.5030266573755369],
AbsoluteThickness[1.6]],
Directive[
RGBColor[0.28026441037696703`, 0.715, 0.4292089322474965],
AbsoluteThickness[1.6]]}, "DomainPadding" -> Scaled[0.02],
"RangePadding" -> Scaled[0.05]},
PlotRange->{{Automatic, Automatic}, {Automatic, Automatic}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.02],
Scaled[0.02]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{3.874809243955184*^9, 3.874856586093244*^9,
3.9769453448886623`*^9, 3.976946731655747*^9, 3.97694750895682*^9},
CellLabel->"Out[7]=",ExpressionUUID->"81f89cf2-fdea-4cf4-bf5a-9ad47eaa8992"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"ps", "[",
RowBox[{"[", "2", "]"}], "]"}]], "Input",
CellChangeTimes->{{3.874809384872851*^9, 3.874809385762994*^9}, {
3.97694669997218*^9, 3.9769467135575123`*^9}},
CellLabel->"In[8]:=",ExpressionUUID->"a3ef92f0-fbd5-469a-947b-e93a159845b6"],
Cell[BoxData[
FractionBox["1",
SqrtBox[
RowBox[{"9", "+",