Abstract: Nuclear data are fundamental data for nuclear energy, nuclear engineering, nuclear technology applications and nuclear science research, and have always received national attention. In particular, the importance of neutroninduced fission nuclear reaction data is unquestionable and the theoretical evaluation of nuclear data is required for a long period of time. The reaction mechanism of neutron induced fission reaction is very complicated and needs to be described by several different nuclear reaction models, such as optical potential model and compound nuclei model and so on. Nowadays, the most famous international nuclear reaction program systems are TALYS and EMPIRE. In China, the program FUNF is a fission nuclear reaction program system with our own independent intellectual property rights, which was developed by Prof. Zhang Jingshang of the China Nuclear Data Center. Each of these program systems need more than 40 parameters to be adjusted. Therefore, it is necessary to optimize a large number of parameters at the same time in order to correctly describe the nuclear fission reaction data. In this work, the FUNF program system was used to calculate the fission reaction data of neutron induced 238U. The fast neutron reaction data of fissile materials with incident energy from about 1 keV up to 20 MeV were calculated by FUNF code. It consists of the spherical optical model, the HauserFeshbach model, and the angular momentum dependent exciton model to describe the emission from compound nucleus to the discrete levels of the residual nuclei in preequilibrium process, meanwhile equilibrium reaction process was described by the HauserFeshbach model with width fluctuation correction. Based on the FUNF nuclear fission reaction model, the parameters of the FUNF program system was adjusted by the advanced optimization method MINUIT. MINUIT is a numerical minimization software library in the Fortron by CERN staff Fred James in the 1970 s. Then, combined with the parallel calculation method of MPI (message passing interface), preliminary calculations for the neutron induced fission reaction of 238U were performed. The results show that the FUNF parameters obtained by the MINUIT optimization method can well describe the fission reaction data of 238U. Next, the procedure will continue to be optimized to further improve the efficiency of the tuning parameters and the accuracy of the parameters. And the application to the relevant fission nuclear reactions will be developed to improve the accuracy of our fission nuclear reaction data.