% v-n fast subsystem for phase plane analysis and bifurcation diagrams
%initial conditions 
v0=-30;n0=0.1;
% Conductances in pS; currents in fA; Ca concentrations in uM; time in ms

[t,Y]=ode45(@hw7_1_fun,[0 4000],[v0 n0]);
v=Y(:,1);
n=Y(:,2);
figure(1)
plot(v,n,'-')
axis([-80 -10 -.05 .15])