Vectorize a system of ODE's and arrays

Question

I have a simple SIR model. I have the parameter values gamma, beta, N and mu which are currently defined as single numeric values as shown below.

    gamma =