load_file("nrngui.hoc") xopen("fig4A-model.hoc")