Parse Lorentz cone constraints from formula

[RussTedrake]

Now we can write, e.g.

prog.AddConstraint(y[0] >= np.linalg.norm(x))

and get a Lorentz cone constraint back out.

+@hongkai-dai for feature review, please.

---

This change is