Suppose you have samples of uniform interval of a sinusoidal signal that is some sum of sinusoids. Forget about aliasing, as the signal is bandlimited. If we have infinite number of samples at sampling rate right at Nyquist rate, we can technically infer the whole bandlimited function and of course differentiate. We cannot do that for finite number of samples.
Thus, my question is how far can numerical differentiation be accurate for each frequency - that is deviation from ideal frequency response - given $n$ samples of uniform interval. One can assume that the filter being used is FIR filter.
Parks-McClellan is often used to obtain an optimal filter, so my question can be rephrased as how required maximum frequency response error magnitude (relative to ideal frequency response of a differentiator) corresponds to upper bound on number of samples/taps/order required for the derived filter.