clc;
Vdd=20;
Id=0.01;
Rd=1000;
Vgs=Vdd-Id*Rd;
disp('V',Vgs,"Vgs=")//The answers vary due to round off error