1. simulators.py optimization is 10 times faster due to moving lambdify
out of the loop.
2. Added control variables to simulators.