% Creates the function odef(x,y)
% Input: x (a real number), y (a real number)
% Output: dy 

function dy = odef(x,y)
dy = 0.25*y*log(1000/y);