global m n nt co T cL ome
m=11;n=5;nt=1;co=[5 1;1 2;2 3;2 3;2 3;3 4;3 4;3 4;4 5;4 5;4 5];
ome=1e4;
T=[5];cL=zeros(n,1);
om=2*pi*[1:0.05:30]*1e6;
sol=[];
x0=zeros(m+n+nt,1);
for jom=1:length(om)
 ome=om(jom);
 [x fval info]=fsolve(@felet3,x0);
 sol=[sol,x]; x0=x;
endfor
ii=sol(1:m,:);
v=sol(m+1:m+n,:);
iA=sol(m+n+1:m+n+nt,:);
result=[om' abs(sol)' angle(sol)'];
save -ascii result3.txt result
figure 1
plot(om/(2*pi),real(1./ii(1,:)),"linewidth",2,om/(2*pi),imag(1./ii(1,:)),"linewidth",2)
legend