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