> (2/3)^(-3/2);
                                      1/2  1/2
                                 3/4 3    2
    > a := I+2^(1/2)+3^(1/2):
    > a := expand(a^2);
                                 1/2        1/2      1/2  1/2
                   a := 4 + 2 I 2    + 2 I 3    + 2 2    3
    > combine(a,power);
                                 1/2        1/2      1/2
                        4 + 2 I 2    + 2 I 3    + 2 6
    > rsolve( {y(n) = n*y(n-1) + (n-1)!, y(0) = 1}, y(n) );
                   GAMMA(n + 1) (Psi(n + 1) + gamma + 1)
    >  isolve(2*x^2-7*y^2+5*z^2);
               _N3 %1                _N3 %3                _N3 %2
     {z = ----------------, x = ----------------, y = ----------------}
          igcd(%1, %2, %3)      igcd(%1, %2, %3)      igcd(%1, %2, %3)
                               2                    2
    %1 :=               - 5 _N1  - 4 _N1 _N2 + 2 _N2
                                   2        2
    %2 :=                     5 _N1  + 2 _N2
                             2                     2
    %3 :=               5 _N1  - 10 _N1 _N2 - 2 _N2