function diffeq_osc

close all
tic
%initial conditions

%global y

y0(1)=-56;
y0(2)=-66;
y0(3)=.9;
y0(4)=.2;

tspan=[0 1000];

%current injection
%i_electrode4uA/cm^2  

%use vector of values for step injections

%options -- set carefully  -- consider odephas2 as output
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-4 1e-4],'stats','on','outputfcn','odeplot');

%dy=hh_solve(t,y,mem_capacitance,g_l,e_l,g_k,e_k,g_na,e_na,i_electrode,cell_area)

% w_vec=[.01 1];

% amp_ss_vec=[];

% for i=1:length(w_vec)
    [t,y] = ode15s(@osc_solve,tspan,y0);
%     amp_ss_vec(i)=max(v(ceil(.8*length(v)):end));
% end

plot(y,'-o')
size(y)
y(1),y(2),y(3),y(4)
% figure,plot(amp_ss_vec)
toc