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