Le forum de XCAS

Xcas: un logiciel libre de calcul formel
Nous sommes actuellement le Sam Nov 22, 2014 12:44 pm

Heures au format UTC




Publier un nouveau sujet Répondre au sujet  [ 16 messages ]  Aller à la page Précédent  1, 2
Auteur Message
 Sujet du message: Re: 0.9.7
MessagePublié: Lun Avr 09, 2012 10:54 am 
Hors-ligne

Inscrit le: Ven Aoû 28, 2009 3:34 pm
Messages: 848
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


Haut
 Profil  
 
Afficher les messages publiés depuis:  Trier par  
Publier un nouveau sujet Répondre au sujet  [ 16 messages ]  Aller à la page Précédent  1, 2

Heures au format UTC


Qui est en ligne ?

Utilisateurs parcourant actuellement ce forum : Yahoo [Bot] et 1 invité


Vous ne pouvez pas publier de nouveaux sujets dans ce forum
Vous ne pouvez pas répondre aux sujets dans ce forum
Vous ne pouvez pas éditer vos messages dans ce forum
Vous ne pouvez pas supprimer vos messages dans ce forum
Vous ne pouvez pas insérer de pièces jointes dans ce forum

Rechercher pour:
Sauter vers:  
Powered by phpBB® Forum Software © phpBB Group
Traduction réalisée par Maël Soucaze © 2009 phpBB.fr