### Python Package
import cr_mech_coli as crm
# Define Settings
config = crm.Configuration()
ag_set = crm.AgentSettings()
# Generate positions and agents
positions = crm.generate_positions_old(4, ag_set,config)
agents = [crm.RodAgent(
pos=p,
vel=p*0.0,
**ag_set.to_rod_agent_dict()
) for p in positions
]
# Run simulation and plot individual results
cells = crm.run_simulation_with_agents(agents, config)
crm.store_all_images(cells, config.domain_size)