A Study of Numerical Method for Solution of Multistep Parameter

This paper considers the numerical solution of delay differential equations. The predictor–corrector scheme based on generalized multistep methods are implemented in variable order variable step size techniques. The formulae are represented in divided difference form where the integration coefficients are computed by a simple recurrence relation. This representation produces simpler calculation as compared with the modified divided difference form, but no sacrifice is made in efficiency and accuracy of the method. Numerical results prove that the method is reliable, efficient and accurate. The P- and Q-stability regions for a fixed step size of the predictor–corrector scheme are illustrated for various orders.