needsPackage "NumericalSchubertCalculus"
setRandomSeed 2

--Problem W11W21^2W31 in G(3,7)
 --a problem with 4 solutions

print("Solving problem W11W21^2W31 in G(3,7)");

SchPblm = randomSchubertProblemInstance(
  {{3, 1},{2, 1},{2, 1},{1, 1}},3,7);
time S = solveSchubertProblem(SchPblm, 3,7);
assert all(S,s->checkIncidenceSolution(s, SchPblm))

 end
 ------

restart
 load"NumericalSchubertCalculus/EXA/ProblemsG37/W11W21^2W31-G37.m2"

