My Current Project
 Given a planar curve which begins and ends on an axis (in this case the z axis)
 numerically generate the initial data for a nonrotating vacuum spacetime which has as a marginally trapped surface, the surface of revolution of the above planar curve.

Some comments
 The spacetime generated by the evolution of the initial data contains an apparent horizon, which is not necessarily the specified marginally trapped surface. Because it contains an apparent horizon, it will automatically contain an event horizon.
 The ultimate code will be generic enough to produce most apparent horizons, excluding ones made from closed planar curves. Those would be an obvious and simple enough extension.
