wolffd@0: function pot = mpot(members, sizes, logp, mu, Sigma)
wolffd@0: % MPOT Make a moment Gaussian potential.
wolffd@0: % pot = mpot(members, sizes, logp, mu, Sigma)
wolffd@0: %
wolffd@0: % All params default to 0 if omitted.
wolffd@0: 
wolffd@0: n = sum(sizes);
wolffd@0: if nargin < 3, logp = 0; end
wolffd@0: if nargin < 4, mu = zeros(n,1); end
wolffd@0: if nargin < 5, Sigma = zeros(n,n); end
wolffd@0:   
wolffd@0: pot.domain = members;
wolffd@0: pot.sizes = sizes;
wolffd@0: pot.logp = logp;
wolffd@0: pot.mu = mu;
wolffd@0: pot.Sigma = Sigma;zeros(n,n);      
wolffd@0: pot = class(pot, 'mpot');