Block predictor-corrector schemes for the parallel solution of ODEs

A fourth-order block method based on the composite Simpson rule is developed for the parallel solution of ordinary differential equations. Like the block scheme based on the composite Trapezoidal Rule, its principal error term is linear in the block size while the increased order and stability allow a modest increase in parallelism without further computational complexity. Numerical results confirm the enhanced properties of the higher-order method.