%diodepara;
I=[11.324,9.332,7.341,5.353,3.368];R=1000;I=I/R;
V=[676.432,669.241,658.851,647.498,632.338];V=V/R;
y=log(I);VT=25.9*0.001;
plot(V,y);title('ln(ID) Current vs V');xlabel('VD');ylabel('ln(ID)');
p=polyfit(V,y,1);m=p(1,1);b=p(1,2);
IS=exp(b)
N=1/(m*VT)


    
