function G=vrni_g(X,l,beta);
  [t,u]=ode45(@sistem,[0;l],[X(1);X(2);0;0],[],beta);
  G=[u(end,3);u(end,4)-1];
