% OPTFUN vrne vrednost optimizacijske funkcije veriznice
% X je vektor neznank x_i=X(i+1), y_i=Y(n+3+i), i=0,1,n+2
% L je vektor dolzin palic
% M je vektor mas palic
%Z je vrednost funkcije
function [Z]=optfun(X,L,M);
  n=length(L)-1; % stevilo palic je n+1
  Z=0;
  for i=1:n+1;
    Z=Z+M(i)/2*(X(n+3+i-1)+X(n+3+i));
  end
