%RESISISTEM resi sistem enacb za veriznico z uporabo ukaza fsolve
%X=resisistem(X0,ZAC,L,M,barva)
%X je resitev (x in y koordinate)
%X0 zacetni vektor x-ov in y-ov
%X0 ena sama vrstica x-ov in y-ov zaporedoma vkljucno z robnima
%ZAC je matrika robnih tock
%ZAC=[x_zac x_kon;y_zac y_kon]
%L je vektor dolzin
%M je vektor mas
%BARVA je barva risanja
%
function X=resisistem(X0,ZAC,L,M,barva);
  %zapisemo zacetni vektor za xi eta in lambda
  n=length(L)-1;
  X00=ones(1,3*n+1);
  for i=1:n
    X00(i)=X0(i+1)-X0(i);% racunamo priblizke za xi
    X00(n+i)=X0(n+i+3)-X0(n+i+2);% racunamo priblizke za eta
  end
  Xres=fsolve('funsistem',X00,[],ZAC,L,M);
  Xres
  n=length(L)-1;
  X=zeros(2,n+2);%vektor x in y koordinat 
  X(:,[1,n+2])=ZAC;
  for i=2:n+1
    X(:,i)=ZAC(:,1)+sum([Xres(1:i-1);Xres(n+1:n+1+i-2)],2);
  end
  plot(X(1,:),X(2,:),barva);
