Question: why not to collect R(xi) accurately

restart;

T := R(xi)*R(xi) + lambda;

u := A[0] + A[1]*R(xi) + B[1]/R(xi);

d[1] := A[1]*T - B[1]*T/R(xi)^2;

d[2] := 2*A[1]*R(xi)*T - 2*B[1]*T/R(xi) + 2*B[1]*(R(xi)^2 + lambda)*T/R(xi)^3;

expand(((-alpha^2*b^2 + a^2)*alpha^2)/(2*beta)*d[2] + (omega + alpha^2*(alpha^2*l^2 + k^2)/2 - a*C[1]/(-alpha^2*b^2 + a^2))*u[0]/(beta - 2*beta*a^2/(-alpha^2*b^2 + a^2)) + u[0]*u[0]*u[0]);

value(%);

simplify(%);

collect(%, R(xi));

      /      6  4    4      2\      3
 A[1] \-alpha  b  + a  alpha / R(xi) 
 ------------------------------------
             /     2  2    2\        
        beta \alpha  b  + a /        

                  /      6  4    4      2\         
      A[1] lambda \-alpha  b  + a  alpha / R(xi)   
    + ------------------------------------------ + 
                     /     2  2    2\              
                beta \alpha  b  + a /              

                         /                               /     
             1           |/                    B[1] \    |     
   --------------------- ||A[0] + A[1] R(xi) + -----|[0] |beta 
        /     2  2    2\ \\                    R(xi)/    \     
   beta \alpha  b  + a /                                       

                                                  2
   /     2  2    2\ /                    B[1] \    
   \alpha  b  + a / |A[0] + A[1] R(xi) + -----|[0] 
                    \                    R(xi)/    

      1  2  2      6   1 /  2  2    2  2\      4
    + - b  l  alpha  + - \-a  l  + b  k / alpha 
      2                2                        

                                                       \\
      /  1  2  2    2      \      2    2               ||
    + |- - a  k  + b  omega| alpha  - a  omega + a C[1]||
      \  2                 /                           //

                  /      6  4    4      2\
      B[1] lambda \-alpha  b  + a  alpha /
    + ------------------------------------
               /     2  2    2\           
          beta \alpha  b  + a / R(xi)     

            6  4       2         4      2       2     
      -alpha  b  lambda  B[1] + a  alpha  lambda  B[1]
    + ------------------------------------------------
                     /     2  2    2\      3          
                beta \alpha  b  + a / R(xi)