function m = moment(b)

% Input: b kx1 parameter vector
% Output: m nxl matrix of "moment conditions".
% Rows (columns) correspond to observations (moment conditions).

global x

m = [x-b,(x-b).^2-1];
