OK -- if I convert to float it works, sort of. I have to comment out (as you did) one of the constraints. If I leave them all in, it crashes python. And, of course, w/ all the constraints (uncomment the last one), it should work (or at least it does on paper, and in maple). thanks! sd