hold off ;
x = [-2:0.01:2] ;
plot(x,q0(x),'.-b','LineWidth',1.5) ; hold on ;
plot(x,q(x,1),'+-r') ; hold on ;
plot(x,q(x,2),'o-y') ; hold on ;
plot(x,q(x,3),'.-m') ; hold on ;
plot(x,q(x,4),'.-g') ; hold on ;

legend('q0','t=1','t=2','t=3','t=4');
axis([-2,2,-0.1,1.1]);