%FUNSISTEM vrne vrednost funkcije, katere niclo iscemo
%Z=FUNSISTEM(X,L,M);
%Z je vrednost funkcije
%X so verdnosti spremenljivk
%X(1:n+1) ksiji, X(n+2:2*n+2) ete, X(2*n+3:end) lambde
%L je vektor dolzin
%M je vektor mas
%
function Z=funsistem(X,ZAC,L,M);
  n=length(L)-1;
  mi=zeros(1,n);%tabela mi-jev
  for i=1:n
    mi(i)=1/2*(M(i)+M(i+1));
  end
  xi_koncni=ZAC(1,2)-ZAC(1,1)-sum(X(1:n));
  eta_koncni=ZAC(2,2)-ZAC(2,1)-sum(X(n+1:2*n));
  Z=zeros(1,3*n+1);
  for i=1:n-1
    Z(i)=X(2*n+i)*X(i)-X(2*n+i+1)*X(i+1);
  end
    Z(n)=X(3*n)*X(n)-X(3*n+1)*xi_koncni;
  for i=n+1:2*n-1
    Z(i)=X(i+n)*X(i)-X(i+n+1)*X(i+1)+1/2*mi(i-n);
  end
    Z(2*n)=X(3*n)*X(2*n)-X(3*n+1)*eta_koncni+1/2*mi(n); 
  for i=2*n+1:3*n
    Z(i)=X(i-2*n)^2+X(i-n)^2-L(i-2*n)^2;
  end
    Z(3*n+1)=xi_koncni^2+eta_koncni^2-L(n+1)^2;
