Please answer in python

Consider a system of a parent nucleus, P, and a daughter nucleus D, both radioactive. The coupled equations which describe their decays are as follows: dNp = -kpNp(t), dt dND - kpNp(t) kpNp(t). dt = = Write a program to numerically solve the above coupled equations and plot the time dependence of Np and Np for kp 0.5s-1 and kp 5.0s-1, 0.55-1, and 0.05s-1. Assume you start out with the number of parent nuclei equal to 100 and no daughter nuclei. The interval for t is 0