Let f be a function k times continuously differentiable in [a, b], then we will prove that the RAFU method provides a sequence (H-n)(n) defined in [a, b] such that for each j - 0,..., k parallel to f((j)) - H-n((j))parallel to <= [M - m/root n + omega (f((k)), b - a/n)] (b - a)(k-j) being n >= 2, parallel to.parallel to the uniform norm, M and m the maximum and the minimum of f((k)) in [a, b] respectively and omega(f((k)), b-a/n) its modulus of continuity. The called RAFU remainder in Taylor's formula will be presented. The simultaneous approximation problem will be solved from average samples, from linear combinations, from local average samples given by convolution, from approximate values and in the case of non-uniformly spaced data of f((k)). We will also study the numerical differentiation case. Our approach is easily realizable in computations. Some examples will be given.