MATLAB: Could anyone explain to me the method used by the ode23s function

functionmatlab functionnumerical integrationodes-function

Hi,
Could anyone explain to me the method used by the ode23s function? I know it is applied to rigid differential equations as well as other solvers. But I did not find the material that explained in more detail about the numerical method in question. Is it only applied to continuous systems? What's your difference with ode15s?