SPECIAL ALGORITHM FOR THE NUMERICAL SOLUTION OF SYSTEM OF INITIAL VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL EQUATIONS USING BLOCK HYBRID EXTENDED TRAPEZOIDAL RULE OF SECOND KIND
We develop self-starting family of three and five step continuous extended trapezoidal rule of second kind a block hybrid type (BHETR2s) methods through interpolation and collocation procedure. The BHETR2s methods are then used to produce multiple numerical integrators which are each of the same order and assembled into a single block matrix equation. These equations are simultaneously applied to provide the approximate solution for the ordinary differential equations. The stability properties of the methods were investigated and found to be consistent, zero-stable and hence convergent. The block integrators were tested on three numerical initial value problems of ODEs to show accuracy and efficiency.