global a b c r
a= 1; b= 0.5; c= -4/9;

options= odeset('RelTol',1e-8,'AbsTol',1e-8);

% first set of parameters
r= 1;
[t1,p1]= ode45('ode',[0 6*pi],[-0.563 0.869*i],options);

% second set of parameters
r= 0.6;
[t2,p2]= ode45('ode',[0 6*pi],[-0.944 0.658*i],options);

clf;
plot(p1(:,1),'b-');
hold on;
plot(p2(:,1),'r--');
hold off;

