-- -*- M2-comint -*- hash: 2157229915827826967

i1 : variables = declareVariable \ {x,y}

o1 = {x, y}

o1 : List

i2 : params = declareVariable \ {a,b}

o2 = {a, b}

o2 : List

i3 : F = gateSystem(matrix{params}, matrix{variables}, matrix{{a*x*y-1},{x^3+y^2-b}})

o3 = gate system: 2 input(s) ---> 2 output(s) (with 2 parameters)

o3 : GateSystem

i4 : PH = parametricSegmentHomotopy F;

i5 : parameters PH

o5 = {{A , A , B , B }}
        a   b   a   b

o5 : GateMatrix

i6 : (a0,b0) = (1,2); startSolution = point{{1,1}};

i8 : (a1,b1) = (2,1);

i9 : H01 = specialize(PH, matrix{{a0,b0,a1,b1}});

i10 : targetSolution = first trackHomotopy(H01,{startSolution})

o10 = targetSolution

o10 : Point

i11 : assert(norm evaluate(F,matrix{{a1,b1}},matrix targetSolution) < 0.0001)

i12 : 
