For simplicity, we will set the resolutions to δ δ x = h m ω ω δ δ p = h m ω ω so that the probabilities reduce to P ⁡ ⁡ [ x j ] = P ⁡ ⁡ [ p j ] = 1 2 [ erf ⁡ ⁡ ( ( j + 1 2 ) 2 π π ) − − erf ⁡ ⁡ ( ( j − − 1 2 ) 2 π π ) ]