We investigate the angular decay distribution of the four-fold Bc→Ds⁎(→Dsγ)μ+μ−, and Bc→Ds⁎(→Dsπ)μ+μ− decays that proceed through b→sμ+μ− quark level transition. We use the model independent effective Hamiltonian with vector and axial vector new physics operators to formulate the angular observables and study the implications of different latest new physics scenarios, taken from the global fits to all the b→s data, on these observables. We also give Standard Model and new physics predictions of several observables such as differential branching ratios, forward backward asymmetry, longitudinal polarization fraction of Ds⁎, and the unpolarized and polarized lepton flavor universality violating ratios. Future measurements of the predicted angular observables, both at current and future high energy colliders, will add to the useful complementary data required to clarify the structure of new physics in b→sℓℓ neutral current decays.