-
Notifications
You must be signed in to change notification settings - Fork 0
/
midpointEllipse.cpp
87 lines (77 loc) · 1.76 KB
/
midpointEllipse.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
#include<graphics.h>
#include <bits/stdc++.h>
using namespace std;
void midptellipse(int rx, int ry,
int xc, int yc)
{
float dx, dy, d1, d2, x, y;
x = 0;
y = ry;
// Initial decision parameter of region 1
d1 = (ry * ry) - (rx * rx * ry) +
(0.25 * rx * rx);
dx = 2 * ry * ry * x;
dy = 2 * rx * rx * y;
// For region 1
while (dx < dy)
{
// Print points based on 4-way symmetry
cout << x + xc << " , " << y + yc << endl;
cout << -x + xc << " , " << y + yc << endl;
cout << x + xc << " , " << -y + yc << endl;
cout << -x + xc << " , " << -y + yc << endl;
// Checking and updating value of
// decision parameter based on algorithm
if (d1 < 0)
{
x++;
dx = dx + (2 * ry * ry);
d1 = d1 + dx + (ry * ry);
}
else
{
x++;
y--;
dx = dx + (2 * ry * ry);
dy = dy - (2 * rx * rx);
d1 = d1 + dx - dy + (ry * ry);
}
}
// Decision parameter of region 2
d2 = ((ry * ry) * ((x + 0.5) * (x + 0.5))) +
((rx * rx) * ((y - 1) * (y - 1))) -
(rx * rx * ry * ry);
// Plotting points of region 2
while (y >= 0)
{
cout << x + xc << " , " << y + yc << endl;
putpixel(x+xc,y+yc,RED);
cout << -x + xc << " , " << y + yc << endl;
putpixel(-x+xc,y+yc,RED);
cout << x + xc << " , " << -y + yc << endl;
putpixel(x+xc,-y+yc,RED);
putpixel(-x+xc,-y+yc,RED);
cout << -x + xc << " , " << -y + yc << endl;
if (d2 > 0)
{
y--;
dy = dy - (2 * rx * rx);
d2 = d2 + (rx * rx) - dy;
}
else
{
y--;
x++;
dx = dx + (2 * ry * ry);
dy = dy - (2 * rx * rx);
d2 = d2 + dx - dy + (rx * rx);
}
}
}
int main()
{ int gd=DETECT,gm;
initgraph(&gd,&gm,"n");
midptellipse(100, 75, 200, 125);
getch();
return 0;
}