| 
      Tol1 and Tol2 represent the ||p|| norm value for each algorithm with P1 and P2 protocols. ΔP1 Tol1/Tol2 represents the absolute difference between Tol1 and Tol2 for P1,
    whereas Δ P2 Tol1/Tol2 represents the same absolute difference for P2. Δ Tol1 P1/P2 represents the absolute difference between P1 and P2 for Tol1 parameters, whereas
    ΔTol2 P1/P2 represents the same absolute difference for Tol2.
        |     | RK4   | RKF   | TR- BDF2 | Rosen-  brock | IMEX   | JVODE  ADAMS | JVODE  BDF | JLSODE   |  
        | Protocol 1 (single pulse) |  
        | Tol1 | 1.0 | 5.004e-2 | 5.154e-2 | 5.006e-2 | 5.007e-2 | 5.008e-2 | 5.008e-2 | 5.003e-2 |  
        | Tol2 | 1.0 | 5.004e-2 | 0.683 | 5.006e-2 | 5.007e-2 | 5.008e-2 | 5.008e-2 | 5.003e-2   |  
        | ΔP1 Tol1/Tol2 | - | 2.5e-8 | 0.6314 | 9.84e-7 | 9.4e-8 | 1.056e-6 | 1e-6 | 1.2e-8 |  
        | Protocol 2 (repeated pulses) |  
        | Tol1 | 1.0 | 5.001e-2 | 5.647e-2 | 5.018e-2 | 5.019e-2 | 5.037e-2 | 5.037e-2 | 5e-2 |  
        | Tol2 | 1.0 | 5.001e-2 | 3.215 | 5.018e-2 | 5.019e-2 | 5.037e-2 | 5.037e-2 | 5e-2 |  
        | ΔP2 Tol1/Tol2 | - | 2.7e-10 | 3.1587 | 2.8e-8 | 1.08e-7 | 6.4e-8 | 2e-9 | 1e-9 |  
        | ΔTol1 P1/P2 | - | 5.976e-6 | 4.933e-3 | 1.198e-5 | 1.192e-5 | 2.988e-5 | 2.899e-5 | 5.99e-6 |  
        | ΔTol2 P1/P2 | - | 6e-6 | 2.5322 | 1.102e-5 | 1.213e-5 | 3.006e-5 | 3e-5 | 6.01e-6 |  |