Code a numerical differential equation solver and solve the SIR model for N=100 until t=10 with a step size of t=0.1. Please make clear the plot of the solution, its code and, if possible, proof of successful implementation.

The initial terms given were

S₀=0.95N,

I₀=0.05N,

R₀=0;

= 0.6

=0.5

Thank you in advance for your time and help.