%FU je funkcija katere nièlo išèemo
%u je spremenljivka funkcije
%L je vektor dolžin palic dim. n+1 = 2*p+1
%x0 je zaèetek xn1 pa konec verižnice
%cmi je kumulativna vsota mi-jev z niclo na
%prvem mestu
function z=fu(u,L,cmi,x0,xn1);
n=length(L)-1;
p = n/2;
m=cmi(p+1)-cmi;
v = L./sqrt(1+u^2*m.^2);
z = sum(v)-(xn1-x0);