A method is for establishing a nuclear reactor core loading pattern (LP) for fuel assemblies and burnable absorbers (BAs). The method establishes an optimum LP through the steps of: a) providing nuclear data representing fuel assemblies and BAs in a nuclear reactor core; b) depleting the nuclear data to form a reference core depletion; c) incorporating the nuclear data into a system of linear equations of a nuclear design quality flux solution method; d) defining the system of linear equations to include constraints which accurately represent the neutron physics of the reactor; employing the equations as a constraint matrix for a MIP solver to find an optimum core pattern solution; f) repeating steps b) through e) updating the constraints and objective functions to satisfy specified engineering requirements and establish an optimum core loading pattern. An algorithm for deriving the system of equations is also disclosed.