function yp=f(x,y)

yp=x+y;     % yp=y'
            % y'=x-y
%yp=2*x;

%yp=0.0225*y-0.0003*y^2;

% ODE is y'=f(x,y)
% in this example,

% y'=x-y