A Monte Carlo study of the diffusion in Spheromak is described. Using appropriate scattering operator to simulate accumulated small angle Fokker-Plank collisions, the diffusion coefficient in Spheromak is estimated with varying collision frequency. Also the orbit of a particle for a given model magnetic field is computed and verified the possibility of the existence of neoclassical diffusion in Spheromak.