In this paper, 3-D neutron diffusion equation was discretized into the generalized eigenvalue problem of large sparse matrix with finite difference method. Then, keff and power distributation were obtained by nonlinear algorithm JFNK(Jacobian-Free Newton-Krylov). Numerical results demonstrated that JFNK algorithm could get high accuracy for 3-D neutron diffusion equation.