function [] = return_probability_plot(time)

clf
eval([ 'ret_' time ]);
set(gca, 'XLim',[-.7 1.3], 'YLim',[-.5 .6+eps], 'DataAspectRatio',[1 1 1]);

global energy_contour f_text f_radius
if ~isempty(energy_contour),
  line(energy_contour(1,:),energy_contour(2,:), 'Color', [1 0 0]);
end
if ~isempty(f_text),
  text(1.3, .62, f_text, 'FontSize', 7.5, 'FontUnits', 'Points', ...
      'VerticalAlignment', 'baseline', 'HorizontalAlignment', 'right' );
end
if ~isempty(f_radius),
  hold on
  theta = [0:10:360] * pi/180;
  plot(1.3-f_radius*(cos(theta)+1), .6-f_radius*(sin(theta)+1));
  text(1.3-f_radius, .6-f_radius, 'P=1', ...
      'FontSize', 5, 'FontUnits', 'Points', ...
      'VerticalAlignment', 'middle', 'HorizontalAlignment', 'center');
end
hold off
drawnow

%orient landscape;
print('-deps2', [ 'returns_' time ]);
