%solve ODE 
[t,x] = ode23s('p426singlefun',[0;10000],[0.0]);