Matrix Inversion ANA

[churichuro]

This program will invert Matrixes

EXPORT ANA(a)
BEGIN
local i,j,k,r,p,b,n;
r:=1;
n:=rowDim(a);
b:=MAKEMAT(0,n,n);
FOR k FROM 1 TO n DO
p:=a[k,k];
FOR j FROM 1 TO n DO
b[k,j]:= -a[k,j];
b[j,k]:= a[j,k];
END;
b[k,k]:=r;
FOR i FROM 1 TO n DO
IF i=k THEN CONTINUE; END;
FOR j FROM 1 TO n DO
IF j=k THEN CONTINUE; END;
b[i,j]:=(p*a[i,j]-a[i,k]*a[k,j])/r;
END;
END;
a:=b;
r:=p;
END;
return(a/p);
END;