ci2d function

Create 2-dimensional empirical confidence regions