Icing simulations are of importance in the development of aircraft and jet engines. Currently, icing simulations are performed in the Eulerian framework. On the other hand, we employed the MPS method which is a mesh-free method in the Lagrangian framework. In addition, the MPS method can predict liquid interface if its deformation is large. Our final goal is to perform the three-dimensional icing simulation on NACA airfoil. In the present paper, we perform icing of one droplet as for a preliminary. The obtained result indicates that the three-dimensional icing simulation using MPS method reproduces the complex ice shape even in the two-dimensional simulation. If angle of attack is large, splash mass computed by the MPS method shows good agreement with that by the previous studies.