The aim of this study is to examine dust dynamics on a large scale and investigate the coupling of dust with gas fluid in the star formation process. We propose a method for calculating the dust trajectory in a gravitationally collapsing cloud, where the dust grains are treated as Lagrangian particles and are assumed to be neutral. We perform the dust trajectory calculations in combination with non-ideal magnetohydrodynamics simulation. Our simulation shows that dust particles with a size of are coupled with gas in a star-forming cloud core. We investigate the time evolution of the dust-to-gas mass ratio and the Stokes number, which is defined as the stopping time normalized by the freefall time-scale, and show that large dust grains () have a large Stokes number (close to unity) and tend to concentrate in the central region (i.e. protostar and rotationally supported disc) faster than do small grains (). Thus, large grains significantly increase the dust-to-gas mass ratio around and inside the disc. We also confirm that the dust trajectory calculations, which trace the physical quantities of each dust particle, reproduce previously reported results obtained using the Eulerian approach.