function [f] = f(x,y)
% Evaluation of right hand side of a differential
% equation.  

f = 1/y;


