function g = deltav2(ii) global m n ome v0=1.0; RR= 10*[0 10 0 100 0 0 0 10]; LL=1e-3*[0 0 10 10 20 0 40 0]; CC=1e-6*[0 1000000 1 1000000 1000000 1 10 100000000]; 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