The following system is a classic example of stiff ODEs that can occur in the solution of chemical reaction kinetics:

dc₁/dt = – 0.013c₁ – 1000c₁c₃

dc₂/dt = –2500c_2c3

dc₃/dt = –0.013c₁ – 1000c₁c₃ – 2500c₂c₃

Solve these equations from t = 0 to 50 with initial conditions c₁(0) = c₂(0) = 1 and c₃(0) = 0. If you have access to MATLAB software, use both standard (for example, ode45) and stiff (for example, ode23s) functions to obtain your solutions.