SNES with Finite Difference Jacobian Approximations

Dear all,
I was wondering if we could use SNES to solve nonlinear PDEs with Finite Difference Jacobian Approximations. Are there any examples? It is not easy to compute Jacobian matrices analytically in my simulation.
Thanks