Please write a code in python that solves this age class - SEIR model

A meningitides outbreak has been detected in a population of N = 10,000 people. We consider 3 age groups in our multi-class susceptible-exposed-infected-recovered (SEIR) model. Class Age Population % Beta 1 [0,19] 40.2% 0.4 2 [20,59] 50.5% 0.3 3 [60,100] 9.3% 0.1 The governing equations of the i=1,2,3 age groups are given as: dS(i)/dt = -1(i)S(i) dE(i)/dt = 1(i)S(i) - YE(i) di(i)/dt = yE(i) - oI(i) dR(i)/dt = oI(i) Such that 1(i) = (B(i)*I(i))/N(i) and such that S(i), E(i), I(i) and R(i) are the susceptible, exposed, infectious and recovered individuals in the i-th age group. Given the following for parameter values: o= 1/4 #standard infectious time y = 1/6 #standard exposure time Ntotal = 10,000 And the initial conditions: Class S(O) E(0) IO) R(O) 1 4020 2. 0 0 2 5050 4 0 0 3 930 1 0 0 Use Python to write a code that solves the system of equations for the 3 age groups AND plots a case versus time graph to predict the spread of the disease in the given population. A meningitides outbreak has been detected in a population of N = 10,000 people. We consider 3 age groups in our multi-class susceptible-exposed-infected-recovered (SEIR) model. Class Age Population % Beta 1 [0,19] 40.2% 0.4 2 [20,59] 50.5% 0.3 3 [60,100] 9.3% 0.1 The governing equations of the i=1,2,3 age groups are given as: dS(i)/dt = -1(i)S(i) dE(i)/dt = 1(i)S(i) - YE(i) di(i)/dt = yE(i) - oI(i) dR(i)/dt = oI(i) Such that 1(i) = (B(i)*I(i))/N(i) and such that S(i), E(i), I(i) and R(i) are the susceptible, exposed, infectious and recovered individuals in the i-th age group. Given the following for parameter values: o= 1/4 #standard infectious time y = 1/6 #standard exposure time Ntotal = 10,000 And the initial conditions: Class S(O) E(0) IO) R(O) 1 4020 2. 0 0 2 5050 4 0 0 3 930 1 0 0 Use Python to write a code that solves the system of equations for the 3 age groups AND plots a case versus time graph to predict the spread of the disease in the given population