Daniel@0: function pot = upot(domain, sizes, p, u)
Daniel@0: % UPOT Make a discrete utility potential.
Daniel@0: % pot = upot(domain, sizes, p, u)
Daniel@0: %
Daniel@0: % sizes(i) is the size of the i'th domain element.
Daniel@0: % p defaults to all 1s, u defaults to all 0s.
Daniel@0: 
Daniel@0: if nargin < 3, p = myones(sizes); end
Daniel@0: if nargin < 4, u = 0*myones(sizes); end
Daniel@0: 
Daniel@0: pot.domain = domain;
Daniel@0: pot.p = myreshape(p, sizes);
Daniel@0: pot.u = myreshape(u, sizes);
Daniel@0: pot.sizes = sizes(:)';
Daniel@0: pot = class(pot, 'upot');