function y=xSqrt(t,x0)

if (x0 < 0) y=-sqrt(2*t+x0.^2);
else y= sqrt(2*t+x0.^2); end;
