Numerical study of neoclassical plasma pedestal in a tokamak geometry