set title "morphing damped"
set ylabel 'Centre of Mass'
set xlabel 'time'

set term 'png'
set output 'zMorphingDamped.png'

plot "cylindermoveMeshsymplectic/plotfile" u 1:($4/0.0016)   title 'symplectic' with lines lt 1, "cylindermoveMeshOuterCorrectorsYes/plotfile" u 1:($4/0.0016)   title 'Newmark 2 outer' with lines lt 2, "cylindermoveMeshOuterCorrectorsYes4Outer/plotfile" u 1:($4/0.0016)   title 'Newmark 4 outer' with lines lt 3, "cylindermoveMeshOuterCorrectorsNo/plotfile" u 1:($4/0.0016)   title 'Newmark no outer correction' with lines lt 7, "cylindermoveMeshOuterCorrectorsYes4Outer/plotfile" u 1:(exp(-0.5*0.0039/0.03575*$1)*cos(sqrt(69.48/0.03575)*$1))   title 'analytical solution' with lines lt -1,


    


