Le forum de XCAS

Bonjour, j'ai écrit des petits programmes pour dessiner des fractales avec la tortue. Je me suis inspiré d'un programme Python dans lequel on peut faire varier la vitesse de la tortue et la voir dessiner. Je n'ai pas retrouver cette option dans Xcas : j'aurais souhaité montrer à mes élèves comment é...

Juste pour terminer cette discussion, le plus simple, me semble - t - il, c'est encore de demander n puis de calculer le pas. Du coup, on voit bien les résultats obtenus pour x selon que l'on choisit une puissance de 2 ou de 10 et tous les résultats sont cohérents !

Utiliser un while avec des données approximatives c'est dangereux, à cause des erreurs d'arrondi, on risque en effet de faire une itération de plus. Je ne pensais pas que cela poserait problème avec un pas de 0.1. Ensuite plutot que de faire un seq d'appels à Euler, il est beaucoup plus efficace de...

C'est surement la valeur finale de x qui explique le 1er cas. Comme 0.001 est represente en base 2, il y a arrondi/troncature Je me doutais bien qu'il devait y avoir un problème de ce genre... Donc mon conseil c'est d'utiliser une boucle for, d'ajouter un 1e-10 dans le floor pour les erreurs d'arro...

Je viens d'essayer sur mac, je n'observe pas de decalage. Mais ca depend peut-etre de l'architecture. Je pense qu'il serait interessant de renvoyer non pas y mais [x,y] a la fin de la boucle, ca permet de savoir en quel abscisse on a une valeur approchee de y (meme chose pour le 2eme programme avec...

Ce qui m'interpelle le plus : Euldif(f,x,y,xf,h):={ local k,n; n:=floor((xf-x)/h); for(k:=1;k<=n;k++){ y:=y+h*f(x,y); x:=x+h; } retourne evalf(y,6); }:; Ensuite f(x,y):=sin(x,y) Puis seq(Euldif(f,0,3,k,0.1),k=0..2,0.1) donne pour résultat (3.0,3.0,3.02955,3.02955,3.0865,3.16642,3.26183,3.36165,3.451...

Bonjour, je souhaiterais résoudre numériquement une équation différentielle y'=f(x,y) avec les méthodes d'Euler et du point milieu en xf avec (x0,y0) les conditions initiales. Dans un 1er temps, j'ai fixé le nombre d'opération à n, et j'ai écrit une boucle for(k:=1;k<=n;k++) suivi des calculs avec h...

Bonjour, oui effectivement en créant ainsi le point M on peut le voir évoluer dans l'espace et avoir ainsi des informations déjà intéressantes et plus visuelles que la courbe elle même parfois. Je vous remercie de votre aide. J'en profite pour poser une autre question qui n'a rien à voir : pour l'ét...

Bonjour, je me suis mal exprimé mais c'est l'idée que j'avais. J'ai tracé une surface G, j'ai créé deux curseurs u et v puis j'ai tapé M:=element(G,[u,v]) ou M:=element(G,u,v) mais le point M n'est pas sur la surface mais au centre d'un repère 2D et il ne bouge pas lorsque u et v varient. Il y a peu...

Bonjour, Je voudrais savoir si c'est possible de faire apparaître un point mobile sur une surface G donnée par une fonction à deux variables ? J'ai essayé en créant deux paramètres u et v puis en utilisant la fonction élément, mais je perds le graphique 3D et le point apparaît seul sur un graphique ...

Merci beaucoup pour toutes ces informations. Il faut vraiment que je prenne l'habitude de résoudre des systèmes avec le noyau de la matrice effectivement. Quant à la validité de l'algorithme sur des grandes dimensions, je ne doute pas qu'il explose !
François

Bonjour, j'ai écrit un algorithme naïf permettant de déterminer le polynôme minimal d'une matrice. L'idée est de résoudre les n^2 équations à k-1 inconnues suivantes : X[0]*A^0+X[1]*A^1+...+X[k-1]*A^{k-1}+A^k en faisant varier k entre 1 et n^2. Je teste mon programme sur I_2 qui me renvoie -1+x, lor...

[quote="parisse"]C'est lie a la valeur de epsilon par defaut. Si on met epsilon a 1e-12 alors on a la valeur. Je croyais que la valeur d'espilon interviendrait pour un grand nombre d'itération au contraire, lorsque Xcas calcule avec des valeurs proches de 0. Existe - t - il un site ou un livre qui v...

Bonjour Est-il normal que aire(1/(1+x),x=0..1,N,rombergt) n'arrive pas à renvoyer une valeur approchée de ln(2) pour N inférieur à 7 ? J'ai le message suivant affiché " Impossible de touver l'intégrale numérique avec Romberg, on renvoie la dernière ligne calculée des approximations puis la suite des...

Bonjour, Dans l'aide en ligne, les options d'area sont orthographiées rombergt et romberm mais ne sont pas reconnues par Xcas qui propose romberm et rombert. Du coup, taper area(x^2,x=0..1,3,romberm) renvoie quelque chose d'assez improbable mais c'est reconnu (l'option est en rouge) alors qu'il faut...