Le forum de XCAS

D(U):=unapply(U(k+1)-U(k),k) Effectivement, j'avais mal utilisée unapply Attention toutefois, D@@k est très couteux à exécuter, il vaudra mieux faire une boucle et calculer successivement les D(D(...(D())))) dans une boucle que de faire une sum sur des D@@k je confirme ! Merci de la rapidité de vot...

Bonjour, Je souhaiterai à partir d'une série alternée de terme général (-1)^k U(k) l'accélérer avec la transformée d'Euler. Pour cela j'avais pensé définir l'opérateur sur U : D(U)(k):=U(k)-U(k+1) pour pour pouvoir le composer avec @@ et faire calculer à Xcas sum((D@@k)(U)(0)/2^(k+1),k=0..1000) par ...

J'ai accès à la nouvelle page du Capes. Merci bien

Pourtant quand je souhaite aller sur cette partie du forum, j'ai le message suivant :

Vous n’avez pas les permissions appropriées afin de lire les sujets de ce forum.

Peut être me suis-je mal inscrit ?

je me suis mal exprimé, c'est un bug de xcas, vous ne pouvez rien faire qu'attendre une mise à jour (la semaine prochaine). Cela me rassure ! En tout cas je voulais vous remercier pour votre célérité. J'en profite pour vous faire part de mon étonnement concernant la recherche d'information sur le f...

il s'agit d'un bug, dont le correctif est: diff gen.cc gen.cc~ 6898c6898 < return is_zero(a-1); --- > return evalf_double(a,0,context0)._DOUBLE_val==1; 6927c6927 < return is_zero(a+1); --- > return evalf_double(a,0,0)._DOUBLE_val==-1; Que dois-je faire pour corriger le bug ? Comment utiliser le cod...

Bonjour, je cherche à approximer Pi avec un grand nombre de décimal (mettons 500). J'ai donc écrit le calcul de cette suite que je souhaite accélérer ensuite U(n):={ local (c:=0.),(u:=2.),k; for(k:=1;k<=n;k++){ c:=sqrt((c+1)/2);u:=u/c; } return u }:; A priori j'ai configuré avec Digits:= 700 pour av...

A mon avis, le plus simple dans ce cas c'est de chercher le noyau de la matrice dont les colonnes sont les vecteurs de F1 puis de F2 et en derniere colonne le vecteur que l'on cherche a decomposer, puis tester si un des vecteurs de la base renvoyee est de la forme [x y z ... -1], s'il n'y en a pas,...

Bonjour, J'écris un algorithme qui à partir d'une base de F1, d'une base de F2 et d'un élément V de F1+F2 donne des coordonnées de V sur F1 et F2. Si le nombre total de vecteurs des bases de F1 et F2 est inférieur ou égal à N dimension de l'espace vectoriel, en utilsant la fonction rref, je trouve d...

La notation (,,) meme avec des parentheses cree une sequence, pas une liste. Il faut utiliser des crochets comme delimiteur: B:=[E1,...,En], vous pourrez ensuite passer B en 2eme argument sans que le 1er argument et les elements de B ne soient concatenes. Je viens d'essayer et effectivement tout ro...

Bonjour, je souhaite décomposer une matrice P sur une base de matrices donnée sous forme d'une liste B:=(E1,...,En) J'ai écrit une procédure "décompose(P,B)" ayant donc P et B comme arguments et qui en premier calcule la taille de P et la dimension de B. Le problème est que Xcas n'y arrive pas car i...

Bonjour, Merci de votre réponse Je viens d'écrire le programme avec X et k en variable globale : 1) j'obtiens le message suivant : "[[[1,0,0],[0,1,0],[0,0,1]],[[1,-1,-1],[-1,1,-1],[-1,-1,1]],[[1,-1,0],[0,1,-1],[-1,0,1]],[[1,0,-1],[-1,1,0],[0,-1,1]]] size() != [P,B] Erreur: Valeur Argument Incorrecte...

Bonjour, J'ai écrit ce programme et lorsque je le teste ligne de commande par ligne de commande (en variable globale donc) tout fonctionne correctement, et j'obtiens les coordonnées d'une matrice P sur une base de matrice B : decompose(P,B):={ local I,Q,n,d,sys,X,k,l,s; n:=coldim(P);d:=size(B); Q:=e...

parisse a écrit :Pas que je sache, l'algorithme de la décomposition LU (et toute autre variation du pivot de Gauss) fait pratiquement toujours une permutation de ligne, sinon on risque en numérique de faire des erreurs plus grandes et en exact des calculs plus compliqués.

Merci bien

Cela dépend du type de coefficients de la matrice. S'ils sont approchés, le pivot est choisi de module maximal dans la colonne. S'ils sont exacts, c'est le pivot le plus "simple" qui est choisi, pour des entiers cela veut dire le plus petit non nul, pour des polynômes le plus petit degré possible, ...