Variational methods via supersymmetric techniques