Daniel@0: function pot = cpot(members, sizes, g, h, K)
Daniel@0: % CPOT Make a canonical Gaussian potential.
Daniel@0: % pot = cpot(members, sizes, g, h, K)
Daniel@0: %
Daniel@0: % All params default to 0 if omitted.
Daniel@0: 
Daniel@0: n = sum(sizes);
Daniel@0: if nargin < 3, g = 0; end
Daniel@0: if nargin < 4, h = zeros(n,1); end
Daniel@0: if nargin < 5, K = zeros(n,n); end
Daniel@0:   
Daniel@0: pot.domain = members;
Daniel@0: pot.sizes = sizes(:)';
Daniel@0: pot.g = g;
Daniel@0: pot.h = h;
Daniel@0: pot.K = K;
Daniel@0: pot = class(pot, 'cpot');