/* -*-ePiX-*- */
#include "epix.h"
using namespace ePiX;

const double r1(0.5);
const double r2(1.5);
const double ht(0.5);

double f(double t)
{
  return ht + 0.25/((r1 - t)*(r2 - t));
}

int main()
{
  picture(P(-1, -2), P(3, 2), "1.5 x 1.5in");

  begin();
  P pO(0,0), xi(xmax(), 0), eta(0, ymax());

  bold();

  line(P(xmin(), 0), xi);
  line(-eta, eta);
  set_crop();
  plot(f, xmin(), xmax(), 120);

  dashed();
  dash_size(6);
  line(P(r1, ymin()), P(r1, ymax()));
  line(P(r2, ymin()), P(r2, ymax()));
  line(P(xmin(), ht), P(xmax(), ht));

  label(pO, P(-2,-4), "$O$", bl);
  label(xi, P(0,-4), "$X$", b);
  label(eta, P(-2, 0), "$Y$", l);

  label(P(r1,0), P(2, -4), "$x_{1}$", br);
  label(P(r2,0), P(2, -4), "$x_{2}$", br);

  tikz_format();
  end();
}