% point.m 
% Sets the initial conditions for the next integration, and 
% resets to initial time to 0.  In two dimensions the coordinates 
% are obtained using the mouse.  In higer dimensions you will be 
% prompted for the coordinates. 

   if (dimension==2)
     disp('Click on initial point.')
     x=ginput(1);
   end

   if (dimension>2)
     disp('Enter coordinate values (don''t use blanks for zero).')
          for i=1:dimension
          x(i)=input(['x(',int2str(i),') :']);
     end
   end

   if (length(x)~=dimension)
     disp('Wrong dimension. Try again (don''t use blanks for zero).')
   end
   t=0;