maxVds=15;
vds=linspace(0,maxVds,1e3);
vto=2.5;
kp=0.2;
W=1;
L=1;
vgsVector=linspace(2,10,4);

figure; h=subplot(1,1,1); set(h,'FontSize',12,'FontWeight','bold');  hold on; grid on;
xlabel('V_{DS} (Volts)'); ylabel('I_{DS} (mA)');
axis([0,maxVds,-1,6]);
for vgs=vgsVector
    Ids=nmos(vds,vgs,kp,W,L,vto);
    plot(vds,Ids,'b');
    str=sprintf('V_{GS}=%s',num2str(vgs));
    if vgs==2
        h=text(mean(vds)-.5,Ids(end)-.3,str);
    else
        h=text(mean(vds)-.5,Ids(end)+.3,str);
    end
    set(h,'FontWeight','bold');
end


