
%UNTITLED Summary of this function goes here
%   Hodgkin-Huxley
%I=1

function dydt=hoha(t,u)

I=0
dydt=[I-120*u(3).^3*u(4)*(u(1)-115)-36*u(2).^4*(u(1)+12)-0.3*(u(1)-10.6);0.01*(10-u(1))./(exp((10-u(1))./10)-1)*(1-u(2))-(1/8)*exp(-u(1)./80)*u(2);0.1*(25-u(1))./(exp((25-u(1))./10)-1)*(1-u(3))-(4)*exp(-u(1)./18)*u(3);0.07*exp(-u(1)./20)*(1-u(4))-1/(exp((30-u(1))./10)+1)*u(4)];