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