An analytical method for continuous refuelings is developed. The starting point is the first-order nonlinear differential equation, and the technique for solving this is a use of an elliptic function theory. The nature of the neutron flux distribution, eigenvalue, and critical velocity are investigated for the in-out and the out-in continuous refueling. Numerical results for neutron flux distribution are given.