The original equation of the ellipse:
f(x, y) = b2x2 + a2y2 - a2b2 = 0
- Initial point: (0, b)
- Decision parameter: p1k = b2 - a2 b + 1/4 a2
- Next point:
- If f(x, y) < 0:
- The next point will be (Xk+1, yk)
- p1k+1 = p1k + 2 b2 Xk+1 + b2
- Else:
- The next point will be (Xk+1, yk-1)
- p1k+1 = p1k + 2 b2 Xk+1 + b2 - 2 a2 yk+1
- Continue until 2b2 Xk >= 2 a2 yk
- Initial point: (Ending point of Region 1)
- Decision parameter: p2k = (xk + 1/2)2 b2 + (yk - 1)2 a2 - a2b2
- Next point:
- If f(x, y) > 0:
- The next point will be (xk, yk-1)
- p2k+1 = p2k + a2 - 2 a2 yk-1
- Else:
- The next point will be (xk+1, yk-1)
- p2k+1 = p2k + a2 - 2 a2 yk-1+ 2 b2 xk+1