“…Let us consider the following delay differential equation: For sake of comparison of our delayed model to the non delayed model introduced in ( [8]), we are using the same coefficients in our numerical study. Also, as in ( [8]), we take the initial temperature of the slag to be T in = T (0) = 1773K. Thus, in our following numerical investigation, we will consider If the delay τ is larger enough than the step of integration, the above fourth Runge-Kutta method can be extended to a fourth Runge Kutta algorithm with delay by expressive constants K1 , K2 , K3 and K4 with the delayed form of the f (t, T (t − τ ), T (t)).…”