 More specifically, it uses six function evaluations to calculate fourth- and fifth-order accurate solutions.The difference between these solutions is then taken to be the error of the (fourth-order) solution.This error estimate is very convenient for adaptive stepsize integration algorithms.