function y = myfunc3ODE(t,x, f0,dfun)
%
% Homotopy function used in ODE solver.
%
dy = feval(dfun,x);
y  = - dy \f0;
