please provide the code in python

I need a Python code for an age-structured SEIR Infectious Disease Model I just need a basic code where there are 3 age groups: young, middle aged, old The equations that make up the SEIR model are described by four ordinary differential equations as follows: ds = - i(t)Si(t) dt dE = $i(t)S(t) - E:(t) dt dr. = (E;(t) - 1:(t) dt dR = 1:(t). dt Such that every susceptible individual in age group i is subject to an age-specific and time-dependent force of infection: n I;(0) Pi(t) = BMij(t) j=1 Nj(t)' And such that is the probability of transmission and the contact matrix M, where M_ij(t) is the average frequency of contacts that an individual in age group i has with an individual in age group j. I will adjust the values afterwards. I just need a python code with the basic structure for how to set this up with the contact matrix and age structure and everything so that I get outputs that look kinda like this: Suscribe 001- 004- 000- Exposed tool 4001 2004 000- Who Age 0 years 8-10 yws 10 yrs 20 400- 30 1801- Recovered 04- 02- 0.0 100 Time (years) I need a Python code for an age-structured SEIR Infectious Disease Model I just need a basic code where there are 3 age groups: young, middle aged, old The equations that make up the SEIR model are described by four ordinary differential equations as follows: ds = - i(t)Si(t) dt dE = $i(t)S(t) - E:(t) dt dr. = (E;(t) - 1:(t) dt dR = 1:(t). dt Such that every susceptible individual in age group i is subject to an age-specific and time-dependent force of infection: n I;(0) Pi(t) = BMij(t) j=1 Nj(t)' And such that is the probability of transmission and the contact matrix M, where M_ij(t) is the average frequency of contacts that an individual in age group i has with an individual in age group j. I will adjust the values afterwards. I just need a python code with the basic structure for how to set this up with the contact matrix and age structure and everything so that I get outputs that look kinda like this: Suscribe 001- 004- 000- Exposed tool 4001 2004 000- Who Age 0 years 8-10 yws 10 yrs 20 400- 30 1801- Recovered 04- 02- 0.0 100 Time (years)