rsolve

Bugs

Modérateur : xcasadmin

alb
Messages : 1229
Inscription : ven. août 28, 2009 3:34 pm

rsolve

Message par alb » sam. nov. 02, 2013 5:19 pm

salut,
sauf erreur de ma part cette commande:
a:=-1/sqrt(6);
rsolve(u(n+2)=(a*sqrt(3)-sqrt(2)/2)*u(n+1)+a*sqrt(6)/2*u(n),u(n),u(0)=sqrt(2)+sqrt(3),u(1)=3a-1);
devrait renvoyer:
[(sqrt(2)+sqrt(3))*((-(sqrt(2)))/2)^n]
mais on obtient Type Argument Incorrect

j'ai cru que c'était lie au fait que le polynome caracteristique avait une racine double et que dans la solution generale (a+b*n)*r^n on a b=0
mais un exemple similaire marche bien:
rsolve(u(n+2)=4*u(n+1)-4*u(n),u(n),u(0)=3,u(1)=6) renvoie [2^n*3]

parisse
Messages : 5084
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Re: rsolve

Message par parisse » dim. nov. 03, 2013 8:14 am

En effet, petit bug

Code : Tout sélectionner

diff maple.cc maple.cc~
1867c1867
<       return pow(l0,n,contextptr)*sum(r2e(P,v,contextptr),n,remains,contextptr)/r2e(l,v1,contextptr)/r2e(P0d,v,contextptr);
---
>       return pow(l0,n,contextptr)*sum(r2e(P,v,contextptr),n,remains,contextptr)/l/r2e(P0d,v,contextptr);
Par contre, ca renvoie rootof([[1,-1,-19,-25],[1,0,-22,-48,-23]])/4 et pas sqrt(2)+sqrt(3) en facteur, il manque un code de detection de rootof -> somme de racines carrees.

Répondre