% 
% At a Matlab prompt, type the file name of this script to run the program.
%
zeta=.1;eps=.4;b1=.3;b2=.1;d1=.2;d2=.62;s1=0.25;s2=0.25;% Set the parameter
values. tend=200;
% Set the range for t.
%initial=[y(end,1),y(end,2),y(end,3)]; 
initial=[8.8918461349712519e-001  5.2908259272560820e-001  4.4607645052032174e-001]; %for 'blue'.
%initial=[8.8918461349712519e-001  5.29082592e-001  4.4607645052032174e-001];
options=odeset('e','off','RelTol',1e-8,'AbsTol',[1e-16;1e-16;1e-16]);
%
[t,y]=ode15s('fdchain',[0 tend],initial,options,zeta,eps,b1,b2,d1,d2);
%
% Plot the graphs.
%
figure;
plot3(y(:,1),y(:,3),y(:,2));
view(50,20);
text(.7,-.2,-.1,'x');
zlabel('y');
text(1.2,.4,-.1,'z');
% 
