Numerical solution of first order initial value problem using 7-stage tenth order Gauss-Kronrod-Lobatto IIIA method