Page 1 sur 1
un argument particulier de rsolve
Publié : mer. sept. 10, 2014 4:27 pm
par alb
rsolve(u(n+1)=u(n)/(1+u(n)),u(n),u(0)=1) renvoie Unable to isolate rsolve_x0 in seqsolve([rsolve_x0/(rsolve_x0+1),[rsolve_x0,n],[rsolve_x0]]
Est-ce incurable ?
Re: un argument particulier de rsolve
Publié : jeu. sept. 11, 2014 12:32 pm
par parisse
Pour les homographiques, Xcas ne sait pour l'instant faire que s'il y a 2 points fixes distincts r1 et r2, car alors (un-r1)/(un-r2) est geometrique. Il y a une methode simple si la racine est double?
Re: un argument particulier de rsolve
Publié : jeu. sept. 11, 2014 6:49 pm
par xavier
bonjour Bernard
parisse a écrit : Il y a une methode simple si la racine est double?
Oui, la suite 1/(un-r) est arithmétique : voir
ici page 4.
Merci pour tout le travail que tu fais sur Xcas que j'utilise tès souvent avec mes étudiants de cpge ATS.
A bientôt
Re: un argument particulier de rsolve
Publié : jeu. sept. 11, 2014 6:56 pm
par alb
avec en prime les valeurs interdites du premier terme.
Re: un argument particulier de rsolve
Publié : jeu. sept. 11, 2014 7:12 pm
par parisse
peut-etre qu'on peut s'en affranchir, si on accepte d'avoir un terme infini?
Re: un argument particulier de rsolve
Publié : jeu. sept. 11, 2014 7:15 pm
par alb
oui.
Re: un argument particulier de rsolve
Publié : jeu. sept. 11, 2014 8:18 pm
par alb
on peut rêver ?
rsolve(u(n+1)=(4+u(n))/(1+u(n)),u(n),u(0)=-365/182); renvoie [(-1458*(-1/3)^n-2)/(729*(-1/3)^n-1)]
ne pourrait-on pas imaginer un programme
qui afficherait/repondrait:
Attention expression valable pour n strictement inferieure à 6
[(-1458*(-1/3)^n-2)/(729*(-1/3)^n-1)]
ou bien qui repondrait
[6>n,(-1458*(-1/3)^n-2)/(729*(-1/3)^n-1)]
au passage je remarque que solve(729*(-1/3)^n-1,n) renvoie une liste vide.
Edit
c'est peut-etre plus simple de laisser l'utilisateur se debrouiller avec la validite de la formule.
Re: un argument particulier de rsolve
Publié : ven. sept. 12, 2014 5:23 am
par parisse
si j'implémente ce genre de choses, ce sera plutot un warning, pour que le résultat soit facile à exploiter.
Re: un argument particulier de rsolve
Publié : sam. sept. 13, 2014 6:37 pm
par parisse
Ca devrait etre bon sur les versions instables (sans warning...)
Re: un argument particulier de rsolve
Publié : lun. sept. 15, 2014 8:09 pm
par alb
c'est bon
Re: un argument particulier de rsolve
Publié : ven. mai 08, 2015 12:25 pm
par alb
il semble subsister un petit bug dans le cas d'une racine double, peut-etre une erreur sur le premier terme de la suite auxiliaire ?
rsolve(u(n+1)=9/(6-u(n)),u(n),u(0)=1) renvoie [-3/(n-1)] au lieu de [(6*n+3)/(2*n+3)]
si on pose v=1/(u-3) alors v(0)=-1/2 et v(n)=-1/3*n-1/2 et non -1/3*n-2
Re: un argument particulier de rsolve
Publié : ven. mai 08, 2015 1:19 pm
par parisse
Vu, je viens de corriger sur mon source...
Re: un argument particulier de rsolve
Publié : ven. mai 08, 2015 2:40 pm
par alb
c'est encore moi !
est-ce le meme bug pour:
rsolve(u(n+1)=1/(2-u(n)),u(n),u(0)=0) qui devrait retourner [n/(n+1)]
Re: un argument particulier de rsolve
Publié : ven. mai 08, 2015 3:10 pm
par parisse
non, c'est un autre bug dans la decomposition d'expression symbolique en numerateur/denominateur, mais c'est facile a corriger.