0.9.7

Nouveautes concernant Xcas.
News about Xcas

Modérateur : xcasadmin

alb
Messages : 1238
Inscription : ven. août 28, 2009 3:34 pm

Re: 0.9.7

Message par alb » lun. avr. 09, 2012 10:54 am

En attendant le retour de bernard parisse je me suis essayé à une animation un peu gourmande en mémoire

Code : Tout sélectionner

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))}:;
On exécute en tapant Animation(1/20) par exemple
RQ Pour bien voir les cas b=-2 et b=2 on peut pendant l'animation faire M>Animation>Pause et cliquer dans le graphique pour avancer pas à pas

Répondre