Added a system.py to contain simulator functions
added a function to process analytical function directly
implemented the ode system and the analytical function for erf(t)