I'm looking to find out the depletion width in a MOSFET including incomplete ionization. A closed form equation is not necessary. Any numerical method to figure it out using integration would be good enough and very helpful. I just want to compare it to the complete ionization approximation.