function dy = riccati(x,y)
% dy = riccati(x,y)
% Example of derivative function
% for solving differential equation
dy = y.^2 - x;