function thetabx = thetabx(x,m)
  mm = double(m);
  thetabx = 1/(x+mm);
