function rep_pot=replicate_pot(table,over,dimsizes)
%absorption_pot replicates pot table over the dimensions in over
if ~isempty(over)
rep_pot=squeeze(repmat(table,[ones(1,length(size(table))) dimsizes(over)]));

index=1:length(dimsizes);
ind=mysetdiff(index,over);
order=[ind over];
id=eye(length(dimsizes));
permatrix=id(order,:);
perm=[1:length(dimsizes)]*permatrix;
rep_pot=permute(rep_pot,perm); %permute dimensions
%put back relevant dimension of size one
rep_pot=reshape(rep_pot,dimsizes);
else

    rep_pot=table; 
    rep_pot=reshape(rep_pot,dimsizes);%put back relevant dimension of size one
end