A method for modeling a three dimensional shape of object using a level set solution on a partial differential equation derived from a Helmholtz reciprocity condition is provided. The method includes the steps of: a) inputting an image pair satisfying a Helmholtz reciprocity condition; b) performing an optical correction and simultaneously performing a geometric correction; c) performing a camera selection to select cameras capable of seeing a point (X, Y, Z), and defining and calculating a cost function by the Helmholtz reciprocity condition; d) calculating a speed function of a PDE that minimizes the cost function obtained in the step c); and e) generating a three dimension mesh model from a set of points configuring the object surface provided from the step d), and deciding a final three dimension mesh model by comparing cost function values.