Applications of Fréchet derivative emerge in the sensitivity analysis of matrix functions. Our work extends the generalized complex step approximation using the complex computation f (A + e iθ hE) as a tool to matrix case, and combines it with finite difference formula to estimate the Fréchet derivative. We provide numerical results for the approximation to the first and the second order Fréchet derivative of the matrix exponential and matrix square root.