Development of Multistep and Degenerate Variational Integrators for Applications in Plasma Physics