a=0.0001
b=1.
E0=1.433268e-6
E1=1.492525e-7
E2=1.50796e-8
R0(x)=abs(1-16.*x/E0)
R1(x)=abs(1-16.*x/E1)
R2(x)=abs(1-16.*x/E2)
RR0(x)=10*abs(16.*x/E0)
RR1(x)=10*abs(16.*x/E1)
RR2(x)=10*abs(16.*x/E2)
set border 1
set noxtics; set noytics; set xtics nomirror; set ytics nomirror
set key below
set xlabel "alpha"
set ylabel ""
set title "nergie en fonction de alpha\n pour t/L = 1/100 et N = 4"
set terminal x11 0
plot [a:b] \
'4a/cy-100.alpha' using ($1):(RR0($6)) t "10 X |E_c - E_nc|/E" with l lw 2, \
'4a/cy-100.alpha' using ($1):(R0($4)) t "|1 - E_c/E|" with l lw 2
set title "nergie en fonction de alpha\n pour t/L = 1/100 et N = 8"
set terminal x11 1
plot [a:b] \
'8a/cy-100.alpha' using ($1):(RR0($6)) t "10 X |E_c - E_nc|/E" with l lw 2, \
'8a/cy-100.alpha' using ($1):(R0($4)) t "|1 - E_c/E|" with l lw 2
set title "nergie en fonction de alpha\n pour t/L = 1/100 et N = 16"
set terminal x11 2
plot [a:b] \
'16a/cy-100.alpha' using ($1):(RR0($6)) t "10 X |E_c - E_nc|/E" with l lw 2, \
'16a/cy-100.alpha' using ($1):(R0($4)) t "|1 - E_c/E|" with l lw 2

#set title "nergie en fonction de alpha\n pour t/L = 1/100 et N = 32"
#set terminal x11 3
#plot [a:b] \
#'24a/cy-100.alpha' using ($1):(RR0($6)) t "10 X |E_c - E_nc|/E" with l lw 2, \
#'24a/cy-100.alpha' using ($1):(R0($4)) t "|1 - E_c/E|" with l lw 2

set title "nergie en fonction de alpha\n pour t/L = 1/1000 et N = 4"
set terminal x11 4
plot [a:b] \
'4a/cy-1000.alpha' using ($1):(RR1($6)) t "10 X |E_c - E_nc|/E" with l lw 2, \
'4a/cy-1000.alpha' using ($1):(R1($4)) t "|1 - E_c/E|" with l lw 2
set title "nergie en fonction de alpha\n pour t/L = 1/1000 et N = 8"
set terminal x11 5
plot [a:b] \
'8a/cy-1000.alpha' using ($1):(RR1($6)) t "10 X |E_c - E_nc|/E" with l lw 2, \
'8a/cy-1000.alpha' using ($1):(R1($4)) t "|1 - E_c/E|" with l lw 2
set title "nergie en fonction de alpha\n pour t/L = 1/1000 et N = 16"
set terminal x11 6
plot [a:b] \
'16a/cy-1000.alpha' using ($1):(RR1($6)) t "10 X |E_c - E_nc|/E" with l lw 2, \
'16a/cy-1000.alpha' using ($1):(R1($4)) t "|1 - E_c/E|" with l lw 2

#set title "nergie en fonction de alpha\n pour t/L = 1/1000 et N = 32"
#set terminal x11 7
#plot [a:b] \
#'24a/cy-1000.alpha' using ($1):(RR1($6)) t "10 X |E_c - E_nc|/E" with l lw 2, \
#'24a/cy-1000.alpha' using ($1):(R1($4)) t "|1 - E_c/E|" with l lw 2

set title "nergie en fonction de alpha\n pour t/L = 1/10000 et N = 4"
set terminal x11 8
plot [0.000001:1] \
'4a/cy-10000.alpha' using ($1):(RR2($6)) t "10 X |E_c - E_nc|/E" with l lw 2, \
'4a/cy-10000.alpha' using ($1):(R2($4)) t "|1 - E_c/E|" with l lw 2
set title "nergie en fonction de alpha\n pour t/L = 1/10000 et N = 8"
set terminal x11 9
plot [a:b] \
'8a/cy-10000.alpha' using ($1):(RR2($6)) t "10 X |E_c - E_nc|/E" with l lw 2, \
'8a/cy-10000.alpha' using ($1):(R2($4)) t "|1 - E_c/E|" with l lw 2
set title "nergie en fonction de alpha\n pour t/L = 1/10000 et N = 16"
set terminal x11 10
plot [1e-10:1] \
'16a/cy-10000.alpha' using ($1):(RR2($6)) t "10 X |E_c - E_nc|/E" with l lw 2, \
'16a/cy-10000.alpha' using ($1):(R2($4)) t "|1 - E_c/E|" with l lw 2

#set title "nergie en fonction de alpha\n pour t/L = 1/10000 et N = 32"
#set terminal x11 11
#plot [a:1] \
#'24a/cy-10000.alpha' using ($1):(RR2($6)) t "10 X |E_c - E_nc|/E" with l lw 2, \
#'24a/cy-10000.alpha' using ($1):(R2($4)) t "|1 - E_c/E|" with l lw 2
