Space‐charge, field, and potential distributions in inhomogeneously doped crystals are determined by a nonlinear differential equation which requires numerical computation in order to be solved. By means of numerical computation an attempt is made to obtain the potential, field, and space charge throughout the crystal. For exponential impurity distributions bounding curves can be obtained which allow one to obtain the complete distributions. A few examples are computed which demonstrate the difficulties encountered in obtaining solutions.