function G=dUdx(x,U)
if x==0
    G=[U(2); gavu1(U(1))/3];
else
    G=[U(2); gavu1(U(1))-2/x*U(2)];
end