In classical/Newtonian mechanics, the doppler shift (for light) can be expressed as:
$$ \frac{f_r-f_s}{f_s}=\frac{1+\beta}{1-\beta}-1 $$
In relativity the doppler shift can be expressed as:
$$ \frac{f_r-f_s}{f_s}=\sqrt{\frac{1+\beta}{1-\beta}}-1 $$
In the Newtonian limit where $ \beta\ $is small, we should be able to approximate both of the above expressions in terms of first order Taylor expansions. If we do so, we get the following approximations (obtained from WolframAlpha):
$$ \frac{f_r-f_s}{f_s}=2\beta\ (Newtonian) $$ $$ \frac{f_r-f_s}{f_s}=\beta\ (Relativitic) $$
However, in the Newtonian limit these approximations should be the same (i.e. we expect the relativistic doppler shift equation to give the same results as the Newtonian doppler shift equation when $ \beta\ $is small). Why do these last two equations not match? Where is the flaw in my logic?
Definitions:
$ f_s $: Source frequency
$ f_r $: Received frequency
$ \beta=v/c $