load_file("init.hoc") atol_ = 1e-3 cvode.atol(atol_) cvode_active(1) run()