% Full beta-cell model
%initial conditions 
v0=-60;n0=0.0001;s0=.4;
% Conductances in pS; currents in fA; Ca concentrations in uM; time in ms

[t,Y]=ode45(@hw7_2_fun,[0 50000],[v0 n0 s0]);
v=Y(:,1);
n=Y(:,2);
s=Y(:,3);
figure(2)
subplot(2,1,1)
plot(t,v,'-')
xlabel('t')
ylabel('V')
subplot(2,1,2)
plot(t,s)
xlabel('t')
ylabel('s')