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