The meshing of the domain has long been the major bottleneck in performing the finite element analysis. Research efforts which are so-called meshfree methods have recently been directed towards eliminating or at least easing the requirement for meshing of the domain. In this paper, a new meshfree method for solving nonlinear boundary value problem, based on the bubble packing technique and Delaunay triangle is proposed. The method can be efficiently implemented to the problems with singularity by using formly distributed nodes.