[en] Complete text of publication follows. In order to explore the borders of island of inversion we investigated the nuclei ^30,31Na and ^33,34Mg by inelastic proton scattering in inverse kinematics. A ⁴⁰Ar primary beam of 94 MeV/nucleon energy with 60 pnA intensity was transported to a ¹⁸¹Ta production target of 0.5 mm thickness. The RIPS fragment separator analyzed the momentum and mass of the reaction products. The ^30,31Na and ^33,34Mg beams were produced in two individual runs with different settings of Bρ values. In the first run of ^30,31Na, the secondary beam included neutron-rich O, F, Ne and Na nuclei with A=Z ∼ 3 while mainly Mg and Al isotopes were mixed in the second run of ^33,34Mg. The total intensity was about 100 particle/s (pps) for both runs, while the ^30,31Na/^33,34Mg intensities reached 8, 6, 3 and 2 pps, respectively, on average. The identification of incident beam species was performed on an event-by-event basis by means of energy loss, time-of-flight (TOF) and magnetic rigidity (Bρ). The secondary beam hit a liquid hydrogen target. The mean energy of the isotopes in the target was around 50 MeV/nucleon. Two PPACs at F3 upstream of the target monitored the position of the incident particles. The reaction products and scattered particles were detected and identified by a PPAC and a silicon telescope. The Z identification was performed by TOF-energy loss method where the TOF was taken between the PPACs upstream and downstream of the secondary target. The isotope separation was done by use of the ΔE-E method. The de-exciting γ rays emitted by the inelastically scattered nuclei were detected by the DALI2 setup consisting of 146 NaI(Tl) scintillators surrounding the target (see, e.g., Figure 1). To determine the cross sections of the production of the γ rays in proton inelastic scattering the peak positions determined were fed into the detector simulation software GEANT4 and the resultant response curves plus smooth polynomial backgrounds were used to analyze the experimental spectra. From a distorted wave analysis of the cross sections, we derived 'matter' deformation parameters (βM), which are consistent with the charge deformations determined from Coulomb excitation experiments. Based on the 'matter' and charge deformation parameters, the neutron deformation parameters could also be extracted. The results show that all these nuclei are largely deformed; the deformation of the proton and neutron distributions are similar and cannot be distinguished at the present experimental uncertainties. (author)