Function: algisinv
Section: algebras
C-Name: algisinv
Prototype: iDGGD&
Help: algisinv({al},x,{&ix}): tests whether x is invertible in al (Hamilton
 quaternions if omitted) and sets ix to the inverse of x.
Doc: Given an element $x$ in \var{al} (Hamilton quaternions if omitted), tests
 whether $x$ is invertible, and sets $ix$ to the inverse of $x$.
 \bprog
 ? A = alginit(nfinit(y), [-1,1]);
 ? algisinv(A,[-1,1]~)
 %2 = 0
 ? algisinv(A,[1,2]~,&ix)
 %3 = 1
 ? ix
 %4 = [Mod(Mod(-1/3, y), x^2 + 1), Mod(Mod(2/3, y), x^2 + 1)]~
 @eprog

 Also accepts square matrices with coefficients in \var{al}.
