A nonlinear programming method is used for finding an equitable optimal fair division of the unit interval [0, 1) among n players. Players' preferences are described by nonatomic probability measures μ 1 , . . . , μ n with density functions having piecewise strict monotone likelihood ratio property. The presented algorithm can be used to obtain also an equitable ε-optimal fair division in case of measures with arbitrary differentiable density functions. An example of an equitable optimal fair division for three players is given.