Please see run.jl
in order to run the simulations shown in the Overleaf document. Here's an overview of the code:
NeuralPopSoma{Area,Type}
is a method that creates a point neuron model in regionArea
(can be eitherCortex
orThalamus
at this point) of typeType
(Inhibitory/Excitatory
)- Two
NeuralPopSoma
can be grouped into oneNeuralArea
, for which the functionPar
computes all impermeants and leak conductances - The method
BioNM
constructs the neural mass that couples multiple areas, viaBioNM.conn
. Note thatBioNM.areas
is a vector of typeNeuralArea
. The first element of this vector receives stimulation - Change the hyperparameters using the method
HyperParam
solve(nm::BioNM)
,plot(nm::BioNM)
,plot_syn(nm::BioNM,1)
solves the system, plot ion dynamics and plot synaptic currents respectively.