plotstyle = 'o';
Frame = 10;
plotframe1

rhotrue
axis([-15 5 0 6])

print slowshock1.eps
query

axis([-15 5 5.55 5.75])
print slowshock2.eps