Sistema de trazado de ecuaciones diferenciales en Python

Acabo de comenzar a usar Python para el dibujo científico para trazar una solución numérica de ecuaciones diferenciales. Se como usarscipy.odeint para resolver y trazar ecuaciones diferenciales simples, pero no tengo idea de los sistemas de ecuaciones diferenciales. ¿Cómo puedo trazar el siguiente sistema acoplado?

N' = a * N - (C/(1+C)) * b * N
C' = (C/(1+C)) * N - C + 1

a = 4
b = 7
N(0) = 100
C(0) = 5

Respuestas a la pregunta(1)

Su respuesta a la pregunta