calculation of force requires Quantum optics toolbox.
The simulation is based on QOtoolbox example xprobmulti.m.

parameter:
Rabi freq = 0.5Gamma
Delta=-3Gamma
trap=851nm