Numerical method for finding the convex hull of an inclusion in conductivity from boundary measurements

We consider the 2D inverse conductivity problem for conductivities of the f orm gamma = 1 + chi(D)h defined in a bounded domain Omega subset of R-2 wit h C-infinity boundary partial derivative Omega. Here D subset of Omega and h is an element of L-infinity (D) are such that gamma has a jump along part ial derivative D. It was shown by Ikehata (Ikehata M J. Inverse and Ill-Pos ed Problems at press) that the Dirichlet-Neumann map determines the indicat or function I-omega(tau, t) that can be used to find the convex hull of D. In this paper we find numerically the indicator function for examples with constant h and recover the convex hull of D.