Abstract
The quantum behavior of electrons in bilayer graphene with applied magnetic fields is addressed. By using second-order supersymmetric quantum mechanics the problem is transformed into two intertwined one dimensional stationary Schr\"odinger equations whose potentials are required to be shape invariant. Analytical solutions for the energy bound states are obtained for several magnetic fields. The associated spectrum is analyzed, and the probability and current densities are determined.