function [hh,H] = toystrict(m,n)
hh = [];
for i=n-1:-1:1
  hh(end+1) = (m^i-1)/m^i;
end
H = hh2H(hh);
return
H = 1;
for i=1:n-1
  x = 1;
  for j=1:i
    x = x * hh(j);
  end
  H = H+x;
end


return