Le forum de XCAS

parisse a écrit :j'ai aussi fait le laplacien du coup. A peine teste (j'ai juste verifie curl(grad()) et div(curl())==0 en spheriques), donc pas encore commente.

Merci tout plein!

parisse a écrit :Du coup il va falloir aussi faire divergence et curl j'imagine?

Ben... ce serait bien, merci!
odv

parisse a écrit :Il n'y a rien pour l'instant. Je vais ajouter grad(G,[r,t,f],coords=sphere);

Merci
odv

Bonjour,
Je ne trouve pas comment faire un gradient en coordonnées sphérique.
En maple, on faisait:
grad(G,[r,t,f],coords='spherical');
C'est pour avoir les coefficients 1/r et 1/r/sin(theta) devant les coeff theta et phi, bien sûr.
Merci,
odv

parisse a écrit :A mon avis le plus simple c'est d'introduire un facteur d'echelle, appelons-le h, faire un subst par exemple avec 2 variable subst(expression,[N1,N2],[h*n1,h*n2]) puis travailler sur l'expression en puissances de h.

OK, merci beaucoup!
odv

Oui, avec la syntaxe suivante de series, series(1/(1+x+y),[x,y],[0,0],5) Je vais l'ajouter dans l'aide. OK, merci... Si je peux abuser, je redirige un peu. J'ai une très grande expression, avec 8 coefficients différents, notons les Ni. Chaque terme en contient un certain nombre au numérateur et au ...

Bonjour,
Y a-t-il moyen de faire un développement limité à plusieurs variable, par exemple
1/(1+x+y)=1-x-y...
Dans les anciens maples, on écrivait qqcho comme mtaylor(1/(1+x+y),{x,y},2)...

Merci,
odv

je vous conseille plutot de multiplier par la puissance negative la plus grande et prendre le coefficient shifte de cette puissance. Par exemple a:=series(sqrt(x^2+1)-x,x=inf,5,polynom); coeff(a*x^3,2) donne le coefficient de 1/x de a. Je vais reflechir a etendre coeff pour accepter un argument neg...

Bonjour, J'ai une expression est une somme de fraction, et je souhaiterais récupérer tous les termes qui sont en, par exemple, 1/N1^2. J'ai essayé coeff et collect, et cela ne fonctionne pas... me semble-t-il. J'ai une solution batarde qui consiste à substituer 1/N1 par iN1, et à travailler avec cel...

parisse a écrit :En effet, il y a un petit bug en mode maple car l'ordre des arguments de map n'est pas le meme qu'en mode xcas, je le corrige.
Sinon, la commande alternative qui fonctionne serait plutot
xd:='+'(seq(simplify(op(k,x)),k=1..nops(x)))
ou
xd:=sum([seq(simplify(op(k,x)),k=1..nops(x))])

Ca marche!
Merci

Bonjour à tous, et désolé si la réponse existe quelque part... Je l'ai cherchée et pas trouvée. Je souhaite simplifier séparément une somme de termes. Dans les anciens "maple", on faisait cela avec map(simplify, expression) Visiblement, ceci ne fonctionne pas sur Xcas. J'ai essayé aussi de faire une...