In this paper, a novel numerical technique, the first-order Smooth Composite Chebyshev Finite Difference method, is presented. Imposing a first-order smoothness of the approximation polynomial at the ends of each subinterval is originality of the method. Both round-off and truncation error analyses of the method are performed beside the convergence analysis. Diffusion of oxygen in a spherical cell including nonlinear uptake kinetics is solved by using the method. The obtained results are compared with the existing methods in the literature and it is observed that the proposed method gives more reliable results.