introduce some basic parasitic simulation using a numeric integration library (odeint) and a matrix library (Eigen)