clear,clc,clf
global A B
A=1;B=1.8;
x0=[1 0]';
[t,y]=ode45('brusselator',[0 100],x0);
figure(1),plot(t,y,'linewidth',2)
%figure(2),plot(y(:,1),y(:,2),'linewidth',2)
N=length(t);
figure(2)
p1=animatedline;p1.LineWidth=1;p1.Color='blue';
for i=1:N
 addpoints(p1,y(i,1),y(i,2));
 drawnow limitrate
 pause(0.02)
end