-- -*- M2-comint -*- hash: 13249572421015205798

i1 : G = graph({0,1,2,3}, {{0,1},{1,2},{2,3},{0,3}});

i2 : R = oscRing(G, Reduced => true);

i3 : I = oscSystem(G, R);

o3 : Ideal of R

i4 : Jac = oscJacobian(G, R);

             4      4
o4 : Matrix R  <-- R

i5 : realsols = findRealSolutions I;
warning: some solutions are not regular: {9, 10, 14, 18, 19}

i6 : select(realsols, S -> isStableSolution(Jac, S))

o6 = {{1, 1, 1, 0, 0, 0}}

o6 : List

i7 : 
