function Z = bowl(X,A)

Z = diag(X*A*X');
