global m n nt co T cL ome
m=8;n=7;nt=2;co=[1 2;2 4;4 3;4 6;4 5;3 5;6 5;6 7];
ome=1e4;
T=[1 7];cL=zeros(n,1);
k=1:320; ## 160;
om=100*((2.0)^(1/32)).^k;
sol=[];
x0=zeros(m+n+nt,1);
for jom=1:length(om)
 ome=om(jom);
 ##res=felet(x0)
 [x fval info]=fsolve(@felet2,x0);
 ##info
 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 result.txt result
figure 1
loglog(om/(2*pi),abs(ii),"linewidth",2)
legend
figure 2
loglog(om/(2*pi),abs(v(3:6,:)),"linewidth",2)
legend
figure 3
semilogx(om/(2*pi),angle(ii),"linewidth",2)
legend
figure 4
semilogx(om/(2*pi),angle(v(3:6,:)),"linewidth",2)
legend