%solve ODE 
[w,x] = ode23s('cdp4fpd2fun',[0 100],[0;1]); 
plot(w,x(:,1),w,x(:,2));
xlabel('W');
ylabel('X or Y')
legend('X','Y'); 