# define a function containing the derivatives, i.e. the ODE def radioleak(x,t,params): r,beta = params derivs = r - beta * x return derivs #