NOTE: S= y.T[0] I= y.T[0] R= y.T[2]
These are matrix obtained by transposeing the resultant matrix of sim function and odeint. For short cut(.T)method was used.