y[x_] := (a =
First[FindRoot[
ns - 1 == (-8 n^2 - 4 n)/(
3 ((-2 n)/3)^(
1/\[Alpha]) ((-3^(1 - \[Alpha])/2^(1 - \[Alpha]))^(1/(
2 \[Alpha] - 1)) *n - (2 \[Alpha] (n + 1) na)/(
2 \[Alpha] - 1))^((2 \[Alpha] - 1)/\[Alpha])) /. {ns ->
0.96, \[Alpha] -> x, na -> 55}, {n, -2}]]; 68 + n /. a)
FunctionRange[y[x], x, x]

那东西不是给数值函数用的。
自己写一个采样。

FunctionRange（至少目前）是处理不了黑箱函数的。一个更简单的例子（改编自文档）：
func = FunctionInterpolation[Sin[x^2 - 1] - 1/(x^2 + 2), {x, 0, 1}]
FunctionRange[{func[x], 0 <= x <= 1}, x, y]
(* FunctionRange::nmet: Unable to find the range with the available methods. *)