Solve[{x[n] == omega + a*x[n-1] + gamma*sqrt[x[n-1]]}, x]