设函数y=1/√ax²+ax+1 的定义域为R,求实数a的取值范围