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

double f(double u)
{
  double temp=sqrt(1-exp(-2*u));

  return -temp + atanh(temp);
}

double g(double u)
{
  return exp(-u);
}

double fm(double u)
{
  return -f(u);
}

int main()
{
  bounding_box(P(-2,-2),P(2,2));
  unitlength("1in");
  picture(4,4);

  begin();
  grid();

  domain R(P(0,0), P(2,1), mesh(24, 36), mesh(24, 36));
  frame F(E_3, E_1, E_2);

  camera.at(P(0,-10,5));

  fill();
  surface_rev(fm, g, R, F);
  surface_rev(f, g, R, F);

  end();
}