load xyin;
%y0 = [.1 .2 .2]; 
y0=xyin;
options = odeset('RelTol',.0000001,'AbsTol',.000000001)
%
delta = .5;                             % set parameter
tspan = [0:.5:5000];                    % set time interval for calculation
st = size(tspan)                        % number of time steps
p = .01.*randn(1,st(2));                % random perturbation
%
% Call the differential equation solver
[tt,yy] = ode45('dwm',tspan,y0,options,delta,p);    % solve the differential equations over the time interval
%
plot3(yy(:,1),yy(:,2),yy(:,3));
%
% save end state for possible restart
sy = size(yy);
xy = yy(sy(1),:)';
save xyout xy -ascii