#script to draw the mode structure for a chain of three coupled masses

set terminal postscript eps color enhanced size 14cm,12cm \
font "Times-Italic,24pt" dashlength 2.0
set output 'modestructure.eps'

set multiplot
set border 0
set xtics (1, 2, 3) axis offset 0,.3
unset ytics
unset xlabel
unset ylabel
set arrow 1 from 0,0 to 4.5 filled
set arrow 2 from 0,-1.2 to 0,1.3 filled
set label 1 "a_j" at -.35,1
set label 2 "j" at 4.2,-.3
set xrange [0:4]
set yrange [-1:1]
pcol='blue'

# plot w r=1
set origin 0,0.75
set size 0.5,0.24

set label 3 "r=1" at .3,1.2

plot sin(pi*x/4) lt 2 lw 0.5 lc rgb 'red' notitle, \
     '-' using ($1):(sin(pi*$1/4)) pt 7 ps 1.2 lc rgb pcol notitle
1
2
3
e

# plot w r=2
set origin 0,0.5
r=2
set label 3 "r=2" at .3,1.2

plot sin(r*pi*x/4) lt 2 lw 0.5 lc rgb 'red' notitle, \
     '-' using ($1):(sin(r*pi*$1/4)) pt 7 ps 1.2 lc rgb pcol notitle
1
2
3
e

# plot w r=3
set origin 0,0.25
r=3
set label 3 "r=3" at .3,1.2

plot sin(r*pi*x/4) lt 2 lw 0.5 lc rgb 'red' notitle, \
     '-' using ($1):(sin(r*pi*$1/4)) pt 7 ps 1.2 lc rgb pcol notitle
1
2
3
e

# plot w r=4
set origin 0,0
r=4
set label 3 "r=4" at .3,1.2

plot sin(r*pi*x/4) lt 2 lw 0.5 lc rgb 'red' notitle, \
     '-' using ($1):(sin(r*pi*$1/4)) pt 7 ps 1.2 lc rgb pcol notitle
1
2
3
e

# plot w r=5
set origin 0.5,0.75
r=5
set label 3 "r=5" at .3,1.2

plot sin(r*pi*x/4) lt 2 lw 0.5 lc rgb 'red' notitle, \
     '-' using ($1):(sin(r*pi*$1/4)) pt 7 ps 1.2 lc rgb pcol notitle
1
2
3
e

# plot w r=6
set origin 0.5,0.5
r=6
set label 3 "r=6" at .3,1.2

plot sin(r*pi*x/4) lt 2 lw 0.5 lc rgb 'red' notitle, \
     '-' using ($1):(sin(r*pi*$1/4)) pt 7 ps 1.2 lc rgb pcol notitle
1
2
3
e

# plot w r=7
set origin 0.5,0.25
r=7
set label 3 "r=7" at .3,1.2

plot sin(r*pi*x/4) lt 2 lw 0.5 lc rgb 'red' notitle, \
     '-' using ($1):(sin(r*pi*$1/4)) pt 7 ps 1.2 lc rgb pcol notitle
1
2
3
e

# plot w r=8
set origin 0.5,0
r=8
set label 3 "r=8" at .3,1.2

plot sin(r*pi*x/4) lt 2 lw 0.5 lc rgb 'red' notitle, \
     '-' using ($1):(sin(r*pi*$1/4)) pt 7 ps 1.2 lc rgb pcol notitle
1
2
3
e