Question: fix the collect

restart;
T := mu+lambda*H(xi)+(v-1)*H(xi)^2;
                                                2
               mu + lambda H(xi) + (v - 1) H(xi) 
u[0] := a[0]+a[1]*(d+H(xi))+a[2]/(d+H(xi))+a[3]*(d+H(xi))^2+a[4]/(d+H(xi))^2;
                                a[2]                      2
    a[0] + a[1] (d + H(xi)) + --------- + a[3] (d + H(xi)) 
                              d + H(xi)                    

             a[4]    
       + ------------
                    2
         (d + H(xi)) 
diff(u[0], xi);
                            / d        \
                       a[2] |---- H(xi)|
        / d        \        \ dxi      /
   a[1] |---- H(xi)| - -----------------
        \ dxi      /                2   
                         (d + H(xi))    

                                                 / d        \
                                          2 a[4] |---- H(xi)|
                           / d        \          \ dxi      /
      + 2 a[3] (d + H(xi)) |---- H(xi)| - -------------------
                           \ dxi      /                 3    
                                             (d + H(xi))     
collect(%, diff(H(xi), xi));
/           a[2]                               2 a[4]   \ / d       
|a[1] - ------------ + 2 a[3] (d + H(xi)) - ------------| |---- H(xi
|                  2                                   3| \ dxi     
\       (d + H(xi))                         (d + H(xi)) /           

   \
  )|
   /
d[1] := (a[1]-a[2]/(d+H(xi))^2+2*a[3]*(d+H(xi))-2*a[4]/(d+H(xi))^3)*T;
 /           a[2]                               2 a[4]   \ /  
 |a[1] - ------------ + 2 a[3] (d + H(xi)) - ------------| \mu
 |                  2                                   3|    
 \       (d + H(xi))                         (d + H(xi)) /    

                                  2\
    + lambda H(xi) + (v - 1) H(xi) /
diff(d[1], xi);
/       / d        \                                / d        \\ 
|2 a[2] |---- H(xi)|                         6 a[4] |---- H(xi)|| 
|       \ dxi      /          / d        \          \ dxi      /| 
|------------------- + 2 a[3] |---- H(xi)| + -------------------| 
|              3              \ dxi      /                 4    | 
\   (d + H(xi))                                 (d + H(xi))     / 

  /                                 2\   /           a[2]    
  \mu + lambda H(xi) + (v - 1) H(xi) / + |a[1] - ------------
                                         |                  2
                                         \       (d + H(xi)) 

                             2 a[4]   \ /       / d        \
   + 2 a[3] (d + H(xi)) - ------------| |lambda |---- H(xi)|
                                     3| \       \ dxi      /
                          (d + H(xi)) /                     

                     / d        \\
   + 2 (v - 1) H(xi) |---- H(xi)||
                     \ dxi      //
collect(%, diff(H(xi), xi));
//   2 a[2]                  6 a[4]   \ /                 
||------------ + 2 a[3] + ------------| \mu + lambda H(xi)
||           3                       4|                   
\\(d + H(xi))             (d + H(xi)) /                   

                  2\   /           a[2]                         
   + (v - 1) H(xi) / + |a[1] - ------------ + 2 a[3] (d + H(xi))
                       |                  2                     
                       \       (d + H(xi))                      

        2 a[4]   \                           \ / d        \
   - ------------| (lambda + 2 (v - 1) H(xi))| |---- H(xi)|
                3|                           | \ dxi      /
     (d + H(xi)) /                           /             
d[2] := ((2*a[2]/(d+H(xi))^3+2*a[3]+6*a[4]/(d+H(xi))^4)*(mu+lambda*H(xi)+(v-1)*H(xi)^2)+(a[1]-a[2]/(d+H(xi))^2+2*a[3]*(d+H(xi))-2*a[4]/(d+H(xi))^3)*(lambda+(2*(v-1))*H(xi)))*T;
//   2 a[2]                  6 a[4]   \ /                 
||------------ + 2 a[3] + ------------| \mu + lambda H(xi)
||           3                       4|                   
\\(d + H(xi))             (d + H(xi)) /                   

                  2\   /           a[2]                         
   + (v - 1) H(xi) / + |a[1] - ------------ + 2 a[3] (d + H(xi))
                       |                  2                     
                       \       (d + H(xi))                      

        2 a[4]   \                           \ /                 
   - ------------| (lambda + 2 (v - 1) H(xi))| \mu + lambda H(xi)
                3|                           |                   
     (d + H(xi)) /                           /                   

                  2\
   + (v - 1) H(xi) /

eq := (2*k*k)*w*beta*d[2]-(2*alpha*k*k)*d[1]-2*w*u[0]+k*u[0]*u[0];
   2        //   2 a[2]                  6 a[4]   \ /  
2 k  w beta ||------------ + 2 a[3] + ------------| \mu
            ||           3                       4|    
            \\(d + H(xi))             (d + H(xi)) /    

                                 2\   /           a[2]    
   + lambda H(xi) + (v - 1) H(xi) / + |a[1] - ------------
                                      |                  2
                                      \       (d + H(xi)) 

                             2 a[4]   \                          
   + 2 a[3] (d + H(xi)) - ------------| (lambda + 2 (v - 1) H(xi)
                                     3|                          
                          (d + H(xi)) /                          

   \ /                                 2\            2 /    
  )| \mu + lambda H(xi) + (v - 1) H(xi) / - 2 alpha k  |a[1]
   |                                                   |    
   /                                                   \    

         a[2]                               2 a[4]   \ /  
   - ------------ + 2 a[3] (d + H(xi)) - ------------| \mu
                2                                   3|    
     (d + H(xi))                         (d + H(xi)) /    

                                 2\       /    
   + lambda H(xi) + (v - 1) H(xi) / - 2 w |a[0]
                                          |    
                                          \    

                          a[2]                      2
   + a[1] (d + H(xi)) + --------- + a[3] (d + H(xi)) 
                        d + H(xi)                    

         a[4]    \     /                            a[2]   
   + ------------| + k |a[0] + a[1] (d + H(xi)) + ---------
                2|     |                          d + H(xi)
     (d + H(xi)) /     \                                   

                     2       a[4]    \  
   + a[3] (d + H(xi))  + ------------|^2
                                    2|  
                         (d + H(xi)) /  
value(%);
   2        //   2 a[2]                  6 a[4]   \ /  
2 k  w beta ||------------ + 2 a[3] + ------------| \mu
            ||           3                       4|    
            \\(d + H(xi))             (d + H(xi)) /    

                                 2\   /           a[2]    
   + lambda H(xi) + (v - 1) H(xi) / + |a[1] - ------------
                                      |                  2
                                      \       (d + H(xi)) 

                             2 a[4]   \                          
   + 2 a[3] (d + H(xi)) - ------------| (lambda + 2 (v - 1) H(xi)
                                     3|                          
                          (d + H(xi)) /                          

   \ /                                 2\            2 /    
  )| \mu + lambda H(xi) + (v - 1) H(xi) / - 2 alpha k  |a[1]
   |                                                   |    
   /                                                   \    

         a[2]                               2 a[4]   \ /  
   - ------------ + 2 a[3] (d + H(xi)) - ------------| \mu
                2                                   3|    
     (d + H(xi))                         (d + H(xi)) /    

                                 2\       /    
   + lambda H(xi) + (v - 1) H(xi) / - 2 w |a[0]
                                          |    
                                          \    

                          a[2]                      2
   + a[1] (d + H(xi)) + --------- + a[3] (d + H(xi)) 
                        d + H(xi)                    

         a[4]    \     /                            a[2]   
   + ------------| + k |a[0] + a[1] (d + H(xi)) + ---------
                2|     |                          d + H(xi)
     (d + H(xi)) /     \                                   

                     2       a[4]    \  
   + a[3] (d + H(xi))  + ------------|^2
                                    2|  
                         (d + H(xi)) /  
expr := simplify(%);
Error, (in simplify) too many levels of recursion
temp := algsubs(d+H(xi) = freeze(d+H(xi)), numer(expr));
                              expr
thaw(collect(temp, freeze(d+H(xi)))/denom(expr));
                              expr
collect(%, H(xi));