A robust multi-step predictor-corrector method for solving fractional order biological models