function dnei = lci_neuronal(t,dt,input,nei,n0,T_e,T_i,ce,ci,ceu,ciu)

tindex = round(t*1/dt+1,0); 

if nei(1) < n0
    nei(1) = n0;
end
if nei(2) < n0
    nei(2) = n0;
end

dnei=[0;0];

dnei(1) = (1 - nei(1) -  ci*(nei(2)-1) + ceu*input(tindex))/T_e;
dnei(2) = (1 - nei(2) +  ce*(nei(1)-1) + ciu*input(tindex))/T_i;

end