Please provide a working code in python. Will rate thank you!

Python code for simulating an age structured SIR model Prompt: Assume we have a population with 10,000 people and with two age classes; kids and adults. Kids are f_1=25% of the population, and adults are f_2=75% of the population. Further, assume that kids make 27 contacts each day, 2/3 of which are with other kids. Assume that adults make 15 contacts each day. The contact matrix is: C = (18 12 The system of ordinary differential equations for a SIR model with age structure includes the contact matrix, and looks like this: ds; 65; C;;!;/N; dt j dl; dt = +35; CijIj/N; 1 dR; dt = +y1; such that i and j are the indices of the age classes, and N_j is the population in age group j. The quantity beta is the probability of transmission on contact. The quantity gamma is the recovery rate. Note that in this simple model we assume that beta and gamma are the same for all age classes. The reproduction number, RO=1.5, and 1/gamma=3 days, and beta=0.05. Your objective is to write a code in Python that numerically solves the system of ODE's above, for the age-structured SIR model with the two age classes and plot a case versus time graph. Python code for simulating an age structured SIR model Prompt: Assume we have a population with 10,000 people and with two age classes; kids and adults. Kids are f_1=25% of the population, and adults are f_2=75% of the population. Further, assume that kids make 27 contacts each day, 2/3 of which are with other kids. Assume that adults make 15 contacts each day. The contact matrix is: C = (18 12 The system of ordinary differential equations for a SIR model with age structure includes the contact matrix, and looks like this: ds; 65; C;;!;/N; dt j dl; dt = +35; CijIj/N; 1 dR; dt = +y1; such that i and j are the indices of the age classes, and N_j is the population in age group j. The quantity beta is the probability of transmission on contact. The quantity gamma is the recovery rate. Note that in this simple model we assume that beta and gamma are the same for all age classes. The reproduction number, RO=1.5, and 1/gamma=3 days, and beta=0.05. Your objective is to write a code in Python that numerically solves the system of ODE's above, for the age-structured SIR model with the two age classes and plot a case versus time graph