An optimization method, based on linear programming, has been developed that determines the enrichment distribution in a multizone reactor such that the desired power distribution is achieved at beginning-of-life conditions, as well as the poison distribution during reactor operation such that the reactor is kept critical and the desired power distribution is retained. Results of the numerical application of the method are presented.