%
%
%Solution to Assignment 3, part 2, exercise 1
%
%F. Gabbiani 11/03/08
%
%

v_e = 100; %excitatory reversal potential in mV

%Eqs. constants taken from the appendix of lecture g3
a_0 = 1 + 2/9;
b_0 = 1 + 1/9;
c_0 = 1/9;

%distal excitation
g_ed = 0.01:0.01:1000;

figure;

%
%proximal case
%cycle over the other inhibitory conductance values given by Vu and Krasne
%
g_id_vect = [0.0 0.2 0.5 1.0 2.0 5.0];

for i = 1:length(g_id_vect)
  g_id = g_id_vect(i);

  %proximal membrane potential following the formula given in the 
  %appendix of lecture g3
  v_pi1 = (c_0*g_ed*v_e).*(1./(a_0 + b_0*(g_ed+g_id) + g_ed*g_id));
  semilogx(g_ed,v_pi1);
  hold on;
end;

xlabel('relative excitatory strength');
ylabel('depolarization (mV)');

%
%Repeat, for distal case
%

figure;


for i = 1:length(g_id_vect)
  g_id = g_id_vect(i);
  v_p = (c_0*g_ed*v_e).*(1./(a_0 + b_0*(g_ed+g_id)));
  semilogx(g_ed,v_p);
  hold on;
end;

xlabel('relative excitatory strength');
ylabel('depolarization (mV)');