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);
ii=0*ones(m,1);
v=0*ones(n,1);
iA=0*ones(nt,1);
x0=[ii;v;iA];
res=felet(x0)
[x fval info]=fsolve(@felet,x0)