[xv yv nv nc con R]=circuito(50,50,1,1);
a=matVtoI(nv,nc,con,R);
atil=contorno(a,50,50);
[s50 l50]=eig(atil);
d=diag(l50);
[sd id]=sort(d);
sd(541)
sd(542)
id(542)
y=s50(:,46);
y2=reshape(y,50,50);
contourf(y2)
pause
surf(y2)
pause
y=s50(:,id(543));
y2=reshape(y,50,50);
contourf(y2)
pause
surf(y2)
pause
y=s50(:,id(544));
y2=reshape(y,50,50);
contourf(y2)
pause
surf(y2)
pause
y=s50(:,id(545));
y2=reshape(y,50,50);
contourf(y2)
pause
surf(y2)
pause
y=s50(:,id(548));
y2=reshape(y,50,50);
contourf(y2)
pause
surf(y2)
pause
y=s50(:,id(558));
y2=reshape(y,50,50);
contourf(y2)
pause
surf(y2)
pause
y=s50(:,id(588));
y2=reshape(y,50,50);
contourf(y2)
pause
surf(y2)
pause