function homeq1
    hold on
    grid on
    axis([-4,4,-10,10]);
    [x0,y0]=ginput(1);
    plot(x0,y0,'ro')
    
    f=dsolve('y*Dy=-x','y(a)=b','x');
    d=sqrt(x0^2+y0^2);
    s=-d:d/100:d;
    g=subs(f,{'a','b','x'},{x0,y0,s});
    plot(s,g)
   
end

