%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% MorrisLecar
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Name        MorrisLecar
Coordinates V,n
Parameters  phi,gCa,V3,V4,Iapp

CM=5                                                
gK=8                                                
gL=2                                                
VCa=120                                             
VK=-80                                              
VL=-60                                              
V1=-1.2                                             
V2=18                                               
minf=(1+tanh((V-V1)/V2))/2                          
taun=1/cosh((V-V3)/(2*V4))                          
ninf=(1+tanh((V-V3)/V4))/2                          
V'=(-gCa*minf*(V-VCa)-gK*n*(V-VK)-gL*(V-VL)+Iapp)/CM
n'=phi*(ninf-n)/taun                                

parametry pro bif. diagram:
phi=1/15
gCa=4
V4=17.4

V3=12 nebo i jinak
Iapp= ... ? zkuste 0-120

bif. analýza podle Tsumoto: Bifurcations in Morris-Lecar neuron model
web typ I a II: https://neuronaldynamics.epfl.ch/online/Ch4.S4.html