for Frame = [0 1 2 3 4 5 10 15 20 25 30];
   plotframe1
   str = ['print decay2N' num2str(Frame) '.eps'];
   eval(str)
   query
   end