next up previous contents index
suivant: Application à monter: Activité précédent: L'algorithme   Table des matières   Index

Traduction Xcas

On écrit la fonction decomp qui va décomposer selon l'algorithme la fraction frac. Cette fonction va renvoyer la liste lres égale à $ \tt [Q_0,..Q_n]$ avec $ \tt0<Q_0<,..<Q_n$ et $ \tt frac=1/Q_0+..+1/Q_n$.
Attention frac=b/a et donc fxnd(frac)=fxnd(b/a)=[b,a].
decomp(frac):={
local a,b,l,q,r,lres;
l:=fxnd(frac);
b:=l[0];
a:=l[1];
q:=iquo(a,b);
r:=irem(a,b);
lres:=[];
while (r!=0) {
lres:= concat(lres, q+1);
b:=b-r;
a:=a*(q+1);
q:=iquo(a,b);
r:=irem(a,b);
}
lres:=concat(lres,q);
return lres;
}


Documentation de giac écrite par Renée De Graeve