A general fractional Taylor formula and its computation for insufficiently smooth functions are discussed. The Aitken delta square method and epsilon algorithm are implemented to compute the critical orders of the local fractional derivatives, from which more critical orders are recovered by analysing the regular pattern of the fractional Taylor formula. The Richardson extrapolation method is used to calculate the local fractional derivatives with critical orders. Numerical examples are provided to verify the theoretical analysis and the effectiveness of our approach.