Initial conductance density values:

GD_gsubNa=3p3_gL=fitted.mat - initial conductance values for Figure 4
GD_gsubNa=3p5_gsubK=p7_gM=1p75_gAHP=p5_gL=2_stdev=0.1_300mdls.mat - initial conductance values for Figure 5-10

Istim (fast fluctuating signal):

fast_signal_mu_40_std=10p5_length=1p5sec.mat - the Istim (fast fluctuating signal) used for measuring firing rate

Target ranges:

output_target_values.mat - target ranges for outputs— firing rate, energy efficiency, input resistance, etc.

Simulations:

gridsearch.m
feedback_control.m
feedback_control_dg.m
P_control.m - Runs proportional feedback control without noise
P_control_dg.m - Modified from P_control.m to add noise (dg)
ML_HH_adapt.m - Calculates V based on the input (gi’s) and returns spike timing
ML_HH_adapt_ver2.m - Modified from ML_HH_adapt.m to calculate FR and energy consumption rate
Rheo_ML_HH_JY.m - Modified from ML_HH_adapt.m to calculate rheobase
Vrest_Rinput_ML_HH.m - Modified from ML_HH_adapt.m to calculate Vrest and input resistance
calc_EE.m - Modified from ML_HH_adapt.m to calculate energy efficiency
calc_fmin.m - Modified from ML_HH_adapt.m to calculate minimum frequency
find_inter.m - Returns the intersection of two surfaces
round_JY.m - Round to two decimal places

Figures:

figure3E_H.m - Reproduces Figure3E-H, but can be modified to reproduce contour plots and iso-output1 contour vs. output2 plots in Figure 3A-D
figure4.m - Reproduces Figure 4
figure5.m - Reproduces Figure 5D, which is essentially equal to a combination of Figure 5A-C
figure6-10.m - Reproduces Figure 6, but can be modified to reproduce 3D plots of iso-output surfaces, conductance trajectories and correlation matrices in Figures 7-10