Inductive coupling of magnetization dynamics in ferromagnets with the diamagnetic response of superconductors leads to a number of notable phenomena that may become useful for applications in magnonics. In magnetostatically coupled ferromagnet/superconductor bilayers spin waves propagate with substantially enhanced phase velocities. The enhancement occurs due to diamagnetic screening of stray fields at the interface by Meissner currents and, accordingly, affecting the actual effective magnetic field that acts on the magnetic moments1,2. Capability of superconductors to screen the alternating magnetostatic stray fields or the external field can also be employed for development of artificial metamaterials with periodically modulated spin-wave dispersion, known as magnonics crystals3. Hybrid magnonics crystals exhibit wide forbidden bands for spin wave propagation at the GHz range and a rich dependence of the band structure on geometry of the system. All these phenomena can be considered numerically using modified micromagnetic simulations1-3.