> (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