function oblika_vrvi(lambda0,mu0,parametri);
  %parametri=[L,H,ro,w,g];
  L=parametri(1);
  H=parametri(2);
  beta=parametri(5)/(parametri(4)^2*H);
  l=L/H;
  %resevanje z vgrajeno funkcijo fsolve
  X=fsolve(@vrni_g,[lambda0;mu0],optimset('fsolve'),l,beta);
  lambda=X(1)
  mu=X(2)

  %[lambda,mu]=newton(@vrni_gj,@vrni_g,lambda0,mu0,l,beta,1e-3)
  [t,u]=ode45(@sistem,[0;l],[lambda;mu;0;0],[],beta);
  plot(u(:,3)*H,u(:,4)*H,'b');
  axis equal;
  grid on;
