%    dx/ds for particle paths
     function xdot=EX413(t,x)
     xdot=zeros(2,1);
     xdot(1)=x(1)/(1+t);
     xdot(2)=x(2);
