/* -*-ePiX-*- */
/* anim.xp, April 02, 2003 */
#include "epix.h"
using namespace ePiX;

P F(double s, double t)
{ 
  return (P(0.1*s, 0.025*t) + (1/(0.01 + s*s + t*t))*P(-t, s)); 
}

main() 
{
  unitlength("1pt");
  picture(100, 100);
  bounding_box(P(-3, -3), P(6, 6));

  for (int i=0; i<6; ++i)
    {
      begin();
      crop();
      dart_field(F, P(0.8*x_min, 0.8*y_min), P(0.8*x_max, 0.8*y_max), 8, 8);

      bold();
      rect(P(x_min, y_min), P(x_max, y_max));

      for (int j=0; j<7; ++j)
	{
	  rgb(0.25 - 0.05*j, 1 - 0.1*j, 0.2*j);
	  ode_plot(F, P(0.85-0.025*j, 0), 6*i + 1, 90);
	}

      std::cout << "\n\\end{picture}";
      if (i == 2)
	std::cout << "%\n\\vspace*{2ex}";
    }
}