/* -*-ePiX-*- */
#include "epix.h"
using namespace ePiX;
const double C1(cosh(1));

const double mass(0.1);

P f(double x, double y)
{
  return P(x + 0.5*mass*pow(y/x, 2), y*(1 - mass/x));
}

P f(const P& arg) { return f(arg.x1(), arg.x2()); }

/*
// inverse map
P g(double u, double v)
{
  double y(Atan2(v, u+C1));
  return P(acosh(sqrt((u+C1)*(u+C1) + v*v)*pow(Cos(y), 3)), y/0.75);
}

P g(const P& arg) { return g(arg.x1(), arg.x2()); }
*/

domain R(P(0.5, -2), P(2,2), mesh(4, 5), mesh(4, 60));

int main()
{
  picture(P(0,-3), P(6,3), "3.5 x 3.5in");

  begin();

  label(P(0.5*(xmin()+xmax()), ymin()), P(0,-8), "\\textsc{Fig.~14.}", b);

  P ptP(1.125, -0.4), ptQ(0.8, 0.4);

  bold();

  plot(f, R);

  dashed();
  dash_size(6);
  line(f(ptP), f(ptQ));

  font_size("footnotesize");
  label(f(ptP), P(2,0), "$P$", r);
  label(f(ptQ), P(2,0), "$Q$", r);

  tikz_format();
  end();
}