Finally got this shit to give feasible trajectories.  Problem was using >= instead of <= convention for solver.  I think its correct now.  Currently comparing sequence counts with theoretical bounds.