Download Pid: Toolbox !!hot!!

T = ct.feedback(C * G) t, y = ct.step_response(T, time_vector=np.linspace(0, 40, 1000))