Constrained summations

Bert to write this up.

A[] = whatever()