par=[0.3,2.1,2];%parametri sistema so [visinska razlika,dolzina vrvi,razpon]
%/ paramaters of the sistem are [hight difference, length of wire,
%distance]
[X,k]=newton('sistem','odvod',[1;1],1e-6,50,par);%resimo sistem enacb / solving system
a=X(1);%resitvi / the solution
b=X(2);
t=linspace(0,par(3));
x=linspace(b-par(3),b);%obmocje, kjer risemo verizico 
plot(t,a*cosh(x/a)-a*cosh((b-par(3))/a));%narisemo veriznico / plotting the catenary