clear; clc; close;

yos = 20*10^(-6);
Idss = 10*10^(-3);
Idq = 10*10^(-3);
Vp = -4;
Vdsq = Vp/4;
Rg = 10*10^(6);
Av = -8;
Vdd = 20;

Vgsq = Vp/4;
Id = Idss*(1-(Vgsq/Vp))^2;

Rs = -Vgsq/Id;

gmo = 2*Idss/abs(Vp);
gm = gmo*(1-(Vgsq/Vp));
Rd = -Av*(1+gm*Rs)/gm;

disp(Rs,'Rs(ohms) = ');
disp(Rd,'Rd(ohms) = ');