function g = deltav3(ii)
global m n ome
v0=1.0; RR= 1000*[0 0 4.97 0 0 3.338 0 0 2.702 0 0]; big=1e12;
LL=1e-6*[0 4.49 0 12.5 0 0 2.26 0 0 0.522 0]; CC=1e-12*[0 43.9 big big 22.9 big big 30.3 big big 57.1];
for k=1:m
 if (k!=1) 
  Z(k)=RR(k)+i*ome*LL(k)-i/(ome*CC(k));
  g(k)=-Z(k)*ii(k);
 else
  g(k)=v0;
 endif
endfor
endfunction