Abstract: A program for solving multi-group neutron diffusion eigenvalue and transient problems using JFNK method was developed based on MOOSE platform. Albedo boundary condition and vacuum boundary condition were realized in this program. The program was verified by TWIGL benchmark, and the simulation results are in good agreement with the reference. In this program, full implicit Newton method was used to solve eigenvalue problem, and the results show that Newton method can greatly reduce the nonlinear steps of iteration and speed up convergence compared with classical PI (power iteration) method. Spatial-temporal dynamics was used to solve the transient problem, and any point in the domain with the change of time could be simulated. Backward Euler difference was adopted in time discretization, and the simulation results at time step 0.01 s and 0.05 s are matched very well with the reference at time step 0.001 s, so the convergence and accuracy can be ensured with a bigger time step in the program.