I once used the following MAGMA-code (just for finite groups, I don't know if it works for general reductive groups).
intrinsic Presentation(R::RngInvar) -> SeqEnum
{A presentation of the invariant ring R.}
fund := FundamentalInvariants(R);
prim := PrimaryInvariants(R);
sec := IrreducibleSecondaryInvariants(R);
invar := prim cat sec;
P := PolynomialRing(BaseRing(R), #fund);
A := Algebra(R);
invarpres := [];
for f in invar do
b,g := HomogeneousModuleTest(fund,[R!1],f);
Append(~invarpres, g[1]);
end for;
rel := RelationIdeal(R);
phi:=hom<A->P|invarpres>;
return ideal<P|[phi(r) : r in Basis(rel)]>;
end intrinsic;
Put this into a file, say "pres.m" and then you can do:
> Attach("pres.m");
> G:=MatrixGroup<2,Rationals() | [-1,0,0,-1] >;
> G;
MatrixGroup(2, Rational Field)
Generators:
[-1 0]
[ 0 -1]
> R:=InvariantRing(G);
> Presentation(R);
Ideal of Polynomial ring of rank 3 over Rational Field
Order: Lexicographical
Variables: $.1, $.2, $.3
Homogeneous
Basis:
[
$.1*$.3 - $.2^2
]
This is an example of the $A_2$-singularity