Computational model and files to recreate the results from the paper "A Computational Model of the Mammalian External Tufted Cell" https://doi.org/10.1016/j.jtbi.2018.10.003
Author: Ryan Viertel
data = ET(input,sampling_rate);
input: input vector, if no input then just a vector of zeros
sampling_rate: rate at which the input vector should be sampled. 1000 for milisecond
returns the following struct:
data.T - time vector data.X - ODE variables at each time step
data.events - list of spike events
data.which - event type
data.current - system currents
input = zeros(1,5000);
data = ET(input,1000);
plot(data.T,data.X(:,1))
The code used to sample the model throughout parameter space to determine stability and investigate the effect of model parameters on model output is found in the ME-PCM directory
The ODE file to recreate the bifurcation diagram is found in the xpp directory