function y = hhsolve(t, state)

%This code solves the Hodgkin-Huxley equations.  The input
%is a vector of time instants <t> = [t0 t1 t2 t3 ...]  and a
%columnx vector <state> = [v0;m0;n0;h0].  The output is a
%matrix whose rows are the solutions at the time instants
%given in <t>.

%%% These default parameters are Hodgkin and Huxley's original values using
%%% today's definition of membrane potential.

[t,y] = ode45(@hhfield, t, state);
