Draw the following figure(1)

#include<stdio.h>
#include<graphics.h>
struct point
{
 int x,y;
}pol[3],polnew[3],circ,insid[12];
int n,r;
float trans[3][3],result[3][1],xandy[3][1];
void plot(int,int,int,int);
void boundaryfill4(int,int,int,int);
void midpoint(int xc,int yc,int r)
{
 int x=0;
 int y=r;
 int p=1-r;
 plot(xc,yc,x,y);
 while(x<y)
 {
  x++;
  if(p<0)
   p+=2*x+1;
  else
  {
   y--;   
   p+=2*(x-y)+1;
  }
  plot(xc,yc,x,y);
 }
}
void plot(int xc,int yc,int x,int y)
{
 putpixel(xc+x,yc+y,10);
 putpixel(xc-x,yc+y,10);
 putpixel(xc+x,yc-y,10);
 putpixel(xc-x,yc-y,10);
 putpixel(xc+y,yc+x,10);
 putpixel(xc-y,yc+x,10);
 putpixel(xc+y,yc-x,10);
 putpixel(xc-y,yc-x,10);
}
void drawpol()
{
 int i; 
 for(i=0;i<2;i++)
  line(pol[i].x,pol[i].y,pol[(i+1)].x,pol[(i+1)].y);
}
void drawnewpol()
{
 int i; 
 for(i=0;i<2;i++)
  line(polnew[i].x,polnew[i].y,polnew[(i+1)].x,polnew[(i+1)].y);
}
void mulmat()
{
 int i,j,k;
 for(i=0;i<3;i++)
 {
  for(j=0;j<1;j++)
  {
   result[i][j]=0.0;
   for(k=0;k<3;k++)
    result[i][j]+=trans[i][k]*xandy[k][j];
  }
 }
}
void transformpoint(int x,int y)
{
 xandy[0][0]=x;
 xandy[1][0]=y;
 xandy[2][0]=1;
 mulmat();
}
void transformpoly()
{
 int i; 
 for(i=0;i<3;i++)
 {
  xandy[0][0]=pol[i].x;
  xandy[1][0]=pol[i].y;
  xandy[2][0]=1;
  mulmat();
  polnew[i].x=result[0][0];
  polnew[i].y=result[1][0];
 }
}
void makeIdentity()
{
 int i,j;
 for(i=0;i<3;i++)
 {
  for(j=0;j<3;j++)
   trans[i][j]=0.0;
  trans[i][i]=1.0;
 }
}
void rotate(int rx,int ry,int an)
{
 float ang;
 ang=(float)(an*(3.14/180));
 trans[0][0]=cos(ang);
 trans[0][1]=-sin(ang);
 trans[0][2]=rx*(1-cos(ang))+ry*sin(ang);
 trans[1][0]=sin(ang);
 trans[1][1]=cos(ang);
 trans[1][2]=ry*(1-cos(ang))-rx*sin(ang);
 
}
void boundaryfill4(int x,int y,int new,int bdry)
{
 int current=getpixel(x,y);
 if((current!=bdry)&&(current!=new))
 {
  putpixel(x,y,new);
  boundaryfill4(x+1,y,new,bdry);
  boundaryfill4(x-1,y,new,bdry);
  boundaryfill4(x,y+1,new,bdry);
  boundaryfill4(x,y-1,new,bdry);
 }
}
int main()
{
 int i,j,angle=0,gd=DETECT,gm=VGAMAX,ch;
 int color[]={1,14,3,4,5,6,7,8,9,11,12,13};
 printf("Enter the orgin ");
 scanf("%d%d",&circ.x,&circ.y);
 printf("Enter the radius ");
 scanf("%d",&r);
 do{
  printf("\nMENU\n-------\n1.Enter inside points to fill\n
  2.Take the inside points automatically");
  printf("\nEnter ur choice  ");
  scanf("%d",&ch);
 }while(ch>2&&ch<1);
 if(ch==1)
 {
  printf("\nEnter the inside point of circle ");
  scanf("%d%d",&insid[0].x,&insid[0].y);
  for(i=1;i<12;i++)
  {
   printf("\nEnter the inside point %d ",i);
   scanf("%d%d",&insid[i].x,&insid[i].y);
  }
 }
 initgraph(&gd,&gm,"NULL");
 setcolor(10);
 midpoint(circ.x,circ.y,r);
 makeIdentity();
 j=360/11;
 pol[0].x=circ.x;
 pol[0].y=circ.y+r;
 pol[1].x=pol[0].x;
 pol[1].y=pol[0].y+50;
 rotate(circ.x,circ.y,j);
 transformpoint(pol[0].x,pol[0].y);
 pol[2].x=result[0][0];
 pol[2].y=result[1][0];
 rotate(circ.x,circ.y,360/22);
 transformpoint(pol[1].x,pol[1].y);
 pol[1].x=result[0][0];
 pol[1].y=result[1][0];
 if(ch==2)
 {
  insid[0].x=circ.x;
  insid[0].y=circ.y;
  insid[1].x=(pol[0].x+pol[1].x+pol[2].x)/3;
  insid[1].y=(pol[0].y+pol[1].y+pol[2].y)/3;
 }
 boundaryfill4(insid[0].x,insid[0].y,color[0],10);
 drawpol();
 boundaryfill4(insid[1].x,insid[1].y,color[1],10);
 for(i=0;i<10;i++)
 {
  angle+=j;
  rotate(circ.x,circ.y,angle);
  transformpoly();
  drawnewpol();
  if(ch==2)
  {
   insid[i+2].x=(polnew[0].x+polnew[1].x+polnew[2].x)/3;
   insid[i+2].y=(polnew[0].y+polnew[1].y+polnew[2].y)/3;
  }
  boundaryfill4(insid[i+2].x,insid[i+2].y,color[i+2],10);
 }
 while(!kbhit());
 closegraph();
 return 0;
}

Rotate a triangle along the circumference of a circle, also rotate the triangle with reference to its centroid

Rotate a triangle along the circumference of a circle,  also rotate the triangle with reference to its centroid. 
Use midpoint algorithm for circle drawing. #include<stdio.h> #include<graphics.h> struct point { int x,y; }pol[3],polnew[3],circ,cent; int n,r; float trans[3][3],result[3][1],xandy[3][1]; void plot(int,int,int,int); void midpoint(int xc,int yc,int r) { int x=0; int y=r; int p=1-r; plot(xc,yc,x,y); while(x<y) { x++; if(p<0) p+=2*x+1; else { y--; p+=2*(x-y)+1; } plot(xc,yc,x,y); } } void plot(int xc,int yc,int x,int y) { putpixel(xc+x,yc+y,10); putpixel(xc-x,yc+y,10); putpixel(xc+x,yc-y,10); putpixel(xc-x,yc-y,10); putpixel(xc+y,yc+x,10); putpixel(xc-y,yc+x,10); putpixel(xc+y,yc-x,10); putpixel(xc-y,yc-x,10); } void drawpolygon() { int i; for(i=0;i<3;i++) line(pol[i].x,pol[i].y,pol[(i+1)%3].x,pol[(i+1)%3].y); } void drawnewpolygon() { int i; for(i=0;i<3;i++) line(polnew[i].x,polnew[i].y,polnew[(i+1)%3].x,polnew[(i+1)%3].y); } void mulmat() { int i,j,k; for(i=0;i<3;i++) { for(j=0;j<1;j++) { result[i][j]=0.0; for(k=0;k<3;k++) result[i][j]+=trans[i][k]*xandy[k][j]; } } } void transformpoly() { int i; for(i=0;i<3;i++) { xandy[0][0]=pol[i].x; xandy[1][0]=pol[i].y; xandy[2][0]=1; mulmat(); polnew[i].x=result[0][0]; polnew[i].y=result[1][0]; } } void transformnewpoly() { int i; for(i=0;i<3;i++) { xandy[0][0]=polnew[i].x; xandy[1][0]=polnew[i].y; xandy[2][0]=1; mulmat(); polnew[i].x=result[0][0]; polnew[i].y=result[1][0]; } } void makeIdentity() { int i,j; for(i=0;i<3;i++) { for(j=0;j<3;j++) trans[i][j]=0.0; trans[i][i]=1.0; } } void rotate(int rx,int ry,int an) { float ang; ang=(float)(an*(3.14/180)); trans[0][0]=cos(ang); trans[0][1]=-sin(ang); trans[0][2]=rx*(1-cos(ang))+ry*sin(ang); trans[1][0]=sin(ang); trans[1][1]=cos(ang); trans[1][2]=ry*(1-cos(ang))-rx*sin(ang); } int main() { int i,j,gd=DETECT,gm=VGAMAX; printf("\nEnter the sides of triangle "); for(i=0;i<3;i++) { printf("\nEnter the %d co=ordinates ",i+1); scanf("%d%d",&pol[i].x,&pol[i].y); } printf("Enter the orgin "); scanf("%d%d",&circ.x,&circ.y); cent.x=(pol[0].x+pol[1].x+pol[2].x)/3; cent.y=(pol[0].y+pol[1].y+pol[2].y)/3; r=sqrt((circ.x-cent.x)*(circ.x-cent.x)+(circ.y-cent.y)*(circ.y-cent.y)); initgraph(&gd,&gm,"NULL"); midpoint(circ.x,circ.y,r); drawpolygon(); putpixel(cent.x,cent.y,5); makeIdentity(); i=1; j=1; do{ midpoint(circ.x,circ.y,r); drawnewpolygon(); rotate(cent.x,cent.y,i); transformpoly(); rotate(circ.x,circ.y,j); transformnewpoly(); delay(40); i=(i+10)%360; j=(j+5)%360; cleardevice(); }while(!kbhit()); closegraph(); return 0; }

Clip a polygon and fill the area.

#include<stdio.h> #include<graphics.h> struct point { int x; int y; }wmin,wmax,p1,p2,ipt,pin[50],pout[50]; int np=0,n,gd=DETECT,gm=VGAMAX; void newpin(); void boundaryfill4(int,int,int,int); void drawpolygon(); void showpoints(); void clippolygon() { float m=0.0; int i,j; for(i=0,j=0;i<n;i++) { m=0.0; p1=pin[i]; p2=pin[(i+1)%n]; if(p1.x>=wmin.x&&p2.x>=wmin.x) { pout[j]=p2; j++; } else if(p2.x==wmin.x) { pout[j]=p2; j++; } else if(p1.x<wmin.x&&p2.x>wmin.x||p1.x>wmin.x&&p2.x<wmin.x) { if(p1.x!=p2.x) m=(float)(p1.y-p2.y)/(float)(p1.x-p2.x); ipt.x=wmin.x; ipt.y=p1.y+(wmin.x-p1.x)*m; if(p1.x<wmin.x) { pout[j]=ipt; pout[j+1]=p2; j=j+2; } else { p2=ipt; pout[j]=ipt; j++; } } } n=j; newpin(); outtextxy(220,10,"After left clip"); drawpolygon(); while (!kbhit()); closegraph(); initgraph(&gd,&gm,"NULL"); for(i=0,j=0;i<n;i++) { m=0.0; p1=pin[i]; p2=pin[(i+1)%n]; if(p1.x<=wmax.x&&p2.x<=wmax.x) { pout[j]=p2; j++; } else if(p2.x==wmax.x) { pout[j]=p2; j++; } else if(p1.x>wmax.x&&p2.x<wmax.x||p1.x<wmax.x&&p2.x>wmax.x) { if(p1.x!=p2.x) m=(float)(p1.y-p2.y)/(float)(p1.x-p2.x); ipt.x=wmax.x; ipt.y=p1.y+(wmax.x-p1.x)*m; if(p1.x>wmax.x) { pout[j]=ipt; pout[j+1]=p2; j=j+2; } else { p2=ipt; pout[j]=ipt; j++; } } } n=j; newpin(); outtextxy(220,10,"After right clip"); drawpolygon(); while (!kbhit()); closegraph(); initgraph(&gd,&gm,"NULL"); for(i=0,j=0;i<n;i++) { m=0.0; p1=pin[i]; p2=pin[(i+1)%n]; if(p1.y>=wmin.y&&p2.y>=wmin.y) { pout[j]=p2; j++; } else if(p2.y==wmin.y) { pout[j]=p2; j++; } else if(p1.y<wmin.y&&p2.y>wmin.y||p1.y>wmin.y&&p2.y<wmin.y) { if(p1.x!=p2.x) m=(float)(p1.y-p2.y)/(float)(p1.x-p2.x); ipt.y=wmin.y; if(m!=0) ipt.x=p1.x+(float)(wmin.y-p1.y)*(float)(1/m); else ipt.x=p1.x; if(p1.y<wmin.y) { pout[j]=ipt; pout[j+1]=p2; j=j+2; } else { p2=ipt; pout[j]=ipt; j++; } } } n=j; newpin(); outtextxy(220,10,"After bottom clip"); drawpolygon(); while (!kbhit()); closegraph(); initgraph(&gd,&gm,"NULL"); for(i=0,j=0;i<n;i++) { m=0.0; p1=pin[i]; p2=pin[(i+1)%n]; if(p1.y<=wmax.y&&p2.y<=wmax.y) { pout[j]=p2; j++; } else if(p2.y==wmax.y) { pout[j]=p2; j++; } else if(p1.y>wmax.y&&p2.y<wmax.y||p1.y<wmax.y&&p2.y>wmax.y) { if(p1.x!=p2.x) m=(float)(p1.y-p2.y)/(float)(p1.x-p2.x); ipt.y=wmax.y; if(m!=0) ipt.x=p1.x+(float)(wmax.y-p1.y)*(float)(1/m); else ipt.x=p1.x; if(p1.y>wmax.y) { pout[j]=ipt; pout[j+1]=p2; j=j+2; } else { p2=ipt; pout[j]=ipt; j++; } } } n=j; newpin(); outtextxy(220,10,"clipped Polygon"); showpoints(); drawpolygon(); } void newpin() { int i; for(i=0;i<n;i++) pin[i]=pout[i]; } void showpoints() { int j; printf("Vertex"); for(j=0;j<n;j++) printf("\n%d.(%d,%d)",j+1,pin[j].x,pin[j].y); } void drawpolygon() { int j; setcolor(10); rectangle(wmin.x,wmin.y,wmax.x,wmax.y); setcolor(6); for(j=0;j<n;j++) line(pin[j].x,pin[j].y,pin[(j+1)%n].x,pin[(j+1)%n].y); } void boundaryfill4(int x,int y,int new,int bdry) { int current=getpixel(x,y); if((current!=bdry)&&(current!=new)) { putpixel(x,y,new); boundaryfill4(x+1,y,new,bdry); boundaryfill4(x-1,y,new,bdry); boundaryfill4(x,y+1,new,bdry); boundaryfill4(x,y-1,new,bdry); } } int main() { int c,j; int newcolor,ix,iy; printf("Enter the window minimum coordinates"); scanf("%d%d",&wmin.x,&wmin.y); printf("Enter the window max coordinates"); scanf("%d%d",&wmax.x,&wmax.y); printf("Enter the no of sides in polygon:\n"); scanf("%d",&n); printf("Enter the coordinates(x,y)\n"); for(j=0;j<n;j++) { printf("Enter the coordinate %d\n",j+1); scanf("%d%d",&pin[j].x,&pin[j].y); } printf("\nEnter the fill color "); scanf("%d",&newcolor); printf("\nEnter the inner point "); scanf("%d%d",&ix,&iy); initgraph(&gd,&gm,"NULL"); outtextxy(220,10,"Orginal Polygon"); drawpolygon(); while (!kbhit()); closegraph(); initgraph(&gd,&gm,"NULL"); clippolygon(); boundaryfill4(ix,iy,newcolor,6); while(!kbhit()); closegraph(); return 0; }

Bezier, Bspline, Cubic Spline Curves

#include<stdio.h> #include<graphics.h> #include<math.h> #define ROUND(a)((int)(a+0.5)) int gd=DETECT,gm=VGAMAX; struct point { double x,y,z; }controls[10],output[200]; int outCount; void computeCoefBezier(int n,int *c) { int k,i; for(k=0;k<=n;k++) { c[k]=1; for(i=n;i>=k+1;i--) c[k]*=i; for(i=n-k;i>=2;i--) c[k]/=i; } } void computePointBezier(float u,int nControls,int *c,struct point *calcxyz) { int k,n=nControls-1; float blend; calcxyz->x=0; calcxyz->y=0; calcxyz->z=0; for(k=0;k<nControls;k++) { blend=c[k]*pow(u,k)*pow(1-u,n-k); calcxyz->x+=controls[k].x*blend; calcxyz->y+=controls[k].y*blend; calcxyz->z+=controls[k].z*blend; } } void bezier(int nControls) { int *c=(int *)malloc(nControls*sizeof(int)); int i; computeCoefBezier(nControls-1,c); for(i=0;i<=outCount;i++) computePointBezier(i/(float)outCount,nControls,c,&output[i]); } double blend(int k, int t, int *u, double v) {   double value;   if (t==1)   {     if ((u[k]<=v) && (v<u[k+1]))       value=1;     else       value=0;   }   else   {     if ((u[k+t-1]==u[k]) && (u[k+t]==u[k+1]))       value = 0;     else      if (u[k+t-1]==u[k])       value = (u[k+t] - v) / (u[k+t] - u[k+1]) * blend(k+1, t-1, u, v);     else      if (u[k+t]==u[k+1])       value = (v - u[k]) / (u[k+t-1] - u[k]) * blend(k, t-1, u, v);      else       value = (v - u[k]) / (u[k+t-1] - u[k]) * blend(k, t-1, u, v) +       (u[k+t] - v) / (u[k+t] - u[k+1]) * blend(k+1, t-1, u, v);   }   return value; } void compute_intervals(int *u, int n, int t) {   int j;   for (j=0; j<=n+t; j++)   {     if (j<t)       u[j]=0;     else     if ((t<=j) && (j<=n))       u[j]=j-t+1;     else     if (j>n)       u[j]=n-t+2;   } } void compute_point(int *u, int n, int t, double v,struct point *calc) {   int k;   double temp;   calc->x=0;   calc->y=0;   calc->z=0;   for (k=0; k<=n; k++)   {     temp = blend(k,t,u,v);     calc->x = calc->x + (controls[k]).x * temp;     calc->y = calc->y + (controls[k]).y * temp;     calc->z = calc->z + (controls[k]).z * temp;   } } void bspline(int n, int t) {   int *u=(int *)malloc((n+t+1)*sizeof(int));   double increment,interval;   int i;   compute_intervals(u, n, t);   increment=(double) (n-t+2)/(outCount-1);   interval=0;   for (i=0; i<outCount; i++)   {     compute_point(u, n, t, interval,&output[i]);     interval=interval+increment;   }   output[outCount].x=controls[n].x;   output[outCount].y=controls[n].y;   output[outCount].z=controls[n].z; } void drawCurve() { int i; for(i=0;i<outCount;i++) line(output[i].x,output[i].y,output[i+1].x,output[i+1].y); } void cubic(int n) { int i,j; float u,f1,f2,f3,f4,px,py,px1,py1; initgraph(&gd,&gm,"NULL"); setcolor(2); for(i=0;i<n-1;i++) { for(j=0;j<=40;j++) { u=j/(float)40; f1=(2*u*u*u)-(3*u*u)+1; f2=-(2*u*u*u)+(3*u*u); f3=(u*u*u)-(2*u*u)+u; f4=(u*u*u)-(u*u); px=((float)controls[i].x*f1)+((float)controls[i+1].x*f2)+((float)controls[i].x*f3)+((float)controls[i+1].x*f4); py=((float)controls[i].y*f1)+((float)controls[i+1].y*f2)+((float)controls[i].y*f3)+((float)controls[i+1].y*f4); if(j>0) line(ROUND(px1),ROUND(py1),ROUND(px),ROUND(py)); px1=px; py1=py; } } setcolor(5); for(i=0;i<n;i++) circle(controls[i].x,controls[i].y,2); while(!kbhit()); closegraph(); } int main() { int i,n,ch,t=3; printf("\nEnter the no of controls points"); scanf("%d",&n); printf("\nEnter the cordinates"); for(i=0;i<n;i++) { printf("x%d,y%d=",i,i); scanf("%lf%lf",&controls[i].x,&controls[i].y); } outCount=n*15; while(1) { printf("\nMENU\n-------\n1.Bspline\n2.Bezier\n3.Cubic Spline\n4.EXIT"); printf("\nEnter the choice "); scanf("%d",&ch); switch(ch) { case 1: if(n<3) { printf("\nBspline doesnt give good result"); exit(0); } bspline(n-1,t); break; case 2: bezier(n); break; case 3: cubic(n); break; case 4: exit(0); default:printf("\nInvalid Entry"); } if(ch==1||ch==2) { initgraph(&gd,&gm,"NULL"); setcolor(2); drawCurve(); setcolor(5); for(i=0;i<n;i++) circle(controls[i].x,controls[i].y,2); while(!kbhit()); closegraph(); } } return 0; }

Animation of some real world objects(Wind Mill)

#include<stdio.h> #include<graphics.h> #include<math.h> #define ROUND(a)((int)(a+0.5)) void paraCircle(int,int,int,int,int); void plotellipse(int,int,int,int); void paraCircle(int xc,int yc,int r,int t,int c) { int i,x,y; for(i=0;i<10;i++) { x=xc+ROUND(r*cos(3.1412*(i+t)/180)); y=yc+ROUND(r*sin(3.1412*(i+t)/180)); line(xc,yc,x,y); putpixel(x,y,c); } line(xc,yc,x,y); for(i=120;i<130;i++) { x=xc+ROUND(r*cos(3.1412*(i+t)/180)); y=yc+ROUND(r*sin(3.1412*(i+t)/180)); line(xc,yc,x,y); putpixel(x,y,c); } line(xc,yc,x,y); for(i=240;i<250;i++) { x=xc+ROUND(r*cos(3.1412*(i+t)/180)); y=yc+ROUND(r*sin(3.1412*(i+t)/180)); line(xc,yc,x,y); putpixel(x,y,c); } line(xc,yc,x,y); } void midpointellipse(int xc,int yc,int rx,int ry) { int rx2=rx*rx; int ry2=ry*ry; int tworx2=2*rx2; int twory2=2*ry2; int p,x=0,y=ry,px=0,py=tworx2*y; plotellipse(xc,yc,x,y); p=ROUND(ry2-(rx2*ry)+(0.25*rx2)); while(px<py) { x++; px+=twory2; if(p<0) p+=ry2+px; else { y--; py-=tworx2; p+=ry2+px-py; } plotellipse(xc,yc,x,y); } p=ROUND(ry2*(x+.5)*(x+.5)+rx2*(y-1)*(y-1)-rx2*ry2); while(y>0) { y--; py-=tworx2; if(p>0) p+=rx2-py; else { x++; px+=twory2; p+=rx2-py+px; } plotellipse(xc,yc,x,y); } } void plotellipse(int xc,int yc,int x,int y) { putpixel(xc+x,yc+y,7); putpixel(xc-x,yc+y,7); putpixel(xc+x,yc-y,7); putpixel(xc-x,yc-y,7); } int main() { int theta=0,speed=500; int gd=DETECT,gm=VGAMAX; initgraph(&gd,&gm,"NULL"); midpointellipse(300,350,30,5); line(290,350,300,200); line(310,350,300,200); paraCircle(300,200,50,theta,7); outtextxy(200,30,"Press any key to continue"); getch(); cleardevice(); midpointellipse(300,350,30,5); while(!kbhit()){ setcolor(7); line(290,350,300,200); line(310,350,300,200); paraCircle(300,200,50,theta,7); delay(speed); setcolor(0); paraCircle(300,200,50,theta,0); theta=theta+10; if(speed>=100) speed-=30; } closegraph(); return 0; }

Circumcircle of a polygon

#include<stdio.h> #include<graphics.h> #include<math.h> int main() { int n,i,x[10],y[10],a[10],b[10]; float m[10],c[10],xc,yc,r; printf("Enter the no of sides: "); scanf("%d",&n); for(i=0;i<n;i++) { printf(" Enter the x and y cordinates: "); scanf("%d %d",&x[i],&y[i]); } x[n]=x[0]; y[n]=y[0]; int gd=DETECT,gm=VGAMAX; initgraph(&gd,&gm,"NULL"); for(i=0;i<n;i++) { line(x[i],y[i],x[i+1],y[i+1]); delay(500); } for(i=0;i<n;i++) { a[i]=(x[i]+x[i+1])/2; b[i]=(y[i]+y[i+1])/2; } for(i=0;i<n;i++) {   if(y[i]!=y[i+1]) {   m[i]=-(float)(x[i]-x[i+1])/(y[i]-y[i+1]);   c[i]=b[i]-m[i]*a[i];   }   else   m[i]=0; } i=0; if(y[i]==y[i+1]) { xc=a[i]; yc=(float)(m[i+1]*xc+c[i+1]); } else if(y[i+1]==y[i+2]) {   xc=a[i+1];   yc=(float)(m[i+2]*xc+c[i+2]); } else { xc=(c[i+1]-c[i])/(m[i]-m[i+1]); yc=m[i]*xc+c[i]; } r=sqrt((xc-x[1])*(xc-x[1])+(yc-y[1])*(yc-y[1])); putpixel((int)xc,(int)yc,5); printf("Circumcenter--> (%f, %f)\nCircumradius--> %f",xc,yc,r); setcolor(9); circle((int)xc,(int)yc,r); while(!kbhit()); closegraph(); return 0; }

Bouncing Ball

#include<stdio.h> #include<graphics.h> #include<math.h> #define ROUND(a)((int)(a+0.5)) void paraCircle(int,int,int,int); int base=400; int color=15; void parabola() { int xo=70,x=0,yo=base,theta=0; float vh=10.0,vy=80.0; int vhit=0; int y=0; float t; int dly=40; int xmax=getmaxx(); outtextxy(20,50,"Press any key to throw"); paraCircle(xo,yo,10,theta); line(0,base+10,xmax-20,base+10); line(xmax-20,0,xmax-20,base+10); getch(); do { theta=0; paraCircle(xo,yo,10,theta); for(t=1;;t+=.2) { x=ROUND(xo+(vh*t)); y=ROUND(yo-(vy*t)+5*(t*t)); line(0,base+10,xmax-20,base+10); line(xmax-20,0,xmax-20,base+10); paraCircle(x,y,10,theta); delay(dly); cleardevice(); if(x==xmax-30) { vh=-vh; vhit=1; } if(vhit==0) theta+=10; else  theta-=10; if(y==base) break; } vy-=5; if(vhit==0) vh-=1; else if(vhit==1) vh+=1; dly-=4; paraCircle(x,y,10,theta); xo=x; }while(vy>0&&vh!=0); line(0,base+10,xmax-20,base+10); line(xmax-20,0,xmax-20,base+10); } void paraCircle(int xc,int yc,int r,int t) { int i,x,y; for(i=0;i<90;i++) { x=xc+ROUND(r*cos(3.1412*(i+t)/180)); y=yc+ROUND(r*sin(3.1412*(i+t)/180)); putpixel(x,y,5); } for(i=90;i<180;i++) { x=xc+ROUND(r*cos(3.1412*(i+t)/180)); y=yc+ROUND(r*sin(3.1412*(i+t)/180)); putpixel(x,y,2); } for(i=180;i<270;i++) { x=xc+ROUND(r*cos(3.1412*(i+t)/180)); y=yc+ROUND(r*sin(3.1412*(i+t)/180)); putpixel(x,y,12); } for(i=270;i<360;i++) { x=xc+ROUND(r*cos(3.1412*(i+t)/180)); y=yc+ROUND(r*sin(3.1412*(i+t)/180)); putpixel(x,y,7); } } int main() { int gd=DETECT,gm=VGAMAX; initgraph(&gd,&gm,"NULL"); parabola(); while(!kbhit()); closegraph(); return 0; }