clc;
//e.g 22.19
gm=2*10**-3;
rd=50*10**3;
Rs=1*10**3;
Ro=(Rs/gm)/(Rs+1/gm);
disp('ohm',Ro*1,"Ro=");