function homeq1
    hold on
    grid on
    axis([-4,4,-10,10]);
    [x0,y0]=ginput(1);
    plot(x0,y0,'ro')
    if x0*y0>0
       [T,y]=ode45(@ff,[x0,x0+10],y0);
       [T1,y1]=ode45(@ff,[x0,x0-10],y0);
       plot(T,y)
       plot(T1,y1)
    else text(x0,y0,'no solution')
    end
    function z=ff(x,y)
        z=(y./x).*(1+log(y./x));
    end
end

