set data s l set autoscale xy set nolog set term x11 0 set log y # set xrange [0:100] set xlabel "iteration #" set ylabel "relative error in residual" plot "fort.15" u 1:2 title "actual err",\ "fort.15" u 1:3 title "recursive err" set terminal postscript eps enhanced color solid 22 set output "err-vs-it.eps" replot set term x11 1 set nolog set xlabel "x" #set ylabel "Phi(x,j) at j=1, 2, nx/2, nx-1, nx (slight offsets)" #plot "fort.31" u 1:($2*0.98) title "j=1" ,\ # "fort.32" u 1:($2*0.99) title "j=2" ,\ # "fort.33" u 1:($2*1.0 ) title "j=nx/2" ,\ # "fort.34" u 1:($2*1.01) title "j=nx-1" ,\ # "fort.35" u 1:($2*1.02) title "j=nx-2" kmid=int(sqrt(101)) set ylabel "Phi(x,j) & Phi_true(x,j) at j= nx/2" plot "fort.43" u 1:($2*1.0) title "Phi_True",\ "fort.33" u 1:($2*1.0) title "Phi" # "fort.43" u 1:((cos(pi*$1)/pi + cos(kmid*pi*$1)/(kmid*pi))) \ # title "Phi_cos" ,\ set term x11 2 set nolog set xlabel "iteration #" set ylabel "sin(theta)^2" plot "fort.15" u 1:5 title "sin(theta)^2" # set term x11 1 # set nolog # set xlabel "i" # set ylabel "x(i)-xtrue(i)" # plot "fort.12" u 1:($2-$3) # set term x11 2 # set nolog # set xlabel "i" # set ylabel "solution Psi" # plot "fort.12" u 1:3 title "u_true", \ # "fort.12" u 1:2 title "u"