-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathHFunctionDiscretizer.nb
More file actions
9316 lines (8876 loc) · 472 KB
/
HFunctionDiscretizer.nb
File metadata and controls
9316 lines (8876 loc) · 472 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.1' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 483439, 9308]
NotebookOptionsPosition[ 472226, 9127]
NotebookOutlinePosition[ 472621, 9143]
CellTagsIndexPosition[ 472578, 9140]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell["The dumb approach", "Section",
CellChangeTimes->{{3.8236217782427588`*^9, 3.8236217791973457`*^9}, {
3.823662211381742*^9,
3.823662213306555*^9}},ExpressionUUID->"083f9c3e-4f49-40a3-b312-\
a3f2f1dd0c82"],
Cell[BoxData[
RowBox[{
RowBox[{"HF", "[", "t_", "]"}], ":=",
RowBox[{
SuperscriptBox["\[ExponentialE]",
RowBox[{"-",
FractionBox[
RowBox[{
SuperscriptBox["t", "2"], "+", "1"}], "t"]}]],
FractionBox["1",
RowBox[{"2",
RowBox[{"BesselK", "[",
RowBox[{"1", ",", "2"}], "]"}]}]]}]}]], "Input",
CellChangeTimes->{{3.823619024069336*^9, 3.823619027969645*^9}, {
3.823619100111895*^9, 3.8236191149207067`*^9}},
CellLabel->"In[43]:=",ExpressionUUID->"c729291b-1b6a-45ea-8de9-c4e7a8de9a9b"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"N", "[",
RowBox[{
FractionBox["1",
RowBox[{"2",
RowBox[{"BesselK", "[",
RowBox[{"1", ",", "2"}], "]"}]}]], " ", ",", "32"}], "]"}]], "Input",
CellChangeTimes->{{3.823619142949854*^9, 3.823619160216188*^9}},
CellLabel->"In[45]:=",ExpressionUUID->"d8b88f5e-86c4-4ae7-9d0e-3ce595cb0137"],
Cell[BoxData[
FormBox["3.5748532344428743120454466520221183410295739102951184284105`32.",
TraditionalForm]], "Output",
CellChangeTimes->{3.823619160425475*^9},
CellLabel->"Out[45]=",ExpressionUUID->"ff7f2323-52f4-4ea4-b5e9-f8441c3f734a"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Integrate", "[",
RowBox[{
RowBox[{"HF", "[", "t", "]"}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "\[Infinity]"}], "}"}]}], "]"}]], "Input",
CellChangeTimes->{{3.823618919989259*^9, 3.8236190346473503`*^9}, {
3.8236191637416363`*^9, 3.823619166305882*^9}},
CellLabel->"In[47]:=",ExpressionUUID->"40bb604c-037f-4af2-9fbb-183d427238b8"],
Cell[BoxData[
FormBox["1", TraditionalForm]], "Output",
CellChangeTimes->{{3.823618960916399*^9, 3.823619035088964*^9},
3.823619127258328*^9, {3.823619162659822*^9, 3.823619166763135*^9}},
CellLabel->"Out[47]=",ExpressionUUID->"9af041a5-e46e-40f5-9eb1-3cc67c589e95"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Plot", "[",
RowBox[{
RowBox[{"HF", "[", "t", "]"}], ",",
RowBox[{"{",
RowBox[{"t", ",",
SuperscriptBox["10",
RowBox[{"-", "2"}]], ",", "8"}], "}"}], ",",
RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
CellChangeTimes->{{3.823619018852666*^9, 3.8236190717035713`*^9}, {
3.8236199271656647`*^9, 3.823619941440579*^9}},
CellLabel->
"In[101]:=",ExpressionUUID->"458a9279-ddf2-4bb0-8a3d-227c65f1ebec"],
Cell[BoxData[
FormBox[
GraphicsBox[{{{}, {},
TagBox[
{RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6],
Opacity[1.], LineBox[CompressedData["
1:eJwVl3k0ld/bxg2FTClCEipEEVIa6FxPKiGRk5L4qVB9zUWKJjJkjkyJChWK
hIwhmcmUxEk4T+bx4BgrFe95/9rrs/faa+/7vq91rXVtsnKmXuRgY2Pbzs7G
9v/rf7X2237e9IV3urK+oj/3nsBVmu4fPAMRrqpyi+d78960Y7w1Pr4PQC//
0sa2uLxv4XKbKs/RKHhz29RZLt/V5FpVYmr/KQYSF1viGHbLmqJpL+42H3uM
zEJT2T2bLmrJHwtOVmuOR7K3e1GQw3utPQyXxqgTz9B4a+dm8zy2Ay91Mn7f
sk3EE8cnfcphmgeUxDaFZqQlIbOitnDO0/lA7nCUzI/x5/hq48HhQIs7oFXI
kyuk/BJKV+6/2MxecaAq4PZRbadkzAQkCOjR+g8YmDE7XTNTEDCfTTkSzEFp
U7RxSmam4j0R71yoIkX53+I39m9qr3EjRobMou2mDNUfi+ZxTUOJtdL572EG
FKf4jwr7c9MxdyL/16uNFygL9uol9vNvcFE9nmrI6Uq5q5Vq9FTjLcIvuZSO
VvhQuAQ29DffyMQh2mFtwXORlFqno5peFlnYEDNyMp/7OYXPlH3l+n9ZuH3f
8Es2VxbFCMXN2U+zccolZpW5aTElaqtbrD7lHfTT4nuJmzWUjtUqVv3kOzCs
rH/x9XyhSP4a2X7bMwd9aSH2UzFdlPM9z+dFZHJxi+eqxdnCAcrLOouPGWW5
SDj1IyfUYpIykiUaqHMhD5/cHPvCPs5TlB63UH+w54MxQGs28FuiXLkXJOn+
PB+BzbtcsjhXItf28JDQoQKo+kvVXL3Cj1/GS5mv+wvgPnk0Yj1zLbT2F3po
+xZCoKju+eJ3cXhtdjnUJfseL1L6FB6uk0EVr5LAter3SOJ8KlD0QBY8s4M0
/ktFcDjS5BP/UxEGXQmJyVzF6OToUc7Zq4rwSjM7SmoxtHRvNcWs3oW2dOFd
346WYG7F7GRnxF6IRzX9cx4pQUnJlHxi+wFY3Pav5Qn8gGyT6fC2+YNItDn4
MEmxFO+Oe7x+du0IBgz+nN1fX4pfTS2eflf1oLA7T/ar3Ud8DIl7sOLnceQ/
2PMzb2UZ9O3NztL8jdHhLWLBcboMMzW9NbzWJli8Pl1mmFIGukRGwuUIU0ja
N8vFL5RB/i6ctm0wB+VcetCwTjmMjnV/8OQ4h/MnA6bUH5VDZNrWrzTlAryP
XjTxGi7HmSXtsFA2G7zU1H7fuKcCm6u5ncUfXsLIlr/eFzsqsPLCuYY35xzA
K/59OFuhEsdLNpyWfe4MJf58gyX3SvR3bNssHnEVhmyR2fqfKlHuU3BwwNUV
V+acRR+tr4LErFWGzQ03RIwY3Oq3rcK2Ao/dHk9vILdbsUelqAr9QRbMmSEP
/K7qf113thr05OFfgsOsIgRvytOyq/FWxZ2SMuGFAVOhF/08NQjs/FG87a43
Wsa0nizl18DAuu5war8fdHZ9FecXrAXC85Tk3/uj5I5t9PqLtTD3pjeLpwTi
lVDMg13CddD0VXm4UBMKqbNKfNp2dZisHh0V+xOGqBcV/kbldQjz7/0TsuUh
vDSmvOycPyE//v0j68hInLHQdU1sqMfAtG9bOPtjNCeTUxmbG5DySUPc9WYc
Dk9dcyj2aEClxKdfA0vxUPVOukjb2ohxp0/2K7YlgOfVoim/bxPOBZWZFnK8
wF2WLNZ3NcHxi/6HiPEXmNu/1XjrzmY85z6qm/j9JXqaTupr9zQj+MQXjYWq
FHR7nJbSkmpBkdK+W03taXj0LKm7Qq8FiqEtxf6D6aBWMuL03FoQSl8KTvn9
BnX8PmKmjS3IsNSQTNqRidyETCGXm18QeYfkVR/IRkg1D+crWiuupAz40lry
obWmeFg4rB1L/NzFc3VlKB3t0SosaEec7LEmU8lyoIIrwqKnHR2fr8tNO5dD
25WqlaxGg16hglyuWAV020fDNdpo+PL0uNTSf5UweSy+32x9B2IEVv74JVgD
h03XQ5696ERWyNFxwqgBk7/je7UbO8FXMt4+lNQA59ZyjeG5ThgstxacnGuA
i49Ar4pOF4ram+U8YhvhPvhyd/lIF7j4L9jn0Zvg+/or2a9Mh757629J5RY8
UVNTUyz8gU6vGLpKQCtUhehO7D0/kFDwztUvvRVWNw7TGJw9oJ6Ptn7V3IqF
U0zh6u09SJZd+9VE5CukhHXD3G71wPB3zznK069wDP3pR5PoRZbVwjfh5Dbk
3256WbGjFzpnvbZIVLeBzeFF5VvtXmy/z317dKANUfpG7P52vcjLssrp29KO
Eu7U23uKerElMshVM6Ed/N6nXWPN+tCxbdL/SjANp68qRfg69sHgdM7xxVQa
Es9zZF+514fNZVJ1x6po2IW3k7qv+yCYsaftyF8azP+utPv9qw+Ln0v9Ymy/
Ie163vmzsf0IclT9m72vA3r2IoaS3wZwZlpv4H9Z3zGS2fGukDGAG9GC5o21
33F/7onYKY5BNBb31P0hv6PqjlxvqPIgtkRvMX/G3wkiXOPass8gmD4XTvyy
6cTevDNxfapD2DdkYrKTtwvffksu39UZQrD4ldb6jV24Tum13mAxBOac4+JW
tS7k1toqm/gPoehohpGaaRdUOm+WVdOH8D7qsY1/Yhfkl58MvQoaxuJQjs8d
hW5UH7pw7EjiMAojjkR77e2GTYBcVm/eMNzEkWqg243na956SPQO4ybvyCRx
uRsbZT/yh+wZQf6jD3zPn3ej5D+fqwrHR6B8e9/G+axunM04SquyGsHF1aUe
wh+7EavRkvAvdAScl7y0cjq7IaLXu9N5YAQdbGF8/qvpyAlNfsS3yDq/2h+j
JUkHtdX2b+rqUayUXVH+UYGOcPOZ6p79o9AniyRVD9LB58RpRn04yvKzLTLl
znRwRMp57qaMYVfivvzBSjp81Uwzz1PHkG14YX62kQ6uloAfwZfGYNuiadnZ
zrovwEDfgzEsy9g8VR2mQ8T/HVvYjzFMzrtbruIhESM3oFo0O4Y9SVXuP1aT
EK9ad2GQexwPv1TGhYmRkGT3KNdUHUcLX0rfXXkScrfgPeI1joyAudkEbRKv
1l99Jxw9jom1OaUGeiS2FT7vo7wex4CQ5OV2IxI75lceiv4yDkFR0X23LEho
ODdyam9m4HHdGzXpayQKBZbUHTUYqHvatj/bg8T+Nyo2sfoMRCWlTsh4ksBo
RNWkCwOiRpzqcYEkdK3N/OIrGejv2X5x8QmJevbgvJpvDOR/8hv0TCJhkFgy
OD3OgIqMWMtAMokTdGkdXZEJOMrdKTv+lkTrLePrrgoTaIqW2UV9R8JEwifl
mdYE/NbETmjkkzAzHeKat5lAxD8bvucfSHTPi+2R8ZgA4+Rwl0I5CcsovcvH
Qiew9mNMfXgVCasvb2qT8iZA+RDYsaKBhN0x1wBDrknk/N7p50MjwRh9Wegh
MYkEVW6x9d9JOAXQRl7umESlWYTjwy4SrtX79P6cnoTl09k923tILFjbe8jb
T8Kn5vPa430k3Dmevjb2nERbpJMXdYDEbbDxvkqdRKFL2A+eERJLdLX9X4sn
EX/ohEbRKAmv29Z2S58nMdB/YsxonATHhug4xYFJcAolKdczSPi+r6k3+TWJ
uGbRnYqTJLjO/Fr05J9C0AlDRccpEgELitvTZaaw++vrwzFMErzR5ua0XVPo
8TPOezlNIlQ9NJhdbwrkFfU3sTMkbLZmsNMtp/BfU7m1zywJrQ1NNwqvTaH9
EudNozkSwqsnJiKDpuAlZ6r0j8VjHAI2zolTmDVuKAiZJ1G+oNSpnz8F5r2N
Zkssjh0zOCHfOIW9gUUG1AUSzqRDDXvfFEwdS5sDWKzTGqJF/8l632JO4jmL
Ba1Pm3VzMZH6mubwjMWLFpdE6tcwYVzzYPkOi4dOX/9cIMnETZ15gYMsbj1x
Pyh5KxNr+x62DrPeK9WPORK5kwm3W/H+Lix+fTiF7d4BJpQUfl3qZf03ipJf
7KTLhIvWgbRdLPbaW3Pd4iQT28wygxxY9TrspKnpWzIRflv7YCCrH6ZKQ4w9
tkxsXhfHFczql7b8QqrcNSbExVYJu7D6uUOGy1rYk4mVIhHPKKx+S0iISrEH
MVFqMtzDnCBZ+/LfJ6OY8Jq3WxPImlc3t45RfRoTHInz4g6s+daxn+YtzGNi
WVLtfwXDJHL+XKxOLmOiUX6F3eggiaApP817NCboO5xyOVj6uT4aveDUy0Tu
kw8WYz9Yeu1PzrZgMGHXsO7FezqJfd+qt+7lmEZaPaNrNUufsl/a++QEppFP
zdB4xNKvUMPgU2HxaWxyWFzmaiMxUrpSZEp5Gm89LqQnNrP8I+UIW4rZNEqb
bzDLK0h4J54qjrSZxsWxdoOIjyQc4y5ev+c8DTWVnWkGJSQOP/BjWPhNoyPX
brtXHokZt+oO4exp6D1ghA2nkCCd2yPZS6bxucWGIvCCxCfbQcOpmmk8S9nv
KJVAIvF/K6vru6exnXMT14pHJAyPHMm+xzMDaYmU7b1+JNJFqgOnzs0g5L5g
KO95ln+FL1DS7Gawxn9szzpzEiF8CnM2bjPQ85pvEjzNmh9HsGVn0Axci2vP
fj5GQol5Qr06bwaVlYmHgjRIvG3o7o7jm8U6t5IBnlUksrznVHQKZnFa7oCl
4gs6ZJbkBtkqZiGaUc80e8Lybw/TuOLGWVA+Fc7ciabjypWiFWp9szDRmK5O
9KdDzdLr+waBOZiz28qfsqfj3T5+n2mrOZSu/8cRpUpHHnML7YngPC78c12z
L6cbouv5p5+Jz0Ot0bXmWXo33A/O8SVtngdNZ6rg74tuaEVUHUzWmIeoRsHL
pKhuVKlfzMiwnMejhDuMh27daL2e4lv6dh6HdKX+btToxuRfBfUewwVMm4cG
27/rwtZVO8K3hP2EuUy8+7v4Tuxr5ODJObKIDqEBhXnRDuQ0eE9saP6Ljgah
pufZ7aAqWZtX6i0j4sqqxETbr6jb4F/5/h8bQT0UK4R1XxCPneF3wUE4xtRr
qys2o7SwovjoE05CcVNtuM++elyz8tmeNbOCeCQteqOmugaWaeupK3dyEZS9
sUyjkUrEHo3+ZhnCTVxOMlvlqFgOLp8h4eTPPIQWNVfqdesHVAlvzJxS5iW4
3HYYl4W9h0i7z4PvfnxE2Ke5zJNTeYjdVSabVcNPSBnJqlCW3sFhtHd98GZB
QvPS2Vc39mfhmtb3aurl1UQLt3nawFI6uM+zdXVThYjE7KXAup2vcLvxvnfE
shBxjbLxz/WXL8HH4+5XmriGMBAvWvpdlwTtya4CB6O1RHjL+KM1eIYKhmnI
2OJaQqbJ53Jt7WNYFFFKyDhhovRkULG1egye7YzdcEpXhAjIGt0YTUZg0cOC
dnBMhJjSXtxdtCscCd/3qi8FrCMCvES+qo0Gg9Hv/1Frtyjx86eTJx8rT5z8
cDB1ul2U0N7XX1t83hdWB51KFXzFiJyl8Kk5nXsIjlqpRpMXJ8QyN49eenIH
S8u5R07WiBNdurNRSi4eEBf9Fkh1XU9oWb/8uT/aDQZ/rBXrRCSI/+58Ud57
zwUKwda/PlZKELOPnHRfGzvDUP2Wrb/zBmIjQ7xEO8oOD9L+7a5eLUno0Vp4
I5su4XrjjJrfB0kiy8q+ZauRNdw/UAOLz20kHIYG3NP8z6Oypt67fGkjMezr
dfCqlQXyb+al66VLEaklmmdKBM6g8h+vqq6BNHHGVIqeGW+CE/zSq8WMpAmv
vBWxHqEmICXUJweNpQnNv6rrjnuy8uYe8zc+ptLETaNizTWs/Knqkq7w0Uqa
eLRW3UtguwmeDB3bpOEhTWSsU9nWUnISbp9D1sqmSBNv5GNXGA9TwUEmTs+8
kiZWXCKu63RSEc7IbSlPlyZyskLoh5uoeLOK/uBctjQxLht/zSKHir5Dynzx
JdKEsaqmxIInFYaFTZxrv0oTXQFnlpokqOiu6e3raZcmTqYPzRkJUmHbPl+e
2SFNbA+TuNfJToXvzEav46Q0cV6f/dCqUWNWbnL6GzAqTRywDPjUWGAMXU3v
LlOGNCFiNfQpIt0Y7XoxRfJT0kSVnoPahQRjMC+VulfNSRNKl/9IibHy9123
VtPIn9JEn3zU46WbxuD3HdKwWpQmOjtcbSadjBEXsbhO7Z80UTc57jFkZYyt
SYLzy8vSxKVc47ah08b4Px4yU4M=
"]]},
Annotation[#, "Charting`Private`Tag$616963#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->{All, All},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.05],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}], TraditionalForm]], "Output",
CellChangeTimes->{{3.8236190448156023`*^9, 3.823619071958427*^9},
3.823619122933017*^9, 3.8236191709436502`*^9, {3.8236199376756067`*^9,
3.823619941709908*^9}},
CellLabel->
"Out[101]=",ExpressionUUID->"81bf77e2-82fb-4fd0-b41e-02dd671cda6a"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{
RowBox[{"HF", "[",
FractionBox["x",
RowBox[{"1", "-", "x"}]], "]"}],
FractionBox["1",
SuperscriptBox[
RowBox[{"(",
RowBox[{"1", "-", "x"}], ")"}], "2"]]}], ",",
RowBox[{"{",
RowBox[{"x", ",", "0", ",", "1"}], "}"}]}], "]"}]], "Input",
CellChangeTimes->{{3.823621869635233*^9, 3.823621922360859*^9}},
CellLabel->
"In[169]:=",ExpressionUUID->"21a8fecf-fecd-4a64-b3e9-9b90e0a100a7"],
Cell[BoxData[
FormBox["1.0000000000015423`", TraditionalForm]], "Output",
CellChangeTimes->{{3.823621876090323*^9, 3.8236219226978073`*^9}},
CellLabel->
"Out[169]=",ExpressionUUID->"d19e10d7-30ad-4eb5-b89f-bcd18ab15088"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Plot", "[",
RowBox[{
RowBox[{
RowBox[{"HF", "[",
FractionBox["x",
RowBox[{"1", "-", "x"}]], "]"}],
FractionBox["1",
SuperscriptBox[
RowBox[{"(",
RowBox[{"1", "-", "x"}], ")"}], "2"]]}], ",",
RowBox[{"{",
RowBox[{"x", ",",
SuperscriptBox["10",
RowBox[{"-", "2"}]], ",",
RowBox[{"1", "-",
SuperscriptBox["10",
RowBox[{"-", "2"}]]}]}], "}"}]}], "]"}]], "Input",
CellChangeTimes->{{3.823621636095001*^9, 3.823621662128373*^9},
3.823621961852228*^9},
CellLabel->
"In[170]:=",ExpressionUUID->"631cf00c-b490-4ef5-b187-157764f459d1"],
Cell[BoxData[
FormBox[
GraphicsBox[{{{}, {},
TagBox[
{RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6],
Opacity[1.], LineBox[CompressedData["
1:eJwVl3c8Vf8fx4UrKiPZlJGdFUXFPe9DJYqMhISsMrKTpLKyIzMhe5PsvTfZ
ITvz2vdeuyT1/d3fX+fxfJzP+Xw+5/15vV/v94fXzEH7MTkZGZkYBRnZ/5/N
W7M/F919wX0UIzqvyiMnViiPd2PzA+YptCZuzEwu1u7jPG2FH2Ay8n6uu6XJ
2a/f7ZXb9ocFVrG1W70sVzhw9anvLIPgn+r0J5ifveIy9klNRuc95Kyxt/W/
pr12ZvaCVjNLNLjz//fC9dY7eZWGiZ53j6OhIfzumqZTorxHYoCyblk0DEtf
yBKjL5Rff7gov6H1AVLP7UTLyQzKN43HCzKHxkDiK8au7ysnFWyHaI6sKOKA
tud3vSTLM4W04kpXGe04YI9f6mCy9FQYC3+8/TclDuyT3Lyz/IMUlDSalyOR
eMBjr+Fz1VIU3CQcTI1C42FVyHxrISVXoYD27LTQdDyME9QDvm2VKuDw3bo7
op9gA180kF3VoMDe6/at7uUnyCm/7LHl0Klw97OgWkDXJzhl6bFYU/FNwTd4
pEOLNQGsNVs0CYGTCjXWPopcTxIgupvR6k3posKWilTdclkCiE7cOsCQExQE
hWdkiykSYVIxbvKRxL7Cw+Mhxa+0E8HTLyuMquifQvjyVTHl1ERQeKd8SY+e
CtvRvpLFsJUIYhcphpfd6bFKxcqtr24mwfQ6L1XYHgu2ISFzdjk+Cezjlf9x
xvNgrwVSHmltJUHY7ELd2pAQtvKZOVv9zWSos5uPtvophb30qOWS8KdkcPTo
0kn9LYctus2rFbWVDNEULFfVhlAsC2eT7H98KXBM/6lCpuMt7Ej4ZluZegrs
PLA/FnZFExtJxX3Pxi0FtHRXfyeq62I1X9+d505Pgd3KAdePbsZYup03Dt/7
UqBkPGmws9AC22f55W/wQQpctpdxeLPxFBvyY/odej4VmhJfsKvWOmNv3zvF
8VM9FXa7L4dT/XPDUn+Vz/nslgrmtG40g9Fe2A7kqaxpeirUiwz4vVzww/qW
xbex9KfC9bNVXfR1IVgl0W7t3oNUYGBY2N1si8SSpfye8z6fBkIP7f8YPYnF
NjKLOMjdTYOm6NQN34hkrMc7/b94tzQwOOjN/XIvA6twLPBdWnoa0FlMkJ2i
z8Meulay6/engaR6/z/N1CJsFX45m/Z3Gri7pr1L5CvHupqxyLaeTwcvaa7R
0sYa7KXxm21ud9Ph6D/didypJuyO+nNtiZfpwCnQbHZdtB1b3Joxt5ieDv6r
7t4P3bqxjldH7OP600E0us/U980gVqKQ4u/d3+nAdlY5IcxtBIvnl3lHyZ8B
1vIOnUbaE9i8eDP2mrsZkHQ3TNZbdhZrxRCZ7fAyAyKeVz+kElvECvo3XxbI
yABCJydO7+UKFvdnq3WyPwMip973HWbiselOPNrhvzNgb8EmRNRuG2u2ojF3
kz8T6GJWleL59rE8Rp72f+5mQoKil7Lu+9/YmaGCo6KXmRAQGuXqP/APm6Ay
E/wkIxN0Fe9WHPmSIwYNtOxcA5nAPypcayRDhbBdwmZ/+50JxxkcgoyDTyCj
ubaXA/iz4MuTBfr2DHokmjuhVUEjC1zZk2t9BxkR7Q89Wjsvs+B73z/eBU0W
hOHkn9nsjCw4zdT/8IwlB9LvJWpvNJAFb2Rb6EGEGwn9+eCI8TALxPN0iqX0
+RA126DgLv5suHtfo18tRAA5sVDF5qGRDSZEJvFGVATp0lvNknHPhq/10vcV
AsQR/z7Wy2sZ2SBRrJf89epF5Mb1W61JA9nQx7Z9cZv2EkJe7aqlc5gNZ3jz
HK3U5JAmiaxZGoEccBfHbSyRyyMeGd/tGjVyQB+N5MtAAVHgwBy5uOfAbriR
mdJJJeQw7FKwaGYOdOvYyHzquYFUYSzY5gZyYE/ayaRg4xby4lVU1ofDHFDP
E+r+4XoHkd1uuXRHIBc22GtFSl5qIHtPdlrINHNBxccu9zeZNlIyzatV4Z4L
QRwlByZhOoijttbs08xc6BXat3zgqIeUyIjOq/XlgoFLgTvmmQGyx0SxKL6f
Cwpc//JT+4wQ2Z9TOLqzeUDekb/o0maCuI2VLW/eyIM+uTOuomNmSE1V6Oqg
bR6UT86cOc3/GDmKe7JeHJ0HIVIP+0qKLBHkFeAj6/Ig8BMl78m3NoinIRvx
GS4PMF2o1PtUO6QJu72pc+oz7ItcOh8b5IiQc3dvX770GZJnfkxUODgjN8jS
d1kMP8PH3SUPkQAXxH/+1f6vt5+hZelxqdKIK9LVovNr/PNn6A6LWc58+BI5
kSH+u3r4M1C9ttm5w/caUfOj+hP/5zNkLoUtBYh7Iu+fzB69Op8Pt63L//op
eiODt6r+Gd7JB6G/N01fML1FGEUiyJBn+ZAfSnd1nNEP0TlhQ879KR/ynuvJ
8aABSMyGEiVZaz4Ek1t4HksJQsZ7Oanm1/OBJ7T141PZEISjYO94C+MXCMvq
d3l/PAxJdMw66Wv2Bch62ApoKiKRGS1P2sfBX6Cex+1MvG80wiOjT69c8gV+
l+v9UXGPQcyYLp4WmvwCPkGnNlbiYpGMfZoz1OQFwDwt2OC8HI8sjy4wrYkU
gEMKtfdp80REuKqWpVurANbC+y+E0qQg+e52HCGpBZCwc+g1SZuOEB8qc9l9
LYDCbZ78/foMRArLfe7udgHoXrnuYhuWhTifO+CWZC+E0tnH12x9c5Cy/wZ5
GRQLoTGTVGSS85Cfc7nnt60KYTKYK0pjIR+50uIjMBReCF8rAiL4NQuRWt9L
ItFzhSAxw7Th9LAE+Vfi/+EaZRHICRS7tU6WIsnnn7FcYCuCd3HB0ReflyNo
9KOPnGJF4L7SQKcvVonMU6qxnUKL4KJdLtPysWrE5/mVuKN7RUCW+2Ak6VcN
cn6Zn4NgWQQyNbZ8BvT1yOPOv5z9YUXAfozLTWu0CaG6sp7QkF4EwlOa2YKc
LUhWzujZwsoiuC5Dl8fxvBW5xd6alNxDWv+xPJ32UhuyGlTIHT5bBE+Yc+v/
2nYgQYefUrx2i8A1V0f/IW0XIvo0kNfpeDHcVgwLK+/8itiqmZ7XliwGrmFT
tv6wPoS2Xj1D6XoxHNZKmD9JGUAKxK8JyOgVwyy7+ZHe6iCyRXdGiMmzGHpe
z3Po3hlGIjz/y6aMKgYJZh+y8Y8jiPTWhvB+VjEURLfQ+/76jrh8axMd7S+G
F9flCij3xpEdQfNgs4Vi0Ayj3bIKnkQcXh9bJ+4Xw32Hy08XZKcRGyEk9/jZ
EjgTyKxxcnQWWXk9TR0tVQLXJgaiHi/NIRZD7lY8N0oA5XtQPYefRx69qRS6
+rQEHgQd4KwFcMj00P2Ado8SGLKnpc3QW0IMhPeWtSJLgD1UZlYrbhnRGZbK
sqkugcac5xYf1NeQIeEBzK/eEnD4LrL1pnkd0fCwe/x2rgTq3sZSeN3CI7dF
8vgTjpeCFP7GUkHIJtLloeIrzFkK5w1WKfOztpCbI8uLZRKlUMeM7QyK2EZQ
z/PpffdLgdrA+u6rwF2kYaSZ3MC6FFxmrj80/7CHyIuamC2/LgURywjsm4J9
RPZ7Au9/6aUwpRlMRX/8ACkTlfd+V1kK6yUt2c63fiMXvSbmWHtKQcJSNmf2
wyEidoElRWqnFEZ77v/uN/6L8HqHnTNDyyB2vdSyZewYJI2KexDvlYGzrd+R
RiU5cIn1/nC3LAOHFXUEk0IBLGPUiVHvy4DvnXMH63sMRItlH3GnlYHRvcGJ
09FUcNrnpmF+eRl8d6LE2mUch1PiPpzt02VAcNL8Eb1EA8E+PK+0tsogpvAv
9SnWk3B8vGHyB0U5OLitWl/TOQXkb//E/RQtB/rnFkWCu3RwMO7CJvyyHDTc
Bi3anzCCq8QZt7KQcjBLaJT0kDkDu2+Lx9CUcuBjt12IPs4ERAlizIPOcngQ
Xnm/sIMZnvqG/FyaLAdnHMdl7hIWWJ0Q1XUmlsNj5z6F9kxWWPS1ZH7HXAEl
om4bMVnsYDKJec4qUgHR77yyZcs4YEYyYyRdoQJ4Dj5V7XVzwvjkXFSdeQXc
iX8q3Mh0DvqlDBiJxRXgn7CAkw3nBaOvFvF8nRXg0eKYdaTOB3hTBz696QpY
cf+rrMV4Hk5E+Uo3UVXCFyaTQ71Cfoi7EFazx1kJlDZsty0CBUC4LU5J5GIl
2OvF7ZZaCcKt/QLtyIeV8KiqYGoChGE0tHqyw7ESGtonLjHIicATwTazP36V
oP0t6+xtWVHw1ZtwtiiqBO5Cte1wdTFoqqaIksVUAZZ2VSzAWQo0tOk4n3JU
AfOlxJn/GC7CzDpberJkFYivs6uplV2EI06JUmqDKghh6TAepJGBqx76w5MF
VZAwdedq3u5l6GIxf0jfVgXGtYZUi7WyoF9ot3h9ogo8zDalJ4LkwHXOZzef
ohrU/512fyl3FUqUvpzx1K+G61+Tx8l+KIDiVOWnUrtqsL2xKRg3ioXBZy3n
V31I77mrlzdHESBmjMlofamGW54n6aRcUbhwnFznPHkNMBlwONFUKkFN8qlp
PdYaoKlFZ4o1r8PtK6wWIWI1wBfC8uom8TpYWYu57OvWgJNASOiLyzcho0c3
uvNzDSTw9Cad/nkLZCxMuY6aakBgU36aq1gFWo6eZkiN1kDiEgWa56gKc+Le
ZXFktUDk4eCK/3sbuMI/jzy9XwtOHq4RfzXvwmfhCsMUm1oII6u6qSKtAdea
m3AjnrUw8y/n+sNTmqC/830Pm1cLB28Sv+QcasIHHTJmhv9qIZTp/vm1X9qw
nDarO8lYBxhtyXLZ5/dAdrshNkOwDnwkW1lP7t+D8dA3XFfv1gGvebe47ZEO
cHUenjdLrINOvfZU2gt6YMc8+VisuA5ktRmMwur1oMG8OvtnWx28z6L5N6et
DyZkbhfe4evA+fuqPbv/A8i4ui9dJl8PrvUCQqv0hrAfMOLioVEP6k0Xhrar
DUF5tLRCxbweHtgaxbBaGsGqs/O1H8H1gOjkPprpNgax/E1Fqsl6oBvNoKqY
NoGyc+ua+m4NYDG4dSP1rAVg7L5G8oU0QLfgzXOToRagW5szgk9ugAjGe4rb
/1nAgZ6lvndnAywKO7AOrD0GhQjcozyWRuAOeNI7OWIJobOtqS6ijbD/fDTx
o44VzIinLyJIIzhGGvobjlqBV7ep5fDjRtDPrOw+M2cNrRSz9kdljUDDKMrb
d8oWmLQbijq6GiHm76XIR4m2YJGSuBM+3Qi3hrl7+aXs4O05iVpG8iY4PZQv
bmRoDwK1eCRbpAmELhik0Vk5gvWetfKgWxOINvR/G2B6Bhz1Ya+C3jZBlGNy
R+aLZ9DrV16k9L4JQg/9en5PPQMJ1mOc5emk70XTWUVzXWDvauxmXH8TrOJn
beh0XcHDsyPWnL8ZqFMK1/4evQRJFXwfl2QzWN9n+GD40h3mGBgpRq82w/es
i6tOv9zhepqRnYpGM8y12FX1H7wC6vY9RXH3ZlihozUrY/CAqBPnN34ONINq
+hd2ATVvuDGswlM02Qy70ssvZ2K94ecn+/vWS82AvaTlF7bsDfriNY1Th82Q
R3U9b9bXB85pakU3CbRA+Dn9p9Z9byHngyc2+FUL8Lqdp/kc5g8PjDOdr/u3
wJnXp4OK8f5wQqgn+yi8BVK3trN7bgeAXRULo0NWC1BvozmKJwJBeurL8r1v
LeA2q95nERNE6tOmw84KtcKmk28UZikEpqcMyjIutoKrC36HRjsUjmImxi8o
tEK24ragdFMoKNCO8V7TagX5vd0He8nvoeb3t1LdV61Azea2XC4bDhVDnWNh
A61gFNsr3MwTSfI75SOWyVbIDf9s2mgSCfsq7TxJuFZwhoKZh6mRcKmxxTr/
dysErdZnOfNHQcnn+j9d59tA794jXTHpaCjwLeUmd2sD71Zrxee2MdAP0jeC
fNrAx+h+eFFFDBAOi6wYQtug70vo9IdjH0HMqaD4XFobDG81KZnGfYQ8o9zr
8r1twOp++ifdt1jIlk22dOFph65QtmOmpp8gZTW4cLmrHawvikTc8kuBlH+7
RVVD7YD97rthXZICqUxGJcHT7fDnQMthdjYF0lDJcomtdijSzd1qk0+FjNjh
GlfWDlho5a3w+pUK2SpcHVRPOkBQx5xj2yMdCnK+TAtSdEKR2U13bdI9raCB
ZebgVCe4DB77ePd4NhSOeM52s3TC6pknwpmQDUVk2gv2op1gkSj72qw4G0r0
f65UaHXCvWL/+PNxOVBBDbvKKZ2g2pwkin+ZB/VWgzSW2C5gL61OXvUuAOUU
B25alS7o/s2DU+spgMExusul2l2wk1kTyMdcCAvK6ib/WXaBZ55OF19eIVAL
dFd8jOwC03/7Dg22RaA932retdIFWOG/2/OkPm7lQWW9SORXcBFWc9hTLgPH
CN3hgYSvsH5KnBrvXAaHXfurz7O/QkzsryvcyWVw6solppa6r3DUfq3oykEZ
SLMU2xisfIUdnt4T176Uw+uhPNZ3Ct1w/cWRuT9fJZy+k+SEX+4GRoEOqppr
NfD6h5DFj61uiPD9ETxmUwPLDsW6/YfdgD9m2SL6qQZqotvki+h6QOmR9nPn
oxowm1nHuMj2gMXBO8rA5looc5KL/+PXA8s9Qewj9+tBP3aw5YRAL1jyJdt9
uNcEKUvkTCLmfSAeOvqC7VcruNks7a/Y9sEAW8uhMXMbaG12jmW59sGdpjWy
Jpk2IP8d8ok/uA/unR2J++bYBhanWM/zFPfBuL4CvQy+DYRkLkiz/tcH+Mox
s82VdijwuqeJ+dQPTvxR2eq7nVDPmREyPzQAgrgW3ucPeuFcS52Py48BePVo
OPuNdy94WX13O746ADf3CNyZub2gVE71WPzvAEwn505oHvaCRXZZUt2ZQegU
amblTOyDNZT1BTUyCJdKAotfrPbD7rNpoZTIQZJPXPrL0TMI1JNPggbkv8Fi
okjIo6xh0CTvXAt6MQRpWckXv8VOQOnDjcGOvmHQXZNQ0382Bz/IJdKbz36H
OAFhbloHHAhQ1TGy6Y1C3Jh55PMXqxB4jY5mO28M9o8PffHmJwB1z+cEgdVx
aK2c/mw+uwUZIbSV0zAJT1TnWTde7MLxhNA6rYgpYNpjO3+c9ScEHH/wiK93
Gr7dMf0lGXkA5zNWA5/zzkCMGQhEXP8D0wmtGvwms3BCubPEo+EvFDjNeKmo
zcGjUVOrF7Rk6Gtx7Nmm/Tk4//G/klM3j6EPPioGLkXNAx+WNZNCmxw1u/iH
LypmHmzXG4j6uuSobU9ZPRo3D8aVdwXzDchRr/+Edj8lzUN1f+d9NXNyNMuS
zlg7dx4alFVdnrmQo3tyUzKNDfOgjhFONIghR9+Pu8zGrs2Dr5nK4KtxcrSV
PUdWDRYgy4H2g+d9CtRqS06doLgANFQJtRceUKB0nZ3m728sQFm26u3vhhSo
vstK+KDqAnSWnY7hsaBA8f2CGzo6C3BJLSD1ozMFyuKbmWxkvQB9sud1MWEU
qA0xjdohijT+2ieHvnYKlL5dmpshZgH0A4rKeb9SoGWfWi4Xxy6Aau8jxme9
FOh/Kgtmu4kLkPmMGUM7TIFGp/PWv8hZgPSP5v5n5yjQBv0UJ6/6BaC1Nsgt
+E2Bnm5NnIxYWQDpDMNmD1FK1F7KKgezvgCWqZdZEHFKtDdR2vUlfgGUPOVi
DyUpUX+3ztOm2wtghBngeHqZEj0S31aR+rMAFwOGp0RQSnQ19kZlP90iFDbu
63Ldp0Qb7TeiTsougm/m/QftrylRrulyU88riyDyOXJfwJMSdVf1kty9tght
fw1d33pTopf4mXsmYREW7HW5LvtTonkTCHme6iK4Gc9VOYdRojE3Ih1VjBYh
cfyy6P1UStSO46q6n+8icBppume1UKJiDKd5JAIWQTEn7MPXNkoUj1nbGQ1a
BPzKi6LVDkrUdjs2VjhsEa6p0Exw9FCiNl0Hiz1xi3Dl89Pku8OUqJVrlTtj
4SL069tqbi9QosK24eq1xYtgc9SyP4SjRFdMrXgsyhbhU/Mb66JlSvSJOltH
efUibL0pLTRep0Qf87udftC2CJctjpkEblOiZkNyOckTi7BkORNN+x8lytdF
/0plehHWYtXz8sgw6Hz9ivr2zCJ4MZ8zvE6OQU1yP+4q4RZhrux2oCUGgz7y
+oUsExdhct4vz/0kBjWUqPwuRokDWXUB7a8sGJSLPyznOxUOLHTDS86xYdBp
dstXHjQ46PKteO3AjkEfYlh5B+lwMKr3QpGCC4MaTLvaPmPHwUMKX5ldHgyq
FyxLUSOBg5FFDHFbBIP2s6iGXr6IA7l2w7F1UQyqnP6QtVgGBxcGT32YvYBB
Zeu8LuRcwYEWGeZtnTgGZSF234tRwoH/zrm56xcx6Hv3Hz9O38RB4L0T5zml
MSjV8S3L0Fs4cFY2FSCS+Cc382tfNRw8cR01CriEQUe1HmU46+Kg/6yzjr8c
Br074yRB1McBMqyhKn8Fg3bY+FZZP8TB7Z2E1Q0SV7zN7TUxIc0/QcGFXsOg
Egx1utNmOGh4Np23SOKshP45vcc4iOdOm/CRx6Afynf37trg4GWr46kyBQx6
SonKs9sWB3ov5E4qYzGobz8bjTLJ10arniQPk9hlBXsW64KDxxdevplHMKgO
e+DNCx44oF9t6uhEMWhvZvxAlhcOFte5hWQVMegN6S8P+N6S+OTk8RQSX7oz
ZMceiANO4hg8UcKg+WO4X9HBODjdNhbRQmJ+i1/eDKGk/1/Oe8BxHYOeecMV
czwSByckrqN1JA6mkeR5G42DpKdtD6huYFDyGMW8fzE4mBGNmLlDYnc+nUvu
cTj4/e54bwiJdwqeNOx/wsGLWj/uryS2kX+p4pSEg8o961Gymxh0ofPdED4F
B/d9GYkyJDbQSTK0SsdBW/cPczMSD80VLS9m4iAZEb0cQuLbdq2Oj3JwYOCg
+6iYxC2/vx9O5uGg6njz7DcSX/Nf9dX9ggPNwa5SAolLGP/QDRXi4GpJ7zil
MgYVTaaNUy/BQWr1BTVWEqdd4Dn/tYykF7abJwVIzFEl/eVGJQ4shd3OSJA4
4sZNuaZqHBxz4DeTJjH1N71m+TocjCfl7l0ksZeRzZ3KBhxQfrLsEiPxwdrr
79LNOCh+2TDJR2JH17BHBa04GPCnE2Ei8Sp52ppIBw7mlKuKyEhsElb2LLML
B50vuRxWSfsd5+z8y9ODA6+zCaY9JNbMmQhI6MMBjVOsXy6JLZxaQk8OkvIj
2jD7LYlfXPsc5T6Eg5+cyQIPSBxMER23NoIDjpXVvyIkTup9naw/hgN3zlbx
n6T4l3x4nNk5gQNa4cqSehK3G9/9LDtNWq/zbrAXifFb3JVM8ziQvBwhdEA6
b7Ia6vq3iyS9n9Jazf+/Ht5ut+ws4eBe6My2IYmvMbf2D67j4Lt5+PcCkn7U
Zz6PAAEHw6G9eZokNsmOnizYxME7BqpGAkl/AVefLIfskfyBJjOencTx5Br4
Pz9J+fdSwCiHpN+CHrkdm984GFKY0JEm8YgRzT+VfzhQdWSelANS/AR3KKrI
luCEUZ1ZKUn/fzYnaYQolkAmYZhFlMR8PvnMGOolEFH78IualC+ydz5wupxY
AklOM24nUn6pMnnwLp5aAux/i1YjpPxzzNIQbz69BONBWnohpPxs6N658YZz
CfZdj4yaSfk+FDV1G392CcpYP2qQk3jJsE3zIc8SkBUHayCypHzd/GB4VWAJ
fiKTjtkkvzA4c9VlX2KJdL/t28eS/MZ+mtfd4uISSHXNsOhJYVDvzBNewzJL
8FmcWeapJAbNkZt+V3xlCfxqH9/zJfnV/kPPdDulJWg0Wt96SfI7agGr3Okb
SxBa9zzHRpjkp0TNwju3lsCCnU9FVwiDKnnx1YqoLYGX8FPZcwIYNDyjfWjp
/hK8vtIX4kjyTzHCSXJj6yUw3lZnOkvy4zSEnJbcdgk+wfGJBCYMyhp+wJpl
vwS4GFVnljOkfJZZEic+WwK2g2yWX/QYdMKt4YGHxxIcroyP2NFgUH9Kp+Kk
qCWomhlJkPhDqq/3LeuUYpbgXG+f2u3flKhztlHncuwSVHfSpD/6RYka3b79
QyJpCbR3f8m+2CXVz/DzNI05SzBiNPrlwQYlOs85ajJbvwRu2imWaxOUqK5d
79O3TUuQ9OimWNIYqb43tLgKtS5BStI3G7XvlGiVaWGIQ9cSTBuP3v4wSIm+
zw6s+m9oCc4vxNwZJ9VPBRl5Bp7VJegyNyXjLSHV49spjSaMy5Ah6qPy+S0l
2ib9zXuMaRly/KN5mrwo0R0O8ht3WZfhg7g65cAbSvTuhnmXPNcy+G2wrUy+
oEQxIULfWASXAQa9l/KeUqIufQULvVeXIWS2TERUmxLV1qzHXDVdBntqqoA4
DkrU6wqxs9CctN58vQGehRIt4OEOFnxCmn+xj/zqGUr0xLYXHdPTZSBIzX1t
PEmJNkfcYN10IY0/+4vm7REFKjXcI5wZuAy5p2npZ6dJ/ZTu5J3TRctw+/b9
H4uxFGhKB9uZ+pJlyCzD6tJFk8bL6U1alS/DjuGO0UVSP6bJNmLVVLMMExfe
MZn6UaDhk31+Du3L0DT75M8TJwqU4VFzY9/kMrC9N80/qUKBMlrlyARjVoDK
L2GpeJMcTRtfPrxMvQKyAqOS+DVyVFpVoGX+xApo4IXw53DkqLZomuZVhhV4
dsOu2o7UX0bhP9mtcqxA5PeKhKImcpTJKSxbWWoFZllV7vWEkaPM7q6cFA9X
wKZKvFhHmBw19t6M8zZageXreVL0fKR+NtCKjcxkBVTnerlaOclR2Y8GzH8t
VoDMWA/LSk+O6pch9L/sV8CWwq3Ycf8YGk+kotjwWYEqHy/6taZjKI/5x/Wh
vBWQ0ZlK/Kl9DBVVq6pOO1wBc40LYxPmZChd+PEypr8roFghdchjSIbuDOsW
+P+3AgH2flJmOmRozcP9NBvKVbjGyrgwdIMMVX0qHSJNtwqWL/Oy7QTIUKvg
fOMW3lVYp6b9wtLzH2R1pVAsqKwCOduQ71HLPzh/M0id9+MqSMoaLyKRR5Cr
W8c/j64BI6iQl9z7BQjLh+67s2vQZImciqHdAVPiQOnf4HUgBKvYI9WrYMxy
ZyM5dB3eVvOI8r1fBQOkg08pfB1sh88ts5mvgvb7ugj/D+tAkd8Ui9KugpJE
rgN9yjrQ6Wv7eJutAJ+9zwXe8nX4hRjcoju9DDjCpfQbs+sQ4O7a99KL1Kcx
F00uz69DDsPw1/v6izCNvcAYhFuH5Qiv8LtSizASyuvdt7YONClqRR/nFqBd
nO6R7t46LNzQiw1XIt0b7FY4rGg24F9xst6jU/NgSYiLfCezAVHkPWEbzj+A
K1/Dm0x2Ayyp1hS6JX/AoA3G8fmVDVANMppfx0/DlVVH9UdY0vdpYe27ltNw
YvEWtcytDdgqtPsqaToF+eP7byYNNuAaItOdbzIBO61a1iI+G1Ah6kvmGPMd
snyO6yf5bkD/vcZibcPvYKBYr3wmYAPOcfQ9f8v3HVobhPn/vtuAuCmc22Dh
CERX/zcz+GEDAvsjjnj6h0Gu8IuOW+4GfCKLS5pmH4LX8TRo5+AG1PNFV5DN
9kHkfK+M0/AGnDjHmH/nRR/kCocLcY5ugOLijrwJQx+MVrDQO0xtQLSzUNbK
jV6QGj4/w7K8AWW/quheVnbD0knk9ZM/G5BwZ3RgoaUTjrSPOTL824DOGdas
dctOYIxvM68hwwP+98/fjbSdgArfuUNHhYfVllmH4YcdEH9Dn6OCAQ+4w+Xx
pGPtcPeNcxVGEA/z1mmBzc9awKLtcn6hMB6OwnooUwVa4NXJ38kPLuDh+Ndq
JvXxZsiJ8wz4IoUHi1oeyhNoM5BXhOjel8fDFCN1KD17E1QQs/YyNPFQ38wp
9mu+Fnov26zevYeHPs7qyT/GtbD4Wnz64D4ewnflE0J/1ADDybJWtYd4sBQN
Wbo3Ww02Qs2Re4/xoF6hIXp7vxK4Taakbrwi7f8235Mi6zIgYiemmt6Q5rP6
PvoTUwZ1nGP+WC88EEvzl9+ml4LB2NC0rB8erpR+k7VdKIGYu92BImF4+PtS
B9QciuGxWJdMVgQebNOsfZOYi+HSiY4Zvmg8eMuFa1yoL4Jv7c2XuOLw0KLi
u0CgL4JT2Oo5unQ82PFM1/YHfIEpjsp37zLxsPEf8q1uJR/yDspkqXPw8MBD
Tpbydj7cKisKOfYFDwM5G6fDmT+DD+l2sFeBh4Cnl69d7MqBA/aP4RPdeOhY
UPXniM2Azl/R8rp9JM7O4OngzICY75HLQwN4eHr7A1dwajrIRLxX6B3Bw+/E
FXfV0jSS7/uvNszgwWdT0ldoMwWwy2+j5OfxsNgmI8j2NgVOtXkjVYt4CJ2m
rz1iS4FczzfRxat4mMygbD37IQlwP13QjB087HAcrKcFfoLSEecNnv3/j7d1
r+qPB58Sx5jEX3gQd6jCWLPFA7e9LT7mCA/dFYppzypiwWDJPDaIigD0Hk8+
sXHFgEir6XUqGgJo/dKk2XrzAQ5SHhF9ThIg46A1knExGj4YPrzxmoEANmOj
lgyVUTA4rL1lz0EAWhXMzDffCLhGFGSd5SKAypPwWnfyCMig/oPV4CaAswga
/59POLzEpr+T5CeAR6WnGoEtDPiydwW3xAngWxZVzez1DkKaO9VNpAiQ72vq
TJMdDL+mPrkMShPgFsNDXY7hIOhluNFSJEeAdMcJO+VrgfDcPdrISZEA65zK
8l+U/GAu2sp3/joBuE5MvncI94XbhQqftZQJkKymnlC+8BbO4XAHF+8QwIrp
O7t3tA90qV+O3tEhgItU7SzVBS+QtqKpNdMjAHOrbkvhZU9I8PkxP/SAACLP
6yfEz3qAU6WfZKkxab9sX79doHgNHLyj3c+sSPFyo1O+KuQGvvJ524s2BBjR
JBxvV3wBm/c92HTsCOAgKh5Fa+4KbaRqfsmZAGNY/2XVaheQyDwMyXAhwCJf
1NcDwjOIbewvZXpBAI73g4WTws/Abs/12P4rArQ6XhwpKnKCcbo7wo89CPCx
k/eF2V9HUBLh1vjuRQC2gZfrzRqOwGLcmVDuRwDMyc01hhP20NTJesU1jADB
J+rSnxlag+jChvFyBAFEH0/R1apaQfRRo59uNAEiGI7dsUItwfqi1bBsHAGM
M78mU6k+hpE7CodZnwiQYJLW4WNgAcgTBl7WJALwXCTMuruYA+OnKvtfaST9
jLkb/20yhdflIR8sMwnAonRyk23XBFYGTOrGsgnwprezd0jMBOooaU5U5RMg
6uyXg1hbI3hs/yDTrYIAD9dyay3k9EH6tPy3s9UECPpjJcm/rQv/Srn+ttQS
wKTj8R5D9X2I/T2rQ9tMAFOOm6Ifne6BRUKzV0krAYDDYEDFVBsuQnq+XgcB
7jsGXNQ30oJu3yeUaT0EqFJ6PFm2fRdihFSkbvWT9MRYGL0loQ7m3SKG+EEC
/CfENh7lfQeO6AmlsqMEuNlnPLplrgJdJf2zU+Ok+Ht8VhDcUYYP94tOek8R
YP9K/euayJsg8emZec8cAWY8ohVEqa/DIfZ+mOMiAR5hknOyphShc062lnmZ
AOfWHte0NqBgInjIaLJBAHeWLJocMyyIfZ1CMEQCyM1lXWf3lIffT+tt8rYI
sMarf3m+/CpEFnu17O0TgGCs+hZvJwvGOmbEuANSPj7wUz5GuAQXfl3ngD8E
GPXkekTvIwMHcQLKuL8EyNsbHWUSk4Y2hePOQWREuCnrwsC2LgXhs6uJEhRE
4HNT1tVwlQQjn+6vwxgizMlkfS1TEIefnaG8504SwXGqkV1fUgRabBzUW2mJ
oGKAx4O+EITRar20YiBCSn/C4+YkATAsks6kPUMEgVQne71j/CB8j+lbCTMR
xi9dTnnnwQd7+/tHemxEWMec+b3PxAvNsWPCfzlI89ek8b/05waDmXivWzxE
8C/s0n3RxgGC3q/z8XxE8Cn3sqVOZYPd88bjEQJEKBQb2gtPZ4GmDqCUEyYC
/ewa//xXJgix5pWaFiWCvpZS7fzpM/DgFIWhtzgRurzbDaNengbBQlyAoBQR
EseK5zn/0UPjXvas42UiSCylHD6sOgEhH4NOslwhwuzNI79FE2rQv/ZUrvYa
Ebg4jKKfyVGBwA81cxMsEbh7asVOX6WEHU+JMAxKhP7hM+yrZuTQwMdQm6dE
hHzbpgTVEjJ41769rHGTCNcXtvR1Vf8h/CfLkfjbRDA5/E9q8dsBsvUlxgbU
iSQfmBr/bvsTqdd0i8FpEIEzau4Vs8IeErz7oCVImwjX3qtKfpPfQXRj5IkS
94kg7KDh6WS1hfBdPcsxokeExuh/HGJvCMjm1L+bLw2I8HTUw0nQdR2p85hz
OmdEhDc3HyXPlq8ggbwtia2PiHCfQJfQLruE0D6c8ZI3I4KouJThRzZSyxl9
aF5qQYp3lMOrbNF5hKWf5dYFSyIsetqkvB6aQT4dlxFNtyZCjO/qvebUKYRH
UYOW05YISjFUa3F140iG+9OtSHsiUFSkTm1xjSIiZQHDJ5yIoOX+dZqtfhgp
IKRX+DwjgqSdUHDnx2+IjFBT3OFzImSJ78fyaA4gVSbTr53diNDM6l7bMNSD
KMQfPFp3J8Xr4QxxraALaR5mum72hgi0qSPazdPtiDLtRcFJTyKwG2FKf+u0
It3K6jTaPkQ4GFAfH2FoRjS8rPFffYlwJ0NFX5qpARmp9htQDCDNv3vJ8Ixr
DfJgN7WkOogIm4LBrnUKlciMWMOHiyFEyNS9edfNuAwxezLplvueCMXLo7eN
h4qRleSfD3kjiPCs6jxDME8hYjvBCHFRRNCt5/9mw/sZ2WKU5DsdQwRF/4XA
urhs5LnaHUxQLBHSaeLNT7/JQA79LFfJPhEhQqW/i70pFfFsfNvjlkiEsTLw
O1echFD+Ti7YSiadD4MNraJ7PBIoXRdhlUbSW7uxTH13DEJrO+4yl0GEmase
J3WjopDIzD09/WwinH5EefVuRzjCMssgP5hLhK/2zFvhBSHIJzbxcyr5RLjw
tumqjmsgwqOteqypgAjW4kHEfwO+SOa7xzi5YiLI5qrQt+V6I6Lt3p2FpUS4
GGU36L3yBin4l5gnVEEEvPJ+Xgq8RGSu1IQmVxHhidbQyxal50iV06gjay0R
Rl7VSaRMOiEKn3fuhdUT4XUrNpx2xQ5pxtHJHW8i5Vv0De8LB1aI8rkLHJ4t
RJB7zuqO6bZAevRu/f3ZRvKHOmWjNkVTRDPCfM6+kwipj739hsIMkZFuz9bl
r0TwvORaPR+mhzygTMgy7iVCmckL6X+P7iEz2Kqg0X4iIOW/ZlRbNBCzFyO2
d78RgcjJc/VwSxVZKdrS6BgmQtipiv2yszcR2/VTMsgoEepg1jCWRxHZPi/C
UjFOhHtSerTTBfKIq9HN3+JTRCDjiGJaU5VFDmNMpzN/EMHOH+Rfzl1EPAff
NJ6dI4JTeh2LSbgYQnkiPu3DAslfOD7KmtcIIoHXK/xol0jxeCCp+IacD6F9
M2Tlt0IEqzSTvC9PuZDICuKdv2sk/42k6Ku7x4Kwbp2QfI4nQtHCGS8B19NI
gogQI4FIBNZhwaphyZMIr/n1fYttIvz5UlD6nAGDZCU8Gp/eJf2//KpAK5YM
uTD6qlbnJxG8Xk46eG8cYIvoY5N6D0j5fqqHgxmzi72kWuZ94w9Jr89ajRXn
8dhqn0GLur9EqECchV/0L2GLhv62vSHbBFflsjKL3FlsFt8FAaDYBDQmyG/5
2yg2wVnfj4xqE84+ci+TuDeIjWzxW2qm3gTJBIFJe/FObCBj6c23JzdhYiQF
DdRuwnqYzWXeoNuEq4S8UmJvBdalhJaK6vQmkAkGEHmjC7E25PJPOs9swkZi
Vy4+PxNrom3VEciyCV79gmFUrxOxumkfBG+zb4JIdLUs/48orNpOi/9Jrk2Q
iDcuSvcKxiopbS33ntuEG0kd2sxe3tgrkWdvvecl7dfjS8idGVesxMLtbA3+
TQhB5NuM39th+aXdjp8W2gQDNT5D2gUzLIdPpuWQyCawfArX/z2mh2UYGuqM
EtuEsXIW8ZP16lgqPjLh+5KbYG41EMCuoYQ9chIPZJHeBCPuy8XziCx2p9lg
dezSJlDYv05YzBbFrp4OVImTI623tRmi+IobO2NanmNwbRP2xP7wrWgwYkeK
F6i5sJswb3+RLvmAEtt9jMH6B2xCkHjLW96UXwpNWtivSUqboFh9qxwvsq5Q
kWojYnJzE+zKCsUPn08r5G9/DOJVIcXXPzGgy6xfIU2xfW3h9ibIWT2dPzbY
qBAbsaOaob4JB9MRZW2aRQrv57nzHmtuwtsrE1J3TVMUfC+qnxC6twndKYPS
d5nCFV55u9us3t8EpXMzWgzfPRScv2V35+pvgnE3RavLqK2CFe930acPN+Ez
wbw15paBgk6z5AbBZBPUFIOf+Zy5pCB2bOnp7tNNoNRN1Wqnn5Of9tKw++e9
CfcqVqzRLir5ocHXfY2+m1BYUPFZSnn5WhdPnrh3AOn8vdJ2Kx+3XStroiRS
hJL2L1qyeOzDq2shZDX2NLGb8DpOb+OV4NBVH82V/u54kj402p2ub6VedUth
kgxJJK3/4XzXZQ+HqxaowyZd+iaoclJSzd4hu2oQnqAxmLkJz55Q8s03tlzR
nPtaGJGzCYQyRwNqrM8VZalf9Pc+bwLvj5bCtz7YK1gvfkemgk3ghHctM693
5f4HLXZaiA==
"]]},
Annotation[#, "Charting`Private`Tag$663836#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->NCache[{{
Rational[1, 100],
Rational[99, 100]}, {0., 2.6592757861488967`}}, {{0.01, 0.99}, {0.,
2.6592757861488967`}}],
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.05],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}], TraditionalForm]], "Output",
CellChangeTimes->{3.8236217986755753`*^9, 3.823621962607039*^9},
CellLabel->
"Out[170]=",ExpressionUUID->"4ae2f1dd-2150-4518-8f3f-9e853a8f0a94"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"monitoredFindRoot", "[", "args__", "]"}], " ", ":=", " ",
RowBox[{"Module", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"s", " ", "=", " ", "0"}], ",", " ",
RowBox[{"e", " ", "=", " ", "0"}], ",",
RowBox[{"j", " ", "=", " ", "0"}]}], "}"}], ",", "\[IndentingNewLine]",
RowBox[{"{",
RowBox[{
RowBox[{"FindRoot", "[",
RowBox[{"args", ",", " ",
RowBox[{"StepMonitor", "\[RuleDelayed]",
RowBox[{"s", "++"}]}], ",",
RowBox[{"EvaluationMonitor", "\[RuleDelayed]",
RowBox[{"e", "++"}]}], ",",
RowBox[{"Jacobian", "\[Rule]",
RowBox[{"{",
RowBox[{"Automatic", ",",
RowBox[{"EvaluationMonitor", "\[RuleDelayed]",
RowBox[{"j", "++"}]}]}], "}"}]}]}], "]"}], ",",
RowBox[{"\"\<Steps\>\"", "\[Rule]", "s"}], ",",
RowBox[{"\"\<Evaluations\>\"", "\[Rule]", "e"}], ",",
RowBox[{"\"\<Jacobian Evaluations\>\"", "\[Rule]", "j"}]}], "}"}]}],
"]"}]}]], "Input",
CellLabel->
"In[200]:=",ExpressionUUID->"937865e9-068d-42b0-96ab-ad7e06b99a23"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"leftBinTest", "=", "0.6"}], ";"}], "\[IndentingNewLine]",
RowBox[{"monitoredFindRoot", "[",
RowBox[{
RowBox[{
RowBox[{"Integrate", "[",
RowBox[{
RowBox[{
RowBox[{"HF", "[",
FractionBox["x",
RowBox[{"1", "-", "x"}]], "]"}],
FractionBox["1",
SuperscriptBox[
RowBox[{"(",
RowBox[{"1", "-", "x"}], ")"}], "2"]]}], ",",
RowBox[{"{",
RowBox[{"x", ",", "leftBinTest", ",",
RowBox[{"leftBinTest", "+", "b"}]}], "}"}]}], "]"}], "-",
FractionBox["1", "10"]}], ",",
RowBox[{"{",
RowBox[{"b", ",", "0.1", ",", "0.", ",", "0.5"}], "}"}], ",",
RowBox[{"PrecisionGoal", "\[Rule]", "16"}]}],
RowBox[{"(*",
RowBox[{",",
RowBox[{"Method", "\[Rule]", "\"\<Brent\>\""}]}], "*)"}], "]"}]}], "Input",\
CellChangeTimes->{{3.823622544465562*^9, 3.823622686343278*^9}, {
3.8236227487298098`*^9, 3.823622767274707*^9}, {3.8236228296678047`*^9,
3.823622877530225*^9}, {3.823622911964579*^9, 3.823622914004518*^9}, {
3.8236229575386972`*^9, 3.823622988570241*^9}, {3.823623043955024*^9,
3.823623046082788*^9}, {3.823623079353895*^9, 3.8236231035652733`*^9}, {
3.823623231878303*^9, 3.8236232684739237`*^9}, {3.823623574124254*^9,
3.823623699535852*^9}, {3.823623738851583*^9, 3.8236238846263103`*^9}, {
3.823623929536949*^9, 3.823623939775478*^9}, {3.823623982381501*^9,
3.823623982545405*^9}, 3.8236248458259687`*^9},
CellLabel->
"In[270]:=",ExpressionUUID->"7e8ce4e4-f2af-4b2b-8a57-4a1b35dd787e"],
Cell[BoxData[
FormBox[
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"b", "\[Rule]", "0.03821517455754332`"}], "}"}], ",",
RowBox[{"\<\"Steps\"\>", "\[Rule]", "4"}], ",",
RowBox[{"\<\"Evaluations\"\>", "\[Rule]", "5"}], ",",
RowBox[{"\<\"Jacobian Evaluations\"\>", "\[Rule]", "4"}]}], "}"}],
TraditionalForm]], "Output",
CellChangeTimes->{{3.823622567500145*^9, 3.823622598189794*^9}, {
3.8236226365098753`*^9, 3.8236226607472973`*^9}, 3.823622691901698*^9, {
3.823622749537052*^9, 3.8236227800162497`*^9}, {3.823622834926324*^9,
3.8236229264092216`*^9}, 3.8236229734185257`*^9, 3.8236230056324*^9,
3.823623062995899*^9, {3.8236230952644587`*^9, 3.823623117705291*^9},
3.823623351093824*^9, 3.8236235866235847`*^9, {3.823623684289515*^9,
3.8236237116501904`*^9}, {3.823623801124671*^9, 3.823623873280802*^9}, {
3.823623932944508*^9, 3.823623945672402*^9}, 3.8236239923610477`*^9,
3.8236248616749268`*^9},
CellLabel->
"Out[271]=",ExpressionUUID->"2fd3aee9-60e9-4f73-9f51-c57a127195a6"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"leftBinTest", "=", "0.0"}], ";"}], "\[IndentingNewLine]",
RowBox[{"monitoredFindRoot", "[",
RowBox[{
RowBox[{
RowBox[{"Integrate", "[",
RowBox[{
RowBox[{
RowBox[{"HF", "[",
FractionBox["x",
RowBox[{"1", "-", "x"}]], "]"}],
FractionBox["1",
SuperscriptBox[
RowBox[{"(",
RowBox[{"1", "-", "x"}], ")"}], "2"]]}], ",",
RowBox[{"{",
RowBox[{"x", ",", "leftBinTest", ",",
RowBox[{"leftBinTest", "+", "b"}]}], "}"}]}], "]"}], "-",
FractionBox["1", "10"]}], ",",
RowBox[{"{",
RowBox[{"b", ",", "0.1", ",", "0.", ",", "0.5"}], "}"}], ",",
RowBox[{"PrecisionGoal", "\[Rule]", "16"}]}],
RowBox[{"(*",
RowBox[{",",
RowBox[{"Method", "\[Rule]", "\"\<Brent\>\""}]}], "*)"}], "]"}]}], "Input",\
CellChangeTimes->{{3.823622544465562*^9, 3.823622686343278*^9}, {
3.8236227487298098`*^9, 3.823622767274707*^9}, {3.8236228296678047`*^9,
3.823622877530225*^9}, {3.823622911964579*^9, 3.823622914004518*^9}, {
3.8236229575386972`*^9, 3.823622988570241*^9}, {3.823623043955024*^9,
3.823623046082788*^9}, {3.823623079353895*^9, 3.8236231035652733`*^9}, {
3.823623231878303*^9, 3.8236232684739237`*^9}, {3.823623574124254*^9,
3.823623699535852*^9}, {3.823623738851583*^9, 3.8236238846263103`*^9}, {
3.823623929536949*^9, 3.823623939775478*^9}, {3.823623982381501*^9,
3.823623982545405*^9}, {3.8236242878995733`*^9, 3.823624288015156*^9}, {
3.823624372786334*^9, 3.823624426749715*^9}, 3.823624480639462*^9, {
3.8236246009393387`*^9, 3.8236246097564707`*^9}, 3.823624695575116*^9, {
3.8236247508347282`*^9, 3.823624778121931*^9}},
CellLabel->
"In[264]:=",ExpressionUUID->"7dfc0a9c-7c85-4cb3-b6aa-3a51e64999cc"],
Cell[BoxData[
FormBox[
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"b", "\[Rule]", "0.39600932274980144`"}], "}"}], ",",
RowBox[{"\<\"Steps\"\>", "\[Rule]", "11"}], ",",
RowBox[{"\<\"Evaluations\"\>", "\[Rule]", "18"}], ",",
RowBox[{"\<\"Jacobian Evaluations\"\>", "\[Rule]", "11"}]}], "}"}],
TraditionalForm]], "Output",
CellChangeTimes->{{3.823624404252491*^9, 3.823624412502036*^9},
3.823624594302231*^9, 3.8236246683426847`*^9, {3.823624743783847*^9,
3.823624761136301*^9}, 3.823624828884489*^9},
CellLabel->
"Out[265]=",ExpressionUUID->"3817a3f4-453f-4da9-9ca6-44aec53755b6"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"PreviousBinBoundary", "=", "0."}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"NBins", "=", "10"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"BinLeftEdges", " ", "=", " ",
RowBox[{"{", "0.", "}"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"For", "[",
RowBox[{
RowBox[{"iBin", "=", "1"}], ",",
RowBox[{"iBin", "<", "NBins"}], ",",
RowBox[{"iBin", "++"}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"Print", "[",
RowBox[{"\"\<Computing right edge of bin #\>\"", "<>",
RowBox[{"ToString", "[", "iBin", "]"}], "<>", "\"\<...\>\""}], "]"}],
";", "\[IndentingNewLine]",
RowBox[{"AppendTo", "[", "\[IndentingNewLine]",
RowBox[{"BinLeftEdges", ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"BinLeftEdges", "\[LeftDoubleBracket]",
RowBox[{"-", "1"}], "\[RightDoubleBracket]"}], "+", "b"}], ")"}], "/.",
RowBox[{"(",
RowBox[{"FindRoot", "[",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"Integrate", "[",
RowBox[{
RowBox[{
RowBox[{"HF", "[",
FractionBox["x",
RowBox[{"1", "-", "x"}]], "]"}],
FractionBox["1",
SuperscriptBox[
RowBox[{"(",
RowBox[{"1", "-", "x"}], ")"}], "2"]]}], ",",
RowBox[{"{",
RowBox[{"x", ",",
RowBox[{"BinLeftEdges", "\[LeftDoubleBracket]",
RowBox[{"-", "1"}], "\[RightDoubleBracket]"}], ",",
RowBox[{
RowBox[{"BinLeftEdges", "\[LeftDoubleBracket]",
RowBox[{"-", "1"}], "\[RightDoubleBracket]"}], "+",
"b"}]}], "}"}]}], "]"}], "-",
RowBox[{"1.", "/", "NBins"}]}], ")"}], ",",
RowBox[{"{",
RowBox[{"b", ",", "0.1", ",", "0.", ",", "0.5"}], "}"}], ",",
RowBox[{"PrecisionGoal", "\[Rule]", "16"}]}], "]"}], ")"}]}]}],
"\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]",
RowBox[{"Print", "[",
RowBox[{"\"\<...\>\"", "<>",
RowBox[{"ToString", "[",
RowBox[{"BinLeftEdges", "\[LeftDoubleBracket]",
RowBox[{"-", "1"}], "\[RightDoubleBracket]"}], "]"}]}], "]"}],
";"}]}], "\[IndentingNewLine]", "]"}], ";"}]}], "Input",
CellChangeTimes->{{3.823619426945154*^9, 3.823619618144533*^9}, {
3.823619652597979*^9, 3.823619710583579*^9}, {3.823619752998989*^9,
3.8236198217342978`*^9}, {3.8236218150380297`*^9,
3.8236218162300367`*^9}, {3.82362184956635*^9, 3.823621853555338*^9},
3.8236221830820827`*^9, {3.8236222882597427`*^9, 3.82362230272022*^9}, {
3.823622527333239*^9, 3.823622532513112*^9}, 3.823623138021388*^9,
3.8236240021683407`*^9, {3.8236241512222233`*^9, 3.82362415171529*^9}, {
3.823624280528253*^9, 3.823624281151475*^9}, 3.823624840541318*^9},
CellLabel->
"In[272]:=",ExpressionUUID->"f2898046-4f39-447e-9706-e324e4f4d7d5"],
Cell[CellGroupData[{
Cell[BoxData[
FormBox["\<\"Computing right edge of bin #1...\"\>",
TraditionalForm]], "Print",
CellChangeTimes->{3.8236221848513317`*^9, 3.823622303326539*^9,
3.8236225332077627`*^9, 3.823623143073491*^9, 3.823623415541893*^9,
3.8236240043568983`*^9, 3.823624152663488*^9, 3.823624843408736*^9,
3.823624886471148*^9},
CellLabel->
"During evaluation of \
In[272]:=",ExpressionUUID->"b21e4c43-cc6f-4b4a-8626-4a5002592326"],
Cell[BoxData[
FormBox["\<\"...0.396009\"\>", TraditionalForm]], "Print",
CellChangeTimes->{3.8236221848513317`*^9, 3.823622303326539*^9,
3.8236225332077627`*^9, 3.823623143073491*^9, 3.823623415541893*^9,
3.8236240043568983`*^9, 3.823624152663488*^9, 3.823624843408736*^9,
3.823624933355492*^9},
CellLabel->
"During evaluation of \
In[272]:=",ExpressionUUID->"4a8bd4ff-a0df-4e24-8b43-f3d54fd87ad8"],
Cell[BoxData[
FormBox["\<\"Computing right edge of bin #2...\"\>",
TraditionalForm]], "Print",
CellChangeTimes->{3.8236221848513317`*^9, 3.823622303326539*^9,
3.8236225332077627`*^9, 3.823623143073491*^9, 3.823623415541893*^9,
3.8236240043568983`*^9, 3.823624152663488*^9, 3.823624843408736*^9,
3.823624933357233*^9},
CellLabel->
"During evaluation of \
In[272]:=",ExpressionUUID->"ef25a292-e4de-41b4-b54b-96d40b5946e6"],
Cell[BoxData[
FormBox["\<\"...0.467879\"\>", TraditionalForm]], "Print",
CellChangeTimes->{3.8236221848513317`*^9, 3.823622303326539*^9,
3.8236225332077627`*^9, 3.823623143073491*^9, 3.823623415541893*^9,
3.8236240043568983`*^9, 3.823624152663488*^9, 3.823624843408736*^9,
3.8236249419540367`*^9},
CellLabel->
"During evaluation of \
In[272]:=",ExpressionUUID->"77d8183b-e570-4fc6-bb40-4cf99db21379"],
Cell[BoxData[
FormBox["\<\"Computing right edge of bin #3...\"\>",
TraditionalForm]], "Print",
CellChangeTimes->{3.8236221848513317`*^9, 3.823622303326539*^9,
3.8236225332077627`*^9, 3.823623143073491*^9, 3.823623415541893*^9,
3.8236240043568983`*^9, 3.823624152663488*^9, 3.823624843408736*^9,
3.823624941954699*^9},
CellLabel->
"During evaluation of \
In[272]:=",ExpressionUUID->"2dc4788d-2d6e-4019-9cf5-bfc265214f20"],
Cell[BoxData[
FormBox["\<\"...0.520801\"\>", TraditionalForm]], "Print",
CellChangeTimes->{3.8236221848513317`*^9, 3.823622303326539*^9,
3.8236225332077627`*^9, 3.823623143073491*^9, 3.823623415541893*^9,
3.8236240043568983`*^9, 3.823624152663488*^9, 3.823624843408736*^9,
3.823624951718741*^9},
CellLabel->
"During evaluation of \
In[272]:=",ExpressionUUID->"a35d75c0-7425-4ae6-83ca-a09b2c0b361e"],
Cell[BoxData[
FormBox["\<\"Computing right edge of bin #4...\"\>",
TraditionalForm]], "Print",
CellChangeTimes->{3.8236221848513317`*^9, 3.823622303326539*^9,
3.8236225332077627`*^9, 3.823623143073491*^9, 3.823623415541893*^9,
3.8236240043568983`*^9, 3.823624152663488*^9, 3.823624843408736*^9,
3.823624951719428*^9},
CellLabel->
"During evaluation of \
In[272]:=",ExpressionUUID->"3a23d504-ac4f-4dd7-bd13-9eac99290148"],
Cell[BoxData[
FormBox["\<\"...0.565322\"\>", TraditionalForm]], "Print",
CellChangeTimes->{3.8236221848513317`*^9, 3.823622303326539*^9,
3.8236225332077627`*^9, 3.823623143073491*^9, 3.823623415541893*^9,
3.8236240043568983`*^9, 3.823624152663488*^9, 3.823624843408736*^9,
3.82362496175991*^9},
CellLabel->
"During evaluation of \
In[272]:=",ExpressionUUID->"eea3a992-f7f7-4c5e-961f-ae605a1e63e8"],
Cell[BoxData[
FormBox["\<\"Computing right edge of bin #5...\"\>",
TraditionalForm]], "Print",
CellChangeTimes->{3.8236221848513317`*^9, 3.823622303326539*^9,
3.8236225332077627`*^9, 3.823623143073491*^9, 3.823623415541893*^9,
3.8236240043568983`*^9, 3.823624152663488*^9, 3.823624843408736*^9,
3.823624961760467*^9},
CellLabel->
"During evaluation of \
In[272]:=",ExpressionUUID->"be7f41d7-39cc-4033-86c8-c175b35b41fc"],
Cell[BoxData[
FormBox["\<\"...0.605451\"\>", TraditionalForm]], "Print",
CellChangeTimes->{3.8236221848513317`*^9, 3.823622303326539*^9,
3.8236225332077627`*^9, 3.823623143073491*^9, 3.823623415541893*^9,
3.8236240043568983`*^9, 3.823624152663488*^9, 3.823624843408736*^9,
3.8236249698533697`*^9},
CellLabel->
"During evaluation of \
In[272]:=",ExpressionUUID->"d0b9a781-c269-4b31-b0a7-0700f98d2272"],
Cell[BoxData[
FormBox["\<\"Computing right edge of bin #6...\"\>",
TraditionalForm]], "Print",
CellChangeTimes->{3.8236221848513317`*^9, 3.823622303326539*^9,
3.8236225332077627`*^9, 3.823623143073491*^9, 3.823623415541893*^9,
3.8236240043568983`*^9, 3.823624152663488*^9, 3.823624843408736*^9,
3.823624969853939*^9},
CellLabel->
"During evaluation of \
In[272]:=",ExpressionUUID->"2d9238a7-cd86-47dd-bcc9-541704eed0b4"],
Cell[BoxData[
FormBox["\<\"...0.643488\"\>", TraditionalForm]], "Print",
CellChangeTimes->{3.8236221848513317`*^9, 3.823622303326539*^9,
3.8236225332077627`*^9, 3.823623143073491*^9, 3.823623415541893*^9,
3.8236240043568983`*^9, 3.823624152663488*^9, 3.823624843408736*^9,
3.823624978982524*^9},
CellLabel->
"During evaluation of \
In[272]:=",ExpressionUUID->"1d7a5565-931d-455e-bc59-749324e4506c"],
Cell[BoxData[
FormBox["\<\"Computing right edge of bin #7...\"\>",
TraditionalForm]], "Print",