function neresheni
    hold on
    grid on
    axis([-4,4,-4,4])
    [x0,y0]=ginput(1);
    plot(x0,y0,'r*')
    if x0^2<y0^2
        f=dsolve('x*(Dy)^2-2*y*Dy+x=0','y(a)=b','x');
        s=-4:0.01:4;
        y=subs(f,{'a','b','x'},{x0,y0,s});
        plot(s,-s,s,s,s,y);
    else text(x0+0.1,y0,'No solution')
    end

end

