Update solvers.f90

Added in more accurate intermediate step (i.e., 2-step solver) for the STS method