% Script EllipsePlot
% Plots 100 different ellipses

  t = linspace(0,2*pi,200);
  c = cos(t);
  s = sin(t);
  figure 
  hold on
  for b = 1:-.01:.01
      plot(c,b*s,'r')
      axis equal off manual
      title(sprintf('The ellipse x^2 + (y/b)^2 = 1 with b = %5.2f',b),'Fontsize',16)
      shg
      pause(.1)
  end
  hold off