i1 : S = ZZ/32009[x_0..x_3];  | 
i2 : I = monomialCurveIdeal(S,{1,2,3});
o2 : Ideal of S
 | 
i3 : X = ideal(gens I * random(source gens I,
             S^{1:-min(degrees source gens I)}));
o3 : Ideal of S
 | 
i4 : time h1NormalBundle(I,X)
     -- used 0.00675235 seconds
o4 = 0
       ZZ
o4 : ------module
     32009
 | 
This function might require quite an amount of time to run, depending on the input data.