
load traj30;
load traj30_9;

global scale dist f_text
get_scale;
scale = sqrt(1.37*dist);
f_text = sprintf('f=%g', (scale.^2)/dist);
set_contour(.3);

%%%%% family 2: two plots (T=7.4,7.5)

figure;
ret_7_4;
set(gca, 'XLim',[-.7 1.3], 'YLim',[-.5 .6+eps], 'DataAspectRatio',[1 1 1]);

global energy_contour f_text
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
hold on
plot(x2,y2,'g');
title(sprintf('T=7.4; family 2 (\\tau=%g)', t2))
hold off
print -deps2 over_7_4.eps

figure;
ret_7_5;
set(gca, 'XLim',[-.7 1.3], 'YLim',[-.5 .6+eps], 'DataAspectRatio',[1 1 1]);

global energy_contour f_text
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
hold on
plot(x2,y2,'g');
title(sprintf('T=7.5; family 2 (\\tau=%g)', t2))
hold off
print -deps2 over_7_5.eps

%%%%% ...plus one of family 9 (T=7.4)

figure;
ret_7_4;
set(gca, 'XLim',[-.7 1.3], 'YLim',[-.5 .6+eps], 'DataAspectRatio',[1 1 1]);

global energy_contour f_text
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
hold on
plot(x9,y9,'g');
plot(x2(54:212),y2(54:212),'g-'); % sigh... =(
title(sprintf('T=7.4; family 9 (\\tau=%g)', t9))
hold off
print -deps2 over_7_4_with9.eps

%%%%% family 4: one plot (T=8.0)

figure;
ret_8;
set(gca, 'XLim',[-.7 1.3], 'YLim',[-.5 .6+eps], 'DataAspectRatio',[1 1 1]);

global energy_contour f_text
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
hold on
plot(x4,y4,'g');
title(sprintf('T=8.0; family 4 (\\tau=%g)', t4))
hold off
print -deps2 over_8.eps

%%%%% family 4 'shadow': two plots (T=6.9,7.0)

figure;
ret_6_9;
set(gca, 'XLim',[-.7 1.35], 'YLim',[-.5 .6+eps], 'DataAspectRatio',[1 1 1]);

global energy_contour f_text
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
hold on
plot(xD,yD,'g');
title(sprintf('T=6.9 at E=0.3; family 4 at E=%g (\\tau=%g)', energyD, tD))
hold off
print -deps2 over_6_9.eps

figure;
ret_7;
set(gca, 'XLim',[-.7 1.35], 'YLim',[-.5 .6+eps], 'DataAspectRatio',[1 1 1]);

global energy_contour f_text
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
hold on
plot(xD,yD,'g');
title(sprintf('T=7.0 at E=0.3; family 4 at E=%g (\\tau=%g)', energyD, tD))
hold off
print -deps2 over_7.eps
