File: zeta.dem

package info (click to toggle)
gnuplot 6.0.2%2Bdfsg1-2
  • links: PTS, VCS
  • area: main
  • in suites: sid
  • size: 14,940 kB
  • sloc: ansic: 95,319; cpp: 7,590; makefile: 2,470; javascript: 2,328; sh: 1,531; lisp: 664; perl: 304; pascal: 191; tcl: 88; python: 46
file content (85 lines) | stat: -rw-r--r-- 2,143 bytes parent folder | download | duplicates (3)
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
#
# Riemann zeta fuction
#
# Check support
if (!strstrt(GPVAL_COMPILE_OPTIONS, "+COMPLEX_FUNC")) {
    print "This copy of gnuplot does not provide zeta(s)"
    exit  # return to caller
}
save_encoding = GPVAL_ENCODING
set encoding utf8

set title "Real and imaginary parts of ζ(s) along the critical line real(s) = 1/2"
set grid x y
set xzeroaxis
set sample 999
set key reverse left
set xrange [5:50]
set yrange [-3:4]

plot real(zeta( 0.5 + I*x )) lt 4 lw 3 title "real( ζ(½ + ix) )", \
     imag(zeta( 0.5 + I*x )) lt 6 lw 3 title "imag( ζ(½ + ix) )"

pause -1 "<cr> to continue"
reset

set title "{/Times=16 Riemann zeta function}    {/Times:Italic=16 ζ(s) = {∑@_{/*0.5 {k=1}}^{/*.5 ∞}}  {/:Italic=12 k^{ -s}}"

set samp 201
set iso 201

set pm3d depth base noborder
set pm3d interp 2,2
set pm3d lighting specular 0.1
set xyplane 0
unset colorbox
set key left
set border 895
set grid x y z vertical

set xlabel "Real(s)" offset -1,0 right
set ylabel "Imag(s)"

set xrange [-9 : 9]
set yrange [-9 : 9]
set zrange [0 : 4]

set palette defined (0 "dark-blue", 1 "white")
set zrange [-1.1 : 2.2]
splot real( zeta(x + I*y) ) with pm3d title "real(ζ(s))"
pause -1 "<cr> to continue"

set palette defined (0 "dark-red", 1 "white")
set zrange [-0.55 : 0.55]
splot imag( zeta(x + I*y) ) with pm3d title "imag(ζ(s))"
pause -1 "<cr> to continue"

set zrange [*:*]
unset grid
set view map
set key right opaque box samplen 0
set colorbox
set cblabel "Phase Angle"
set cblabel  offset -1, 0, 0  rotate
set cbrange [ -3.14159 : 3.14159 ] noreverse writeback
set cbtics   ("-π" -pi, "+π" pi)
set pm3d interp 1,1 nolighting corners2color c1
set palette color model HSV defined ( 0 0 1 1, 1 1 1 1 )

splot arg( zeta(x + I*y) ) with pm3d title "arg(ζ(s))"
pause -1 "<cr> to continue"

unset colorbox
set palette cubehelix negative
set pm3d lighting specular 0.1
set yrange [-10:50]
set xrange [-3:3]
set zrange [0:4]
set view 57, 100, 1.2, 0.5
set key nobox noopaque bottom left at screen 0.0, 0.6
set xlabel offset -3
splot 1./abs( zeta(x + I*y) ) with pm3d title "|ζ(s)|^{-1}"
pause -1 "<cr> to continue"

reset
set encoding save_encoding