load_file("init.hoc")
load_file(1, "fig3.ses")
save_window_.vector(nann, &xvec.x[0], &ca_cab[0](.5))
run()

xpanel("Pump")
xlabel("Choose pump rate")
xradiobutton("pump rate = 0", "ke_cab = 0 run()", 1)
xradiobutton("pump rate = 10^-3 cm/sec", "ke_cab = 1e-3*(1+beta_cab) run()")
xradiobutton("pump rate = 10^-2 cm/sec", "ke_cab = 1e-2*(1+beta_cab) run()")
xpanel(550, 2)