function dy=rhs(t,y)
dy = zeros(2,1);    % a column vector
dy(1) = fun_f(y(2),y(1));
dy(2) = fun_g(y(2),y(1));