function dy=hw7_1_fun(t,Y)
dy=zeros(2,1); 
v=Y(1);
n=Y(2);

s=-0.1;

%Ik
gk=2700.0;vn=-16;sn=5.6;vk=-75;taun=20;

% Is
gs=200;vs=-52;ss=5;

% Ica
gca=1000;vca=25;vm=-20;sm=12;

%Ikatp
gkatp=120;
cm=5300;

%Functions
ninf=1./(1+exp((vn-v)/sn));
minf=1./(1+exp((vm-v)/sm));
sinf=1./(1+exp((vs-v)/ss));
ica=gca*minf.*(v-vca);

ik=gk*n.*(v-vk);
is=gs*s.*(v-vk);
ikatp=gkatp*(v-vk);

% Equations
dy(1)=(-ica-ik-is-ikatp)/cm;
dy(2)=(ninf-n)/taun;