% compute id = f(vds) for various values of vgs

clear idnest;

W=0.375;
L=0.25;
vsb = 0;
vds = [0:0.1:2.5];
vgs = [0:0.5:2.5];
idnest = zeros(length(vds), length(vgs));
for i = 1:length(vgs)
   idnest(:,i)=nmos025(W,L,vgs(i),vds,vsb)';
end
plot(vds,idnest,'k-');
xlabel('V_{DS} (V)');
ylabel('I_{D} (A)');