function homeq
    hold on
    grid on
    axis([-4,4,-10,10]);
    [x0,y0]=ginput(1);
    plot(x0,y0,'ro')
    if x0*y0>0
        f=dsolve('Dy=(y/x)*(1+log(y/x))','y(a)=b','x');
        s=-4:0.01:4;
        g=subs(f,{'a','b','x'},{x0,y0,s});
        plot(s,g)
    else text(x0,y0,'no solution')
    end

end

