A double-well potential for levitated particles can be used to generate non-classical states of motion, which are of fundamental interest in quantum optics to generate for mechanical systems with ever increasing mass. The non-classicality of a state can be quantified using the Wigner function, which is possible to obtain through propagating an initial state using the associated time evolution equation. We have found that by using magnetic levitation, it is possible to obtain more favourable parameters than with other levitation systems and would use the compute resources to simulate the dynamics of the Wigner function under these parameter ranges, including realistic decoherence channels.