%
% ssEI.m
%
% Gabbiani & Cox, Mathematics for Neuroscientists
%
% steady-state depolarization as a functino of normalized conductance for 
% one and two excitatory synapses, as well as their linear sum where Vss is
% the syanptic input and ve = 0
% 
% usage: ssEI
%

function ssEI

figure(1)
ce = 0:.1:15;     % voltage range
VCl = -68;

Vss = VCl./(1+ce);      % synaptic input

plot(ce,Vss,'k');
xlabel('c_e','fontsize',14,'fontname','helvetica medium')
ylabel('V_{ss} (mV)','fontsize',14,'fontname','helvetica medium')
box off
set(gca,'fontsize',14)
text(-2, -6, '(A)','fontsize',18,'fontweight','bold')
print -depsc ssEIA

figure(2)
linecolor = {'k' 'r' 'r'};
linetype = {'-' '--' '-'};

for k = 1:3,

    ci = (k-1)*25;

    Vss = ((1 + ci)*VCl)./(1 + ce + ci);      

    plot(ce,Vss,'color',linecolor{k},'linestyle',linetype{k});
    hold on

end

hold off
box off
legend('c_i = 0','c_i = 25', 'c_i = 50','location','best')
xlabel('c_e','fontsize',14,'fontname','helvetica medium')
ylabel('V_{ss} (mV)','fontsize',14,'fontname','helvetica medium')
set(gca,'fontsize',14)
text(-2, -6, '(B)','fontsize',18,'fontweight','bold')
print -depsc ssEIB