We analyze a (possibly degenerate) second order mean field games system of
partial differential equations. The distinguishing features of the model
considered are (1) that it is not uniformly parabolic, including the first
order case as a possibility, and (2) the coupling is a local operator on the
density. As a result we look for weak, not smooth, solutions. Our main result
is the existence and uniqueness of suitably defined weak solutions, which are
characterized as minimizers of two optimal control problems. We also show that
such solutions are stable with respect to the data, so that in particular the
degenerate case can be approximated by a uniformly parabolic (viscous)
perturbation