Draw a polygon. Use Bresenham/Midpoint line drawing algorithm

Write a program to draw a polygon. Use Bresenham/Midpoint line drawing algorithm to draw line.
Note:- Bresenham/Midpoint algorithm used here can be used to draw lines having any slope.

#include<stdio.h>
#include<graphics.h>
#include<math.h>
#define ROUND(a)((int)(a+0.5))
struct point
{
 int x,y;
}pol[10];
int n;
void DDA(int,int,int,int);
void Bresenham(int,int,int,int);
void readpoly()
{
 int i;
 for(i=0;i<n;i++)
 {
  printf("\nEnter the %d co=ordinates ",i+1);
  scanf("%d%d",&pol[i].x,&pol[i].y);
 }
}
void drawpolygonDDA()
{
 int i; 
 for(i=0;i<n;i++)
  DDA(pol[i].x,pol[i].y,pol[(i+1)%n].x,pol[(i+1)%n].y);
}
void drawpolygonBres()
{
 int i; 
 for(i=0;i<n;i++)
 {
  delay(500);
  Bresenham(pol[i].x,pol[i].y,pol[(i+1)%n].x,pol[(i+1)%n].y);
 }
}
int main()
{
 int xa,ya,xb,yb,choice,gd=DETECT,gm=VGAMAX;
 printf("\nEnter the no of sides ");
 scanf("%d",&n);
 readpoly();
 initgraph(&gd,&gm,NULL);
 printf("Press any key to draw traingle using Bresenham Midpoint line 
 drawing algortihm over the triangle drawn using 
 DDA line drawing algorithm");
 drawpolygonDDA();
 getch();
 drawpolygonBres();
 while (!kbhit());
 closegraph();
 return 0;
}
void DDA(int xa,int ya,int xb,int yb)
{
 int dx,dy,steps,k;
 dx=xb-xa;
 dy=yb-ya;
 float x=xa,y=ya,xinc,yinc;
 if(abs(dx)>abs(dy))
  steps=abs(dx);
 else
  steps=abs(dy);
 xinc=dx/(float)steps;
 yinc=dy/(float)steps;
 putpixel(ROUND(x),ROUND(y),5);
 for(k=0;k<steps;k++)
 {
  x+=xinc;
  y+=yinc;
  putpixel(ROUND(x),ROUND(y),5);
 } 
}
void Bresenham(int xa,int ya,int xb,int yb)
{
 int dx=xb-xa,dy=yb-ya,p,twod,twodd;
 float m=2;
 if(dx!=0)
  m=(float)dy/(float)dx;
 if(m>=-1&&m<=1)
 {
  if(xa>xb)
  {
   Bresenham(xb,yb,xa,ya);
   return;
  }
  if(dy<0)
  {
    m=-1;
       dy=-dy;
  }
  else
     m=1;
  putpixel(xa,ya,10);
  p=2*dy-dx;
  twod=2*dy;
  twodd=2*dy-2*dx;
  while(xa<xb)
  {
   xa++;
   if(p<0) 
    p+=twod;
   else
   {
    ya+=m;
    p+=twodd;
   }
   putpixel(xa,ya,10);
  }
 }
 else
 {
  if(ya>yb)
  {
   Bresenham(xb,yb,xa,ya);
   return;
  }
  if(dx<0)
   {
       m=-1;
       dx=-dx;
    }
  else
     m=1;
  putpixel(xa,ya,10);
  p=2*dx-dy;
  twod=2*dx;
  twodd=2*dx-2*dy;
  while(ya<yb)
  {
   ya++;
   if(p<0) 
    p+=twod;
   else
   {
    xa+=m;
    p+=twodd;
   }
   putpixel(xa,ya,10);
  }
 }
}

Window to viewport transformation

#include<stdio.h>
#include<graphics.h>
struct point
{
 int x,y;
} pol[10],wmin,wmax,vmin,vmax;
int n;
float trans[3][3],result[3][1],xandy[3][1];
void readpoly()
{
 int i;
 for(i=0;i<n;i++)
 {
  printf("\nEnter the %d co=ordinates ",i+1);
  scanf("%d%d",&pol[i].x,&pol[i].y);
 }
}
void drawpolygon()
{
 int i; 
 for(i=0;i<n;i++)
  line(pol[i].x,pol[i].y,pol[(i+1)%n].x,pol[(i+1)%n].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<n;i++)
 {
  xandy[0][0]=pol[i].x;
  xandy[1][0]=pol[i].y;
  xandy[2][0]=1;
  mulmat();
  pol[i].x=result[0][0];
  pol[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;
 }
}
int main()
{
 int gd=DETECT,gm=VGAMAX; 
 int choice,dx,dy;
 float sx,sy;
 printf("\nEnter the minx and miny of window ");
 scanf("%d%d",&wmin.x,&wmin.y);
 printf("\nEnter the maxx and maxy of window ");
 scanf("%d%d",&wmax.x,&wmax.y);
 printf("\nEnter the minx and miny of view port ");
 scanf("%d%d",&vmin.x,&vmin.y);
 printf("\nEnter the maxx and maxy of view port ");
 scanf("%d%d",&vmax.x,&vmax.y);
 while(1)
 {
  makeIdentity();
  printf("\nMENU\n-------\n1.POLYGON\n2.EXIT");
  printf("\nEnter the choice ");
  scanf("%d",&choice);
  switch(choice)
  {
   case 1: printf("\nEnter the no of sides ");
    scanf("%d",&n);
    readpoly();
    initgraph(&gd,&gm,NULL);
    outtextxy(200,10,"Window");
    rectangle(wmin.x,wmin.y,wmax.x,wmax.y);
    drawpolygon();
    break;   
   case 2: exit(0);
   default:printf("\nInvalid choice ");
  }
  getch();
  cleardevice();
  outtextxy(200,10,"View Port");
  sx=(float)(vmax.x-vmin.x)/(float)(wmax.x-wmin.x);
  sy=(float)(vmax.y-vmin.y)/(float)(wmax.y-wmin.y); 
  trans[0][0]=sx;
  trans[0][2]=(1-sx)*(float)wmin.x;
  trans[1][1]=sy;  
  trans[1][2]=(1-sy)*(float)wmin.y;
  transformpoly();
  dx=-wmin.x;
  dy=-wmin.y;
  makeIdentity();
  trans[0][2]=dx;
  trans[1][2]=dy;
  transformpoly();
  dx=vmin.x;
  dy=vmin.y;
  makeIdentity();
  trans[0][2]=dx;
  trans[1][2]=dy;
  transformpoly();
  setcolor(RED);
  rectangle(vmin.x,vmin.y,vmax.x,vmax.y);
  drawpolygon(); 
  while (!kbhit());
  closegraph();
 }
 return 0;
}

Analog clock

#include<stdio.h>
#include<graphics.h>
#include<math.h>
#include<time.h>
#define round(a)((int)(a+0.5))
int main()
{
 int gd=DETECT,gm=VGAMAX,sx=0,sy=0,mx=0,my=0,hx=0,hy=0,f=0;
 float hour,min;
 time_t now;
 struct tm *t;
 initgraph(&gd,&gm,"NULL");
 setcolor(2);
 circle(300,300,100);
 setcolor(14);
 circle(300,300,110);
 setcolor(3);
 circle(300,300,120);
 do
 {
  time(&now);
  t=localtime(&now);
  if(f!=0)
  {
   setcolor(BLACK);
   line(300,300,sx,sy);
   line(300,300,mx,my);
   line(300,300,hx,hy);
  }
  hour=(t->tm_hour)%12+((float)t->tm_min/60.0);
  min=t->tm_min+((float)t->tm_sec/60.0);
  sx=300+round(90*cos(3.1412*(270+(6*t->tm_sec))/180));
  sy=300+round(90*sin(3.1412*(270+(6*t->tm_sec))/180));
  mx=300+round(90*cos(3.1412*(270+(6*min))/180));
  my=300+round(90*sin(3.1412*(270+(6*min))/180));
  hx=300+round(50*cos(3.1412*(270+(30*hour))/180));
  hy=300+round(50*sin(3.1412*(270+(30*hour))/180));
  outtextxy(294,380,"6");
  outtextxy(380,300,"3");
  outtextxy(293,220,"12");
  outtextxy(220,300,"9");
  setcolor(YELLOW);
  line(300,300,mx,my);
  setcolor(10);
  line(300,300,hx,hy);
  setcolor(RED);
  line(300,300,sx,sy);
  f=1;
 }while(!kbhit());
 printf("Current time %d:%d:%d",t->tm_hour,t->tm_min,t->tm_sec);
 getch();
 while(!kbhit());
 closegraph();
 return 0;
}

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; }