Text this: Numerical solution of first order initial value problem using 4-stage sixth order Gauss-Kronrod-Radau 1 method