function thetabxlam = thetabxlam(x,m,lam)
thetabxlam = (1.0/(3.0+2.0*lam))*x*(1 + lam*x*x*sqrt(x)) ...
  + (m-1)*(2.0+lam)/(3.0+2.0*lam);
