For numerical solution of laser induced vaporization of metal, we have to vesearch of calcolation method of fluid dynamics including phase transition. Because the vaporization zone is movement as time, so that, this problem is different from general fluid dynamics with a moving zone.This paper considers the energy deposition as zone deposition. According to the conditions of first class phase transition, outside the vaporization zone, we adppt the von Neumann method, and inert the vaporization zone, the another method is used. We have established the program named LHAP -1DVG, in which this method is applied.The results of a lot of calcolations indicate that this method is successful.