Three-and four-derivative Hermite-Birkhoff-Obrechkoff solvers for stiff ODE