a=-1.6
b=-0.6
K=0.5-2*b
K2=0.5-4*b
Z=0.02
E0=1.433268e-6
E1=1.491515e-7
E2=1.50796e-8
g0(x)=abs(1.-16.*x/E0)
g1(x)=abs(1.-16.*x/E1)
g2(x)=abs(1.-16.*x/E2)
gg0(x)=abs(E0-16.*x)
gg1(x)=abs(E1-16.*x)
gg2(x)=abs(E2-16.*x)
set border 3
set noxtics; set noytics; set xtics nomirror; set ytics nomirror
set key below
set xlabel "N"
set ylabel "|1 - E_c/E|"
set title "Convergence t/R = 1/100"
set terminal x11 0
plot [4.:32.] \
'loglog01' using ($1):(g0($8)) t "Min" with l lw 2, \
'loglog01.1' using ($1):(g0($8)) t "Max" with l lw 2, \
set border 31
set noxtics; set noytics; set xtics mirror; set ytics mirror
set ylabel "log(|1 - E_c/E|)"
set xlabel "log(1/N)"
set terminal x11 1
plot [a:b] 2*x+K t "linaire" with l lw 2, 4*x+K2 t "quadratique" with l lw 2, \
'loglog01' using (log10(1./$1)):(log10(g0($8))) t "Min" with l lw 2, \
'loglog01.1' using (log10(1./$1)):(log10(g0($8))) t "Max" with l lw 2, \
set border 3
set noxtics; set noytics; set xtics nomirror; set ytics nomirror
set key below
set xlabel "N"
set ylabel "|1 - E_c/E|"
set title "Convergence t/R = 1/1000"
set terminal x11 2
plot [4.:32.] \
'loglog001' using ($1):(g1($8)) t "Min" with l lw 2, \
'loglog001.1' using ($1):(g1($8)) t "Max" with l lw 2, \
'loglog001.2' using ($1):(g1($8)) t "Min_2" with l lw 2, \
'loglog001.3' using ($1):(g1($8)) t "Max_2" with l lw 2
set border 31
set noxtics; set noytics; set xtics mirror; set ytics mirror
set ylabel "log(|1 - E_c/E|)"
set xlabel "log(1/N)"
set terminal x11 3
plot [a:b] 2*x+K t "linaire" with l lw 2, 4*x+K2 t "quadratique" with l lw 2, \
'loglog001' using (log10(1./$1)):(log10(g1($8))) t "Min" with l lw 2, \
'loglog001.1' using (log10(1./$1)):(log10(g1($8))) t "Max" with l lw 2, \
'loglog001.2' using (log10(1./$1)):(log10(g1($8))) t "Min_2" with l lw 2, \
'loglog001.3' using (log10(1./$1)):(log10(g1($8))) t "Max_2" with l lw 2
set border 3
set noxtics; set noytics; set xtics nomirror; set ytics nomirror
set key below
set xlabel "N"
set ylabel "|1 - E_c/E|"
set title "Convergence t/R = 1/10000"
set terminal x11 4
plot [4.:32.] \
'loglog0001' using ($1):(g0($8)) t "Min" with l lw 2, \
'loglog0001.1' using ($1):(g0($8)) t "Max" with l lw 2, \
'loglog0001.2' using ($1):(g0($8)) t "Min_2" with l lw 2
set border 31
set noxtics; set noytics; set xtics mirror; set ytics mirror
set ylabel "log(|1 - E_c/E|)"
set xlabel "log(1/N)"
set title "Convergence"
set terminal x11 5
plot [a:b] 2*x+K t "linaire" with l lw 2, 4*x+K2 t "quadratique" with l lw 2, \
'loglog0001' using (log10(1./$1)):(log10(g0($8))) t "Min" with l lw 2, \
'loglog0001.1' using (log10(1./$1)):(log10(g0($8))) t "Max" with l lw 2, \
'loglog0001.2' using (log10(1./$1)):(log10(g0($8))) t "Min_2" with l lw 2
