options = odeset('RelTol',1e-8,'AbsTol',[1e-8 1e-8 1e-8]);
tspan = [0:0.01:4000];
tspan(length(tspan))