Hi, I just need help writing a general code in R that could solve the following system of ODE's as a function of time. I will put in the values and such later, I mainly need help with the general set up and structure of how the code should be written in R. Again, I just need a basic set up, so I will put in the values later. I just need to see how it should be done in R. Thanks!

Consider a population aggregated by age into M groups ( i = 1,2,,M). The population, within age group i, is comprised of the following classes, susceptible S_i , asymptomatic infectious Ia_i, symptomatic infectious Is_i, and recovered R_i individuals. The governing equations are:

such that

Si = -lit)Si, in = ali(t) S-710 yl, i =li(t)S-71, R; = (1+1). I (t) = ii Nj + C ij Nj 1,3 = 1,... j=1 Si = -lit)Si, in = ali(t) S-710 yl, i =li(t)S-71, R; = (1+1). I (t) = ii Nj + C ij Nj 1,3 = 1,... j=1