function g = deltav(ii)
global m n ome
v0=1.0; RR= 100*[0 1 1 1 1 1 1 1];
LL=0e-4*[0 1 1 1 1 1 1 1]; CC=1e+16*[0 1 1 1 1 1 1 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