function yprime = odefun(x,y)
% function yprime = odefun(x,y)
% Kartik's MATLAB code
% November 19, 2007
% User file that contains the ODE function
% after writing the 2nd order ODE as a system of 2 first order ODEs

yprime = [y(2); 1/8*(32+2*x^3-y(1)*y(2))];