function yp = singer( x, y )
%
% Scalar 1st order ODE with a movable pole
% y' = x^2 + y^2
%
yp = x^2 + y.^2;
