The electric potential distribution around a dust particle immersed in a magnetised supersonic plasma flow is studied by numerical simulations. It is shown that with increasing magnetisation of plasma, the peak in the wake potential gets smaller and moves upstream. For strong magnetisation, the trailing peak in the potential distribution vanishes and the potential becomes more isotropic. The results agree qualitatively with the linear response approach. The numerical simulations are carried out with a particle-in-cell code.