function JG=vrni_jg(X,l,beta);
  zacetni=[X(1);X(2);0;0;1;0;0;0;0;1;0;0];
  [t,u]=ode45(@sistem12,[0;l],zacetni,[],beta);
  JG=[u(end,7),u(end,11);u(end,8),u(end,12)];
