function [y,Y]=go(x,y,t,tmax,xmax,nx,dt,Dr,alpha,Dl)

nt=0;
ntplot=1;
fig1=figure(1);
lin1=line(x,y,'EraseMode','background');
axis([0 xmax -1 1]);
drawnow;

Y=y';

while t < tmax
Yq=(Y.*Y);
Yc=(Y.*Y.*Y);
rhs=Dr*Y-Yc;
Dlc=3*alpha*spdiags([Yq],0,nx,nx);
Dls=Dl+Dlc;
Y=Dls\rhs;
y=Y';

if rem(nt,ntplot) == 0
delete(lin1);
lin1=line(x,y,'EraseMode','background');
drawnow;
end

nt=nt+1;
t=t+dt;
end

