\
set term svg; \
set out "Plots.svg"; \
set multiplot layout 2,1; \
set title "Hydrogen s-wave shooting to F_E(r_{max})=0"; \
set tics out; \
set xtics 2; \
set mxtics 2; \
set xlabel "r [a_{B}]"; \
set ylabel "F_E(r)"; \
set size 0.7,0.5 ;\
plot [][0:]\
 "data" index 1 with lines title "r_{max}=2" \
,"data" index 3 with lines title "r_{max}=3" \
,"data" index 5 with lines title "r_{max}=4" \
,"data" index 9 with lines title "r_{max}=10" \
,"data" index 9 using 1:3 with lines linetype 2 title "exact"; \
set mytics 2; \
set title "Calculated energy {/:Italic E} as function of {/:Italic r_{max}}"; \
set xlabel "{/:Italic r_{max} [a_{B}]}"; \
set ylabel "E [Hatree]"; \
set size 0.7,0.5 ;\
plot [2-1:10+1][-0.6:0]\
 "data" index 0 with points pointtype 7 linecolor rgb "blue" title "calculated"\
,"data" index 2 with points pointtype 7 linecolor rgb "blue" notitle \
,"data" index 4 with points pointtype 7 linecolor rgb "blue" notitle \
,"data" index 6 with points pointtype 7 linecolor rgb "blue" notitle \
,"data" index 8 with points pointtype 7 linecolor rgb "blue" notitle \
,-0.5 with lines linetype 2 title "exact"\