function xdot = ode1(t,x);
% ODE1: computes xdot.
xdot = x.^2-t.^2./(t.^4+1);