The determination of the differential branching fraction and the first angular analysis of the decay B0s → ϕμ + μ − are presented using data, corresponding to an integrated luminosity of 1.0 fb−1, collected by the LHCb experiment at s√=7 TeV. The differential branching fraction is determined in bins of q 2, the invariant dimuon mass squared. Integration over the full q 2 range yields a total branching fraction of B(B0s→ϕμ+μ−)=(7.07+0.64−0.59±0.71±0.71) × 10−7, where the first uncertainty is statistical, the second systematic, and the third originates from the branching fraction of the normalisation channel. An angular analysis is performed to determine the angular observables F L, S 3, A 6, and A 9. The observables are consistent with Standard Model expectations.