Dariusz Zolna pisze:
> Czegoś w tych wzorach brakuje...

Zabrakło nawiasów po znaku dzielenia wokół dr*dr, ale i tak gdzieś
jeszcze jest błąd, bo wychodzą wartości spoza okręgu.

Było:
double x_1 =((D*dy+sgn(dy)*dx*sqrt(r*r*dr*dr-D*D))/dr*dr) + x;
(itd..)

Jest:
double x_1 =((D*dy+sgn(dy)*dx*sqrt(r*r*dr*dr-D*D))/(dr*dr)) + x;
(itd..)

Pozdrawiam,
Dariusz Żołna