\def{integer a0=random(1,-1)*random(1..5)}
 \def{integer a1=random(1,-1)*random(1..5)}
\def{integer a2=random(1,-1)*random(1..5)}
\def{integer a3=random(1,-1)*random(1..5)}
\def{integer a4=random(-5..5)}
\def{integer a5=random(-5..5)}

\def{integer b1=random(-5..5)}
\def{integer b2=random(1,-1)*random(1..5)}
\def{integer b3=random(-5..5)}
\def{integer b0=random(1,-1)*random(1..5)}

\def{text P= (\a5)*x^5+ (\a4)*x^4+ (\a3)*x^3+(\a2)*x^2+
 (\a1)*x+ (\a0)}
\def{text Q= (\b3)*x^3+ (\b2)*x^2+ (\b1)*x+ (\b0)}


\def{text data= pari((g(P0,Q)=local(n,q,A,B,E,P,EP,LE);
n=poldegree(P0,x); q=poldegree(Q,x);
A=matrix(2*n+1,n+1);P=P0; EP=0;LE=[];
A[1,]=Vec(P); B=concat(vector(n-q),Vec(Q));p=n;
for(i=1, n,if(p>=q,
E=pollead(P,x)/pollead(Q,x); EP=E*x^(p-q)+EP;
LE=concat(LE,[E*x^(p-q)]);
A[2*i,]=-E*concat(vector(n-p),Vec(x^(p-q)*Q));
A[2*i+1,]= A[2*i-1,]+A[2*i,];
P=Pol(A[2*i+1,],x);
p=poldegree(P,x),s=i;break));
[matrix(2*s-1,n+1,i,j,A[i,j]*x^(n-j+1)),s,EP,LE,P]);
R=g(\P,\Q);
 print(poldegree(\P,x)","R[1]","R[2]","R[3]","\Q","R[4] ","\P","divrem(\P,\Q)[2]))}

\def{text n=item(1,\data)}
\def{integer m=2*\n}
\def{text A=wims(declosing item(2,\data))}
\def{text s=item(3,\data)}
\def{text E=htmlmath(item(4,\data))}
\def{text Q=htmlmath(item(5,\data))}
\def{text Liste=wims(declosing item(6,\data))}
\def{text P0=htmlmath(item(7,\data))}
\def{text Rs=htmlmath(item(8,\data))}
\reload{<img src="gifs/doc/etoile.gif" alt="rechargez" width="20" height="20">}
<table width="90%">
<tr align=center> 
\def{text C=row(1,\A)}
<td> P=</td>
\for{j=1 to  \n+1}{
<td><font color=black>
\def{text a= htmlmath(item(\j,\C))} 
 \if{item(\j,\C)=0}{}{\if{ - isin \a}{ \a}{+\a}}

  </font></td>}
 <td colspan="10"> Q=\Q <hr></td>
 </tr>
 <tr align=center>
\for{i =1 to \s-1}{\def{integer iii=2*\i+1}
\def{integer ii=2*\i}
\def{text C=row(\ii,\A)}
\def{text EE=htmlmath(-item(\i, \Liste))}
<tr align=center ><td> \if{ - isin \EE}{<font color=red><b>\EE</b> </font>Q}{
<font color=red><b>\EE</b></font> Q} </td>
\for{j=1 to  \n+1}{
<td ><font color=green>
\def{text a= htmlmath(item(\j,\C))} 
\if{item(\j,\C)=0}{}{
 \if{  - isin \a}{\a}{+\a} <hr>}</font></td>} 
 <td colspan="10"> 
 \if{\i=1}{E=<font color= red><b>\E </b></font>}{} </td>
 </tr>
<tr align=center> 
\def{text C=row(\iii,\A)}
<td> \if{\i<\s -1}{=}{R=}</td>
\for{j=1 to  \n+1}{
<td >\if{\i<\s -1}{<font color=blue>}{<font color=magenta><b>} 
\if{item(\j,\C)=0}{}{
\def{text a= htmlmath(item(\j,\C))} 
 \if{ - isin \a}{\a }{+\a}}
 \if{\i=\s -1}{</b>}</font></td>
 }
 </tr>
}</table>
<center> 
\P0 =<font color=red><b>(\E)</b></font>(\Q)+<font color=magenta><b>
(\Rs)</b></font>
\reload{<img src="gifs/doc/etoile.gif" alt="rechargez" width="20" height="20">}
