Try the following:
G_max_chl = 20;
G_max_glu = 30;
Vm = -70:10:0;
tau_rise_In = 0.15;
tau_decay_In = 0.5;
tau_rise_Ex = 0.23;
tau_decay_Ex = 0.7;
tmax = 15;
EGlu = 0;
EChl = -70;
dt = 0.1;
t = 0:dt:tmax;
niter = ceil(tmax/dt);
m_plot_I = 1 - exp(-t / tau_rise_In);
h_plot_I = exp(-t / tau_decay_In);
m_plot_E = 1 - exp(-t / tau_rise_Ex);
h_plot_E = exp(-t / tau_decay_Ex);
Gi = G_max_chl * (1 - exp(-t / tau_rise_In)) .* (exp(-t / tau_decay_In));
Ge = G_max_glu * (1 - exp(-t / tau_rise_Ex)) .* (exp(-t / tau_decay_Ex));
Best Answer