ID = 0.5K(VGS - VT)2(1+λVDS) ... (1)
Now,
VDS = VD - VS
= (VDD - IDRd) - (IDRs)
= VDD - ID(Rd + Rs)
= 10 - 31k×ID ... (2)
VGS = VG - VS
= VDD×R2/(R1 + R2) - IDRs
= 1.25 - 1k×ID ... (3)
Putting (2) and (3) in (1) you shall get,
ID = 250µ(0.25 - 1k×ID)(1.2 - 620×ID)
ID = 250µ(0.3 - 1355×ID + 620k×ID2)
The above eqn gives 2 values of ID = 8.58 mA and 56.4 µA
For ID = 56.4 µA and VDS = (10 - 31k×56.4µ) V = 8.25 V
For ID = 8.58 mA, VDS becomes negative
Hence, ID = 56.4 µA and VDS = 8.25 V