若函数f(x)=lg(x^2 - 2ax + 3a)的定义域为R,求实数a的取值范围