This returns a list of μIJ(pd)/pd, or μfJ(pd)/pd, for d = 0,...,e.
As d approaches ∞, the sequence of these terms converges to the critical exponent of I, or of f, with respect to J.
i1 : R = ZZ/5[x,y]; |
i2 : I = ideal(x^2,x*y,y^2); o2 : Ideal of R |
i3 : m = ideal(x,y); o3 : Ideal of R |
i4 : criticalExponentApproximation(2,I,m) 4 24 o4 = {0, -, --} 5 25 o4 : List |
i5 : f = x^2 + y^3; |
i6 : criticalExponentApproximation(2,f,m) 3 19 o6 = {0, -, --} 5 25 o6 : List |