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]
rsolve
Modérateur : xcasadmin
Re: rsolve
En effet, petit bug
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.
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);