%RISIVER narise diskretno veriznico
%RISIVER(T1,T2,l,ro,tol,nt,barva)
%T1 in T2 sta krajisci veriznice oblike
%T1=[a,A], T2=[b,B], l je dolzina, ro
%dolzinska gostota, tol toleranca pri iteraciji,
%nt stevilo tock risanja in barva barva izrisa

function risiver(T1,T2,l,ro,tol,nt,barva);
  z=isciz(T1,T2,l,1,tol); %z0=1
  a=T1(1);A=T1(2);b=T2(1);B=T2(2);
  k1=(B-A)/l;
  k2=atanh(k1);
  u=k2-z;
  v=k2+z;
  g=9.8; %tezni pospesek
  C0=(b-a)/(2*z);
  D=(b*u-a*v)/(2*z);
  x=linspace(a,b,nt);
  lambda0=(C0*cosh((a+D)/C0)-A);
  y=C0*cosh((x+D)/C0)-lambda0;
  plot(x,y,barva);
  grid on
  
