Please run the Matlab code below and generate the IV curve for VG = 100 mV. Other parameters will remain the same. Please only upload the figures. Other parameters will remain the same. set(h, linewidth,[2.01) set(gca, Fontsize'.[25] xlabel('D(E)--->') ylabel(' E (V) --->') grid on % Fig. 1.4.6 clear all %Constants (all MKS, except energy which is in eV) hbar=1.055e-34;q=1.602e-19;10=q*q/hbar; %Parameters UO=0.025;kT=0.025;mu=0;ep=0.2: gl=0.005;g2=0.005;g=gl+g2; alphag=1;alphad=0.5; %Energy grid NE=501;E=linspace(-1,1,NE);dE=E(2)-E(1); D=(g/(2*pi))./((E. 2)+((g/2) 2));% Lorentzian Density of states per eV D=D./(dE*sum(D));%Normalizing to one %Bias IV=101:VV=linspace(0,1.IV); for iV=1:IV Vg=0;Vd=VV(IV); %Vd=0;Vg=VV(IV); mul=mu;mu2=mul-Vd:UL=-alphag*Vg)-(alphad*Vd); U=0;%Self-consistent field dU=1; while dU>le-6 f1=1./(1+exp((E+ep+UL+U-mul)./KT)); f2=1./(1+exp((E+ep+UL+U-mu2)./kT)); N(V)=dE*sum(D.*((f1.*gl/g)+(f2.*g2/g))); Unew=U0*N(V);dU=abs(U-Unew); U=U+0.1 *(Unew-U); end I(IV)=dE*10*(sum(D.*(f1-f2)))*(gl*g2/g); end hold on h=plot(VV,N,'b'); %h=plot(VV.I,'b'); set(h,'linewidth', [2.01) set(gca, Fontsize' [25]) xlabel(Voltage (V)--->') %ylabel(' Current (A)--->') ylabel(' Number of electrons --->') grid on