Index exceeds array elements timei

function membPoten = simulateAdEx(params)

%membPoten = simulateAdEx(params)

% Simulate a neuron according to the adaptive exponential model.

%

% INPUT: params - a structure containing simulation parameters,

% with the following required fields:

%

% params.C : capacitance

% params.gl : leak conductance

% params.El : leak voltage

% params.vt : resting membrane potential

% params.delT : spike width

% params.a : resonantor/integrator constant

% params.tauW : adaptation decay time

% params.b : adaptation jump (for updating w)

% params.vreset : voltage reset value

% params.simDur : simulation time in ms

% params.dt : integration time step in ms

% params.input : input current in pA

% params.EC50 : the drug concentration that produces a 50% maximal response

% params.Max : maximal inhibition

% params.H : Hill coefficient

% params.x : Drug concentration

%

% OUTPUT: membranePotential - A vector of the neuron's membrane

% potential in mV.

%

%

% Equations based on:

% https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2798047/

% Read more about this model:

% http://www.scholarpedia.org/article/Adaptive_exponential_integrate-and-fire_model

%

%

% vector to store membrane potential

membPoten = [params.vt; zeros(params.simDur/params.dt,1)];

v = params.vt;

w = 0;

vpeak = 0; % spike threshold

% simulation in loop over time

for timei = 1:params.simDur/params.dt-1

% set input for this time point

Iin = params.input(timei);

% integration

v = v + params.dt*(-params.gl*(v-params.El) + params.gl*params.delT*exp((v-params.vt)/params.delT) + Iin - w)/params.C;

% note: update w by pre-updated v

w = w + params.dt*(params.a*(membPoten(timei) - params.El) - w)/params.tauW;

% update and store variables

if v>=vpeak

v = params.vreset;

w = w + params.b;

else

% no AP, just store v

membPoten(timei+1) = v;

end

end

end % end function

Hello, I am getting an error in line 47 for index exceeds the number of array elements (1). Please help with this problem.

[Walter Roberson]

As outside observers, we have no reason to believe that params.input is any particular size or datatype or range. You did not attach the inputs for the function.

My guess is that params.input is a scalar but the code expects there to be one input for each time step.