Abstarct: Nowadays attempts to detect and achieve buried structures and underground resources have developed widely and geophysics is a means of identifying these structures. The main goal of geophysical methods is incorporating additional information to the pros of inversion to define the characteristics of geological structures as precisely as possible. To reach this goal, 3D modeling is one of the most successful methods of interpreting geophysical data. The outcome of modeling is to have a better understanding of target such as the shape and depth of source. Finding these parameters influences the subsequent decisions directly which can influence the management of major costs. In this study we have tried to use two-dimensional inversion (Geosoft oasis montaj software) and 3D constrained inversion of magnetic data with Li-Oldenburg method (MAG3D software), correctly interpret the results. This algorithm includes a multi-coefficient objective function which is flexible enough to produce different kinds of models which are in good agreement with geological structures. For optimized use of the algorithm at first it was tested on artificial models which were contaminated with noise, then according to the accommodation of the results with the right model it was used for the inversion of real data. In this research, magnetic data of kharanq in Yazd area, were inverted and the final results which have been shown in different views. The results indicate the presence of three zones of mineralization in the North East, South East and West range. West masses larger area than the other two areas and further depth expand its estimated. Depth to the top and down the main mass anomaly approximately 5 and 90 meters in West and 1 and 30 meters in North-East respectively and southeast anomaly has been expanded from the surface to a depth of 10 meters.