I solve problem which appeared for a few profiles. The Newton method was not able to find D inside Dmin, Dmax because numerical step was too large.