An angular analysis of the rare decay $Bs0\to \varphi \mu +\mu -$ is presented, using proton-proton collision data collected by the LHCb experiment at centre-of-mass energies of $7$, $8$ and $13 TeV$, corresponding to an integrated luminosity of $8.4 fb-1$. The observables describing the angular distributions of the decay $Bs0\to \varphi \mu +\mu -$ are determined in regions of $q2$, the square of the dimuon invariant mass. The results are consistent with Standard Model predictions.