%solve ODE 
[w,x] = ode23s('cdp4ffun',[0 100],[0]); 
