tspan = [0 60];
y0 = [0 0.2 0 0.04 0 0 0.1 0.3 0.01 0 0 0 0 0 0 0 0.007 0 0 0];
sol = ode15s(@odefunrevise,tspan,y0);
%sol2 = ode15s(@odetest,tspan,y0);
%y=sol.y;
%x=sol.x;
%plot(log(x),y(1,:));
plot(sol.x,sol.y(7,:));
%plot(log(x),y(3,:));