% right-hand side of simple IVP with exact solution expressible
% in terms Bessel functions

function yDot = f(t,y)

yDot = y^2 + t^2;