Oscillations between active and sterile neutrinos remain as an open possibility to explain some experimental observations. In a four-neutrino mixing scheme, we use the Magnus expansion of the evolution operator to study the evolution of neutrino flavor amplitudes within the Earth. We apply this formalism to calculate the transition probabilities from active to sterile neutrinos taking into account the matter effect for a varying terrestrial density.