Added reference comparison with GNAT and script for extracing percent differences. Also plot them.HEADmaster

5 files changed, 1041 insertions, 8 deletions

diff --git a/Makefile b/Makefile
index 78a410c..259d96a 100644
--- a/Makefile
+++ b/Makefile
@@ -6,13 +6,14 @@ XNAME = monte_carlo.x
 all: clean build run plot
 
 clean:
-	rm -f *.x *.out *.mod
+	rm -f *.x *.out *.mod *.tex *.pdf
 
 build: main.f90
 	$(FC) $(FFLAGS) main.f90 -o $(XNAME)
 
 run: build
 	./$(XNAME)
+	./extract_percent_diff.sh
 
 plot: run
 	gnuplot plot.gnu
diff --git a/extract_percent_diff.sh b/extract_percent_diff.sh
new file mode 100755
index 0000000..cd6e37e
--- /dev/null
+++ b/extract_percent_diff.sh
@@ -0,0 +1,11 @@
+#!/bin/bash
+
+# Get the percent difference between GNAT benchmark and MC calculation
+# Defined as (GNAT - MC)/GNAT * 100
+
+dx=$(cat flux.out | awk '{print $1}')
+mc=$(cat flux.out | awk '{print $2}')
+ref=$(cat test.ref | awk '{print $2}')
+
+paste <(echo "$dx") <(paste <(echo "$mc") <(echo "$ref") | awk '{print ($2 - $1)/$2*100}') > diff.out
+
diff --git a/main.f90 b/main.f90
index cd61bc5..2575d9d 100644
--- a/main.f90
+++ b/main.f90
@@ -4,8 +4,8 @@ program monte_carlo
     implicit none
 
     !> Initialize some variables for neutrons
-    integer(int64), parameter :: n = 1e7
-    integer(int64), parameter :: bins = 800
+    integer(int64), parameter :: n = 4e8
+    integer(int64), parameter :: bins = 1001
     real(real64), parameter :: eps = 1.0e-6_real64
     real(real64), parameter :: pi = 4.0_real64 * atan(1.0_real64)
 
diff --git a/plot.gnu b/plot.gnu
index e001571..80fd5cd 100644
--- a/plot.gnu
+++ b/plot.gnu
@@ -1,12 +1,32 @@
-#set term x11 persist
-set term cairolatex pdf size 4in,2.5in
+#set term x11 1 persist
+set term cairolatex pdf size 5in,3in
 set output "flux.tex"
 
-set key reverse top left
+set key outside top center horizontal
 set grid
 set xlabel "Distance [cm]"
 set ylabel "$\\phi$"
+set border lw 3
 
+# Draw the regions
+
+set lt 9 lw 1 lc rgb "black" dt 2
+set arrow from 2,graph 0 to 2, graph 1 nohead lt 9
+set arrow from 3,graph 0 to 3, graph 1 nohead lt 9
+set arrow from 5,graph 0 to 5, graph 1 nohead lt 9
+set arrow from 6,graph 0 to 6, graph 1 nohead lt 9
+
+
+plot 'flux.out' using 1:($2 + 3*$3):($2 - 3*$3) with filledcurves lc rgb "#CCCCFF" title "$\\pm 3\\sigma$", \
+    'flux.out' using 1:2 with lines lc rgb "red" lw 2 title "$\\overline{\\phi}_{MC}$", \
+    'test.ref' using 1:2 with lines lc rgb "green" lw 1 dt 2 title "$\\overline{\\phi}_{ref}$"
+
+
+#set term x11 2 persist
+set output "error.tex"
+set yrange [-15:15]
+set key off
+
+set ylabel "Percent Difference [\\%]"
+plot 'diff.out' using 1:2 with lines
 
-plot 'flux.out' using 1:($2 + $3):($2 - $3) with filledcurves lc rgb "#CCCCFF" title "$\\pm 1 \\sigma$", \
-    'flux.out' using 1:2 with lines lc rgb "red" title "$\\overline{\\phi}$"
diff --git a/test.ref b/test.ref
new file mode 100644
index 0000000..02b5cc7
--- /dev/null
+++ b/test.ref
@@ -0,0 +1,1001 @@
+0 1
+0.008 1
+0.016 1
+0.024 1
+0.032 1
+0.04 1
+0.048 1
+0.056 1
+0.064 1
+0.072 1
+0.08 1
+0.088 1
+0.096 1
+0.104 1
+0.112 1
+0.12 1
+0.128 1
+0.136 1
+0.144 1
+0.152 1
+0.16 1
+0.168 1
+0.176 1
+0.184 1
+0.192 1
+0.2 1
+0.208 1
+0.216 1
+0.224 1
+0.232 1
+0.24 1
+0.248 1
+0.256 1
+0.264 1
+0.272 1
+0.28 1
+0.288 1
+0.296 1
+0.304 1
+0.312 1
+0.32 1
+0.328 1
+0.336 1
+0.344 1
+0.352 1
+0.36 1
+0.368 1
+0.376 1
+0.384 1
+0.392 1
+0.4 1
+0.408 1
+0.416 1
+0.424 1
+0.432 1
+0.44 1
+0.448 1
+0.456 1
+0.464 1
+0.472 1
+0.48 1
+0.488 1
+0.496 1
+0.504 1
+0.512 1
+0.52 1
+0.528 1
+0.536 1
+0.544 1
+0.552 1
+0.56 1
+0.568 1
+0.576 1
+0.584 1
+0.592 1
+0.6 1
+0.608 1
+0.616 1
+0.624 1
+0.632 1
+0.64 1
+0.648 1
+0.656 1
+0.664 1
+0.672 1
+0.68 1
+0.688 1
+0.696 1
+0.704 1
+0.712 1
+0.72 1
+0.728 1
+0.736 1
+0.744 1
+0.752 1
+0.76 1
+0.768 1
+0.776 1
+0.784 1
+0.792 1
+0.8 1
+0.808 1
+0.816 1
+0.824 1
+0.832 1
+0.84 1
+0.848 1
+0.856 1
+0.864 1
+0.872 1
+0.88 1
+0.888 1
+0.896 1
+0.904 1
+0.912 1
+0.92 1
+0.928 1
+0.936 1
+0.944 1
+0.952 1
+0.96 1
+0.968 1
+0.976 1
+0.984 1
+0.992 1
+1 1
+1.008 1
+1.016 1
+1.024 1
+1.032 1
+1.04 1
+1.048 1
+1.056 1
+1.064 1
+1.072 1
+1.08 1
+1.088 1
+1.096 1
+1.104 1
+1.112 1
+1.12 1
+1.128 1
+1.136 1
+1.144 1
+1.152 1
+1.16 1
+1.168 1
+1.176 1
+1.184 1
+1.192 1
+1.2 1
+1.208 1
+1.216 1
+1.224 1
+1.232 1
+1.24 1
+1.248 1
+1.256 1
+1.264 1
+1.272 1
+1.28 1
+1.288 1
+1.296 1
+1.304 1
+1.312 1
+1.32 1
+1.328 1
+1.336 1
+1.344 1
+1.352 1
+1.36 1
+1.368 1
+1.376 1
+1.384 1
+1.392 1
+1.4 1
+1.408 1
+1.416 1
+1.424 1
+1.432 1
+1.44 1
+1.448 1
+1.456 1
+1.464 1
+1.472 1
+1.48 1
+1.488 1
+1.496 1
+1.504 1
+1.512 1
+1.52 1
+1.528 1
+1.536 1
+1.544 1
+1.552 1
+1.56 1
+1.568 1
+1.576 1
+1.584 1
+1.592 1
+1.6 1
+1.608 1
+1.616 1
+1.624 1
+1.632 1
+1.64 1
+1.648 1
+1.656 1
+1.664 1
+1.672 1
+1.68 1
+1.688 1
+1.696 1
+1.704 1
+1.712 1
+1.72 1
+1.728 1
+1.736 1
+1.744 1
+1.752 1
+1.76 1
+1.768 1
+1.776 1
+1.784 1
+1.792 1
+1.8 1
+1.808 1
+1.816 1
+1.824 0.99999
+1.832 0.99999
+1.84 0.99998
+1.848 0.99997
+1.856 0.99996
+1.864 0.99994
+1.872 0.9999
+1.88 0.99985
+1.888 0.99976
+1.896 0.99962
+1.904 0.9994
+1.912 0.99904
+1.92 0.99846
+1.928 0.99752
+1.936 0.99598
+1.944 0.99343
+1.952 0.98916
+1.96 0.98189
+1.968 0.96919
+1.976 0.94668
+1.984 0.90292
+1.992 0.8212
+2 0.55201
+2.008 0.41275
+2.016 0.38925
+2.024 0.34222
+2.032 0.31611
+2.04 0.28778
+2.048 0.26586
+2.056 0.24515
+2.064 0.22731
+2.072 0.21096
+2.08 0.19634
+2.088 0.18298
+2.096 0.17083
+2.104 0.15971
+2.112 0.14951
+2.12 0.14013
+2.128 0.13149
+2.136 0.12352
+2.144 0.11614
+2.152 0.10932
+2.16 0.10299
+2.168 0.097114
+2.176 0.091658
+2.184 0.086586
+2.192 0.081867
+2.2 0.077475
+2.208 0.073385
+2.216 0.069575
+2.224 0.066026
+2.232 0.062718
+2.24 0.059635
+2.248 0.056764
+2.256 0.054089
+2.264 0.051599
+2.272 0.049282
+2.28 0.047128
+2.288 0.045128
+2.296 0.043272
+2.304 0.041554
+2.312 0.039965
+2.32 0.038499
+2.328 0.037151
+2.336 0.035915
+2.344 0.034787
+2.352 0.033761
+2.36 0.032834
+2.368 0.032001
+2.376 0.031259
+2.384 0.030605
+2.392 0.030037
+2.4 0.029552
+2.408 0.029149
+2.416 0.028826
+2.424 0.028582
+2.432 0.028415
+2.44 0.028325
+2.448 0.028312
+2.456 0.028375
+2.464 0.028513
+2.472 0.028729
+2.48 0.029022
+2.488 0.029392
+2.496 0.029842
+2.504 0.030372
+2.512 0.030984
+2.52 0.03168
+2.528 0.032463
+2.536 0.033334
+2.544 0.034297
+2.552 0.035355
+2.56 0.036511
+2.568 0.037769
+2.576 0.039134
+2.584 0.040611
+2.592 0.042203
+2.6 0.043918
+2.608 0.045761
+2.616 0.047738
+2.624 0.049857
+2.632 0.052125
+2.64 0.054551
+2.648 0.057144
+2.656 0.059913
+2.664 0.06287
+2.672 0.066028
+2.68 0.0694
+2.688 0.072997
+2.696 0.076833
+2.704 0.080925
+2.712 0.08529
+2.72 0.089947
+2.728 0.094915
+2.736 0.10022
+2.744 0.10588
+2.752 0.11192
+2.76 0.11838
+2.768 0.12528
+2.776 0.13266
+2.784 0.14056
+2.792 0.14902
+2.8 0.15807
+2.808 0.16779
+2.816 0.17821
+2.824 0.18941
+2.832 0.20145
+2.84 0.21441
+2.848 0.22838
+2.856 0.24346
+2.864 0.25975
+2.872 0.27739
+2.88 0.29652
+2.888 0.31731
+2.896 0.33995
+2.904 0.36467
+2.912 0.39175
+2.92 0.4215
+2.928 0.45433
+2.936 0.49071
+2.944 0.53125
+2.952 0.57676
+2.96 0.62829
+2.968 0.68734
+2.976 0.75612
+2.984 0.83841
+2.992 0.94037
+3 1.0926
+3.008 1.1167
+3.016 1.0978
+3.024 1.1135
+3.032 1.1004
+3.04 1.1114
+3.048 1.1021
+3.056 1.11
+3.064 1.1033
+3.072 1.109
+3.08 1.1042
+3.088 1.1082
+3.096 1.1048
+3.104 1.1076
+3.112 1.1053
+3.12 1.1072
+3.128 1.1057
+3.136 1.1069
+3.144 1.106
+3.152 1.1066
+3.16 1.1062
+3.168 1.1064
+3.176 1.1064
+3.184 1.1063
+3.192 1.1065
+3.2 1.1062
+3.208 1.1066
+3.216 1.1061
+3.224 1.1067
+3.232 1.1061
+3.24 1.1067
+3.248 1.1061
+3.256 1.1067
+3.264 1.106
+3.272 1.1067
+3.28 1.106
+3.288 1.1067
+3.296 1.106
+3.304 1.1067
+3.312 1.1061
+3.32 1.1067
+3.328 1.1061
+3.336 1.1067
+3.344 1.1061
+3.352 1.1066
+3.36 1.1062
+3.368 1.1065
+3.376 1.1062
+3.384 1.1065
+3.392 1.1063
+3.4 1.1065
+3.408 1.1063
+3.416 1.1064
+3.424 1.1063
+3.432 1.1064
+3.44 1.1063
+3.448 1.1064
+3.456 1.1064
+3.464 1.1064
+3.472 1.1064
+3.48 1.1064
+3.488 1.1064
+3.496 1.1064
+3.504 1.1064
+3.512 1.1064
+3.52 1.1064
+3.528 1.1064
+3.536 1.1064
+3.544 1.1064
+3.552 1.1064
+3.56 1.1064
+3.568 1.1064
+3.576 1.1064
+3.584 1.1064
+3.592 1.1064
+3.6 1.1064
+3.608 1.1064
+3.616 1.1064
+3.624 1.1064
+3.632 1.1064
+3.64 1.1064
+3.648 1.1064
+3.656 1.1064
+3.664 1.1064
+3.672 1.1064
+3.68 1.1064
+3.688 1.1064
+3.696 1.1064
+3.704 1.1064
+3.712 1.1064
+3.72 1.1064
+3.728 1.1064
+3.736 1.1064
+3.744 1.1064
+3.752 1.1064
+3.76 1.1064
+3.768 1.1064
+3.776 1.1064
+3.784 1.1064
+3.792 1.1064
+3.8 1.1064
+3.808 1.1064
+3.816 1.1064
+3.824 1.1064
+3.832 1.1064
+3.84 1.1064
+3.848 1.1064
+3.856 1.1064
+3.864 1.1064
+3.872 1.1064
+3.88 1.1064
+3.888 1.1064
+3.896 1.1064
+3.904 1.1064
+3.912 1.1064
+3.92 1.1064
+3.928 1.1064
+3.936 1.1064
+3.944 1.1064
+3.952 1.1064
+3.96 1.1064
+3.968 1.1064
+3.976 1.1064
+3.984 1.1064
+3.992 1.1064
+4 1.1064
+4.008 1.1064
+4.016 1.1064
+4.024 1.1064
+4.032 1.1064
+4.04 1.1064
+4.048 1.1064
+4.056 1.1064
+4.064 1.1064
+4.072 1.1064
+4.08 1.1064
+4.088 1.1064
+4.096 1.1064
+4.104 1.1064
+4.112 1.1064
+4.12 1.1064
+4.128 1.1064
+4.136 1.1064
+4.144 1.1064
+4.152 1.1064
+4.16 1.1064
+4.168 1.1064
+4.176 1.1064
+4.184 1.1064
+4.192 1.1064
+4.2 1.1064
+4.208 1.1064
+4.216 1.1064
+4.224 1.1064
+4.232 1.1064
+4.24 1.1064
+4.248 1.1064
+4.256 1.1064
+4.264 1.1064
+4.272 1.1064
+4.28 1.1064
+4.288 1.1064
+4.296 1.1064
+4.304 1.1064
+4.312 1.1064
+4.32 1.1064
+4.328 1.1064
+4.336 1.1064
+4.344 1.1064
+4.352 1.1064
+4.36 1.1064
+4.368 1.1064
+4.376 1.1064
+4.384 1.1064
+4.392 1.1064
+4.4 1.1064
+4.408 1.1064
+4.416 1.1064
+4.424 1.1064
+4.432 1.1064
+4.44 1.1064
+4.448 1.1064
+4.456 1.1064
+4.464 1.1064
+4.472 1.1064
+4.48 1.1064
+4.488 1.1064
+4.496 1.1064
+4.504 1.1064
+4.512 1.1064
+4.52 1.1064
+4.528 1.1064
+4.536 1.1064
+4.544 1.1064
+4.552 1.1064
+4.56 1.1064
+4.568 1.1064
+4.576 1.1064
+4.584 1.1064
+4.592 1.1064
+4.6 1.1064
+4.608 1.1064
+4.616 1.1064
+4.624 1.1064
+4.632 1.1064
+4.64 1.1064
+4.648 1.1064
+4.656 1.1064
+4.664 1.1063
+4.672 1.1064
+4.68 1.1063
+4.688 1.1064
+4.696 1.1063
+4.704 1.1064
+4.712 1.1063
+4.72 1.1064
+4.728 1.1063
+4.736 1.1064
+4.744 1.1063
+4.752 1.1064
+4.76 1.1063
+4.768 1.1064
+4.776 1.1063
+4.784 1.1064
+4.792 1.1064
+4.8 1.1064
+4.808 1.1064
+4.816 1.1064
+4.824 1.1064
+4.832 1.1064
+4.84 1.1064
+4.848 1.1064
+4.856 1.1064
+4.864 1.1063
+4.872 1.1064
+4.88 1.1063
+4.888 1.1065
+4.896 1.1062
+4.904 1.1065
+4.912 1.1062
+4.92 1.1066
+4.928 1.1061
+4.936 1.1067
+4.944 1.106
+4.952 1.1068
+4.96 1.1059
+4.968 1.107
+4.976 1.1057
+4.984 1.1072
+4.992 1.1055
+5 1.1075
+5.008 1.1518
+5.016 1.188
+5.024 1.2198
+5.032 1.2487
+5.04 1.2756
+5.048 1.301
+5.056 1.3251
+5.064 1.3481
+5.072 1.3702
+5.08 1.3914
+5.088 1.4119
+5.096 1.4317
+5.104 1.4508
+5.112 1.4693
+5.12 1.4872
+5.128 1.5046
+5.136 1.5215
+5.144 1.5379
+5.152 1.5538
+5.16 1.5693
+5.168 1.5844
+5.176 1.5991
+5.184 1.6133
+5.192 1.6272
+5.2 1.6407
+5.208 1.6539
+5.216 1.6667
+5.224 1.6792
+5.232 1.6914
+5.24 1.7032
+5.248 1.7147
+5.256 1.7259
+5.264 1.7368
+5.272 1.7474
+5.28 1.7578
+5.288 1.7678
+5.296 1.7776
+5.304 1.787
+5.312 1.7962
+5.32 1.8052
+5.328 1.8139
+5.336 1.8223
+5.344 1.8304
+5.352 1.8383
+5.36 1.846
+5.368 1.8534
+5.376 1.8605
+5.384 1.8674
+5.392 1.8741
+5.4 1.8805
+5.408 1.8867
+5.416 1.8926
+5.424 1.8983
+5.432 1.9038
+5.44 1.909
+5.448 1.914
+5.456 1.9188
+5.464 1.9234
+5.472 1.9277
+5.48 1.9318
+5.488 1.9357
+5.496 1.9393
+5.504 1.9427
+5.512 1.9459
+5.52 1.9489
+5.528 1.9516
+5.536 1.9542
+5.544 1.9565
+5.552 1.9585
+5.56 1.9604
+5.568 1.962
+5.576 1.9634
+5.584 1.9646
+5.592 1.9656
+5.6 1.9663
+5.608 1.9668
+5.616 1.9671
+5.624 1.9672
+5.632 1.967
+5.64 1.9666
+5.648 1.966
+5.656 1.9652
+5.664 1.9641
+5.672 1.9628
+5.68 1.9612
+5.688 1.9595
+5.696 1.9574
+5.704 1.9552
+5.712 1.9527
+5.72 1.9499
+5.728 1.9469
+5.736 1.9437
+5.744 1.9402
+5.752 1.9365
+5.76 1.9325
+5.768 1.9282
+5.776 1.9237
+5.784 1.9189
+5.792 1.9139
+5.8 1.9085
+5.808 1.9029
+5.816 1.897
+5.824 1.8907
+5.832 1.8842
+5.84 1.8774
+5.848 1.8702
+5.856 1.8628
+5.864 1.8549
+5.872 1.8467
+5.88 1.8382
+5.888 1.8292
+5.896 1.8199
+5.904 1.8101
+5.912 1.7999
+5.92 1.7892
+5.928 1.778
+5.936 1.7662
+5.944 1.7539
+5.952 1.7407
+5.96 1.727
+5.968 1.7121
+5.976 1.6965
+5.984 1.6785
+5.992 1.6591
+6 1.6339
+6.008 1.6087
+6.016 1.5891
+6.024 1.5708
+6.032 1.5546
+6.04 1.539
+6.048 1.5245
+6.056 1.5104
+6.064 1.497
+6.072 1.484
+6.08 1.4714
+6.088 1.4592
+6.096 1.4473
+6.104 1.4357
+6.112 1.4244
+6.12 1.4133
+6.128 1.4025
+6.136 1.3919
+6.144 1.3816
+6.152 1.3714
+6.16 1.3614
+6.168 1.3516
+6.176 1.3419
+6.184 1.3324
+6.192 1.3231
+6.2 1.3138
+6.208 1.3048
+6.216 1.2958
+6.224 1.287
+6.232 1.2783
+6.24 1.2697
+6.248 1.2612
+6.256 1.2528
+6.264 1.2446
+6.272 1.2364
+6.28 1.2283
+6.288 1.2203
+6.296 1.2124
+6.304 1.2046
+6.312 1.1969
+6.32 1.1892
+6.328 1.1817
+6.336 1.1742
+6.344 1.1668
+6.352 1.1594
+6.36 1.1521
+6.368 1.1449
+6.376 1.1378
+6.384 1.1307
+6.392 1.1237
+6.4 1.1168
+6.408 1.1099
+6.416 1.103
+6.424 1.0963
+6.432 1.0895
+6.44 1.0829
+6.448 1.0763
+6.456 1.0697
+6.464 1.0632
+6.472 1.0568
+6.48 1.0504
+6.488 1.044
+6.496 1.0377
+6.504 1.0314
+6.512 1.0252
+6.52 1.019
+6.528 1.0129
+6.536 1.0068
+6.544 1.0008
+6.552 0.99477
+6.56 0.98881
+6.568 0.98289
+6.576 0.97701
+6.584 0.97116
+6.592 0.96536
+6.6 0.95959
+6.608 0.95386
+6.616 0.94817
+6.624 0.94251
+6.632 0.93688
+6.64 0.93129
+6.648 0.92574
+6.656 0.92021
+6.664 0.91472
+6.672 0.90927
+6.68 0.90384
+6.688 0.89845
+6.696 0.89309
+6.704 0.88776
+6.712 0.88245
+6.72 0.87718
+6.728 0.87194
+6.736 0.86672
+6.744 0.86154
+6.752 0.85638
+6.76 0.85125
+6.768 0.84615
+6.776 0.84107
+6.784 0.83602
+6.792 0.831
+6.8 0.826
+6.808 0.82103
+6.816 0.81608
+6.824 0.81116
+6.832 0.80626
+6.84 0.80138
+6.848 0.79653
+6.856 0.79171
+6.864 0.7869
+6.872 0.78212
+6.88 0.77736
+6.888 0.77263
+6.896 0.76791
+6.904 0.76322
+6.912 0.75855
+6.92 0.7539
+6.928 0.74927
+6.936 0.74467
+6.944 0.74008
+6.952 0.73551
+6.96 0.73097
+6.968 0.72644
+6.976 0.72193
+6.984 0.71744
+6.992 0.71297
+7 0.70852
+7.008 0.70409
+7.016 0.69968
+7.024 0.69528
+7.032 0.69091
+7.04 0.68655
+7.048 0.6822
+7.056 0.67788
+7.064 0.67357
+7.072 0.66928
+7.08 0.66501
+7.088 0.66075
+7.096 0.65651
+7.104 0.65228
+7.112 0.64807
+7.12 0.64388
+7.128 0.6397
+7.136 0.63554
+7.144 0.63139
+7.152 0.62726
+7.16 0.62314
+7.168 0.61903
+7.176 0.61494
+7.184 0.61087
+7.192 0.60681
+7.2 0.60276
+7.208 0.59872
+7.216 0.5947
+7.224 0.5907
+7.232 0.5867
+7.24 0.58272
+7.248 0.57875
+7.256 0.5748
+7.264 0.57085
+7.272 0.56692
+7.28 0.563
+7.288 0.55909
+7.296 0.5552
+7.304 0.55131
+7.312 0.54744
+7.32 0.54358
+7.328 0.53972
+7.336 0.53588
+7.344 0.53206
+7.352 0.52824
+7.36 0.52443
+7.368 0.52063
+7.376 0.51684
+7.384 0.51306
+7.392 0.50929
+7.4 0.50553
+7.408 0.50178
+7.416 0.49804
+7.424 0.49431
+7.432 0.49058
+7.44 0.48687
+7.448 0.48316
+7.456 0.47946
+7.464 0.47577
+7.472 0.47208
+7.48 0.46841
+7.488 0.46474
+7.496 0.46108
+7.504 0.45742
+7.512 0.45377
+7.52 0.45013
+7.528 0.4465
+7.536 0.44287
+7.544 0.43924
+7.552 0.43562
+7.56 0.43201
+7.568 0.4284
+7.576 0.4248
+7.584 0.4212
+7.592 0.41761
+7.6 0.41402
+7.608 0.41043
+7.616 0.40685
+7.624 0.40327
+7.632 0.39969
+7.64 0.39612
+7.648 0.39255
+7.656 0.38898
+7.664 0.38541
+7.672 0.38184
+7.68 0.37828
+7.688 0.37471
+7.696 0.37114
+7.704 0.36758
+7.712 0.36401
+7.72 0.36044
+7.728 0.35687
+7.736 0.3533
+7.744 0.34972
+7.752 0.34614
+7.76 0.34256
+7.768 0.33897
+7.776 0.33537
+7.784 0.33177
+7.792 0.32816
+7.8 0.32455
+7.808 0.32092
+7.816 0.31728
+7.824 0.31364
+7.832 0.30998
+7.84 0.3063
+7.848 0.30261
+7.856 0.2989
+7.864 0.29518
+7.872 0.29143
+7.88 0.28766
+7.888 0.28386
+7.896 0.28003
+7.904 0.27617
+7.912 0.27227
+7.92 0.26832
+7.928 0.26433
+7.936 0.26029
+7.944 0.25617
+7.952 0.25198
+7.96 0.24769
+7.968 0.24328
+7.976 0.23871
+7.984 0.23387
+7.992 0.22863
+8.000 0.22261