-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathplot.gp
54 lines (45 loc) · 2 KB
/
plot.gp
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
set terminal qt font "Open Sans" size 1100,600
set multiplot layout 2,3 rowsfirst
#set xrange [-1:0.95]
set yrange [0:10]
plot "plot.dat" using 1:2 with lines title "1/acos", \
"" using 1:3 with lines title "1/jd marlin", \
"" using 1:5 with lines title "1/approx", \
"" using 1:7 with lines title "acosInv", \
"" using 1:9 with lines title "acosInvQ"
set yrange [-0.1:0.1]
plot 0, \
"plot.dat" using 1:4 with lines title "1/acos(x) - 1/jd marlin", \
"" using 1:6 with lines title "1/acos(x) - 1/approx",\
"" using 1:8 with lines title "1/acos(x) - acosInv", \
"" using 1:10 with lines title "1/acos(x) - acosInvQ"
unset yrange
set xrange [0.995:]
set yrange [10:90]
plot "plot.dat" using 1:2 with lines title "1/acos", \
"" using 1:3 with lines title "1/jd marlin", \
"" using 1:5 with lines title "1/approx", \
"" using 1:7 with lines title "acosInv", \
"" using 1:9 with lines title "acosInvQ"
set xrange [0.998:]
set yrange [-5:5]
plot 0, \
"plot.dat" using 1:4 with lines title "1/acos(x) - 1/jd marlin", \
"" using 1:6 with lines title "1/acos(x) - 1/approx",\
"" using 1:8 with lines title "1/acos(x) - acosInv", \
"" using 1:10 with lines title "1/acos(x) - acosInvQ"
unset yrange
set xrange [0.9996:]
set yrange [10:90]
plot "plot.dat" using 1:2 with lines title "1/acos", \
"" using 1:3 with lines title "1/jd marlin", \
"" using 1:5 with lines title "1/approx", \
"" using 1:7 with lines title "acosInv", \
"" using 1:9 with lines title "acosInvQ"
set yrange [-5:5]
plot 0, \
"plot.dat" using 1:4 with lines title "1/acos(x) - 1/jd marlin", \
"" using 1:6 with lines title "1/acos(x) - 1/approx",\
"" using 1:8 with lines title "1/acos(x) - acosInv", \
"" using 1:10 with lines title "1/acos(x) - acosInvQ"
unset multiplot