Code for finding isodensity contours