att:=affichage=[epaisseur_ligne_3,epaisseur_ligne_3];

Raccord(a,b):={
  local sol,s,L,k,j,eps;
  sol:=solve(x^3+a*x+b);
  s:=dim(sol);
  L:=NULL;
  pour k de 0 jusque s-1 faire
    pour j de 1 jusque 3 faire
      eps:=10^-j;
      L:=L,plot([sqrt(x^3+a*x+b),-sqrt(x^3+a*x+b)],x=sol[k]-eps..sol[k]+eps,att)
    fpour
  fpour
  retourne L
}:;

Cubique(a,b):={
  ClrIO;
  legende([50,50],"a="+round(a,2)),
  legende([50,80],"b="+round(b,2)),
  legende([50,110],"delta=4a³+27b²="+(4a^3+27*b^2)),
  plot([sqrt(x^3+a*x+b),-sqrt(x^3+a*x+b)],x=-10..10,att),
  Raccord(a,b);
}
:;

Animation(p):={local b;animation(seq(Cubique(-3,b),b,-10,10,p))}:;