Le forum de XCAS
http://pcm1.e.ujf-grenoble.fr/XCAS/

0.9.7
http://pcm1.e.ujf-grenoble.fr/XCAS/viewtopic.php?f=8&t=1008
Page 2 sur 2

Auteur:  alb [ Lun Avr 09, 2012 10:54 am ]
Sujet du message:  Re: 0.9.7

En attendant le retour de bernard parisse je me suis essayé à une animation un peu gourmande en mémoire
Code:
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

Page 2 sur 2 Heures au format UTC
Powered by phpBB® Forum Software © phpBB Group
http://www.phpbb.com/