Draw a Sierpinski triangle using fractal geometric method

#include<stdio.h> 
#include<graphics.h> 
struct point {  
    int x,y; 
}pol[3],polnew[3],wmin,wmax,vmin,vmax; 
float sx=0.5,sy=0.5; 
int n=3,max; 
float trans[3][3],result[3][1],xandy[3][1]; 
float minside() 
{  
    int i;  
    float a,s;  
    i=0;
    s=sqrt((pol[i].x-pol[(i+1)%3].x)*(pol[i].x-pol[(i+1)%3].x)+
    (pol[i].y-pol[(i+1)%3].y)*(pol[i].y-pol[(i+1)%3].y)); 
    for(i=1;i<3;i++)  {
        a=sqrt((pol[i].x-pol[(i+1)%3].x)*(pol[i].x-pol[(i+1)%3].x)+
       (pol[i].y-pol[(i+1)%3].y)*(pol[i].y-pol[(i+1)%3].y));
        if(a<s)    
            s=a;  
    }  
    return s; 
}
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(struct point p[]) {  
    int i;   
    for(i=0;i<n;i++)   
        line(p[i].x,p[i].y,p[(i+1)%n].x,p[(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(struct point p1[],struct point p2[]) {  
    int i;
    for(i=0;i<n;i++)  
    {   
        xandy[0][0]=p1[i].x; 
        xandy[1][0]=p1[i].y;
        xandy[2][0]=1;   
        mulmat();   
        p2[i].x=result[0][0]; 
        p2[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 sierp(struct point pp[],int level) {
    struct point p1[3],p2[3],p3[3];
    int dx,dy,i;
    level++;  
    if(level==max) 
        return;  
    makeIdentity();
    trans[0][0]=sx;
    trans[0][2]=(1-sx)*(float)pp[0].x;
    trans[1][1]=sy; 
    trans[1][2]=(1-sy)*(float)pp[0].y;
    transformpoly(pp,p1); 
    drawpolygon(p1);
    makeIdentity(); 
    dx=-p1[1].x;
    dy=-p1[1].y;
    trans[0][2]=dx;
    trans[1][2]=dy;  
    transformpoly(p1,p2); 
    makeIdentity();  
    dx=pp[1].x; 
    dy=pp[1].y; 
    trans[0][2]=dx;
    trans[1][2]=dy;
    transformpoly(p2,p2);
    drawpolygon(p2);
    makeIdentity();
    dx=-p1[2].x; 
    dy=-p1[2].y;
    trans[0][2]=dx;
    trans[1][2]=dy;
    transformpoly(p1,p3);
    makeIdentity();
    dx=pp[2].x;
    dy=pp[2].y; 
    trans[0][2]=dx;
    trans[1][2]=dy;
    transformpoly(p3,p3); 
    drawpolygon(p3);
    sierp(p1,level);
    sierp(p2,level);
    sierp(p3,level); 
} 
int main() 
{  
    int gd=DETECT,gm=VGAMAX;   
    int choice,dx,dy,i;
    float small;
    readpoly();
    initgraph(&gd,&gm,NULL);  
    drawpolygon(pol);  
    setcolor(RED);
    getch(); 
    small=minside(); 
    max=2+(int)small/50; 
    sierp(pol,0);
    while (!kbhit());
    closegraph(); 
    return 0; 
}

Draw a piechart

#include<stdio.h>
#include<graphics.h>
#include<math.h>
#define ROUND(a)((int)(a+0.5))
int p[10];
float trans[3][3],result[3][1],xandy[3][1];
void plot(int,int,int,int);
void boundaryfill4(int,int,int,int);
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),10);
 for(k=0;k<steps;k++)
 {
  x+=xinc;
  y+=yinc;
  putpixel(ROUND(x),ROUND(y),10);
 } 
}
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 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);
 }
}
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 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 gd=DETECT,gm=VGAMAX,n,i,tot=0,r=100,rx,ry,x,y,midx,midy,prx,pry;
 float an1=0.0,an=0.0;
 printf("enter the no. of data");
 scanf("%d",&n);
 do{
  for(i=0;i<n;i++)
  {
   printf("\nEnter the data  ");
   scanf("%d",&p[i]);
   tot+=p[i];
  }
  if(tot!=100)
  {
   printf("Error in data. Enter it again  ");
   tot=0;
  }
 }while(tot!=100);
 initgraph(&gd,&gm,"NULL");
 rx=getmaxx()/2;
 ry=getmaxy()/2;
 x=rx;
 y=r+ry;
 midx=rx;
 midy=(y+ry)/2;
 midpoint(rx,ry,r);
 DDA(rx,ry,x,y);
 for(i=0;i<n;i++)
 {
  makeIdentity();
  an1=an+(float)(p[i]*3.6)/2;
  an+=(float)p[i]*3.6;
  rotate(rx,ry,an);
  transformpoint(x,y);
  DDA(rx,ry,result[0][0],result[1][0]);
  makeIdentity();
  rotate(rx,ry,an1);
  transformpoint(midx,midy);
  boundaryfill4(result[0][0],result[1][0],i+1,10); 
 }
 while(!kbhit());
 closegraph();
 return 0;
}

Draw a Sierpinski triangle

#include<stdio.h>
#include<graphics.h>
#include<math.h>
#define pi 3.142857
#define ROUND(a)((int)(a+0.5))
int max;
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);
 } 
}
float minside(int x1,int y1,int x2,int y2,int x3,int y3)
{
 float a,b,c,s;
 a=sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
 b=sqrt((x2-x3)*(x2-x3)+(y2-y3)*(y2-y3));
 c=sqrt((x1-x3)*(x1-x3)+(y1-y3)*(y1-y3));
 if(a>b)
  s=a;
 else
  s=b;
 if(s>c)
  s=c;
 return s;
}
void triangle(int x1,int y1,int x2,int y2,int x3,int y3,int i)
{ 

 float newx1,newy1,newx2,newy2,newx3,newy3;
 i++;
 if(i==max)
  return;
 else
 {
  DDA(x1,y1,x2,y2);
  DDA(x2,y2,x3,y3);
  DDA(x3,y3,x1,y1);
  newx1=(x1+x2)/2 ;
  newy1=(y1+y2)/2;
  newx2=(x2+x3)/2 ;
  newy2=(y2+y3)/2;
  newx3=(x1+x3)/2 ;
  newy3=(y1+y3)/2;
  triangle(x1,y1,newx1,newy1,newx3,newy3,i);
  triangle(x2,y2,newx1,newy1,newx2,newy2,i);
  triangle(x3,y3,newx2,newy2,newx3,newy3,i);
 }
}
int main()
{
 int gd=DETECT,gm=VGAMAX;
 int x1,y1,x2,y2,x3,y3;
 float small;
 printf("\n Enter the 1st coordinates of the triangle");
 scanf("%d%d",&x1,&y1);
 printf("\n enter the 2nd coordinates of the triangle");
 scanf("%d%d",&x2,&y2);
 printf("\n enter the 3rd coordinates of the triangle");
 scanf("%d%d",&x3,&y3);
 initgraph(&gd,&gm,"null");
 small=minside(x1,y1,x2,y2,x3,y3);
 max=2+(int)small/20;
 triangle(x1,y1,x2,y2,x3,y3,0);
 while (!kbhit());
 closegraph();
 return 0;
}

Draw a polygon. Use Bresenham line drawing algorithm

//Reference:- David F. Rogers "Procedural Elements for 
   Computer Graphics" Second Edition
#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 coordinates ",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 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 i,dx=xb-xa,dy=yb-ya,p,twod,twodd,s;
 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;
  } 
  p=2*dx-dy;
  if(dy<0)
  {
    s=-1;
       dy=-dy;
  }
  else
     s=1;
  putpixel(xa,ya,10);
  for(i=0;i<dy;i++)
  {
   while(p>0)
   {
    xa++;
    p-=2*dy;
    putpixel(xa,ya,10);
   }
   ya+=s;
   p+=2*dx;
   putpixel(xa,ya,10);
  }
 }
 else
 {
  if(ya>yb)
  {
   Bresenham(xb,yb,xa,ya);
   return;
  }
  p=2*dy-dx;
  if(dx<0)
   {
       s=-1;
       dx=-dx;
    }
  else
     s=1;
  putpixel(xa,ya,10);
  for(i=0;i<dx;i++)
  {
   while(p>0)
   {
    ya++;
    p-=2*dx;
    putpixel(xa,ya,10);
   }
   xa+=s;
   p+=2*dy;
   putpixel(xa,ya,10);
  } 
 }
}

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

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

2D Transformation

Write a program to apply combination of transformation, rotation, reflection and shearing) on the following objects. a)Polygon b)Circle #include<stdio.h> #include<graphics.h> int midx,midy; int polx[10],poly[10],n,pox,poy,r,rnew,polxnew[10],polynew[10],poxnew,poynew; float trans[3][3],result[3][1],xandy[3][1]; void screendiv() { int winx,winy; cleardevice(); winx=getmaxx(); winy=getmaxy(); line(winx/2,0,winx/2,winy); line(0,winy/2,winx,winy/2); midx=winx/2; midy=winy/2; } void readpoly() { int i; for(i=0;i<n;i++) { printf("\nEnter the %d co=ordinates ",i+1); scanf("%d%d",&polx[i],&poly[i]); } } void readcircle() { printf("Enter the orgin "); scanf("%d%d",&pox,&poy); printf("Enter the radius "); scanf("%d",&r); } void drawcircle(int *x,int *y,int *rad) { circle(midx+*x,midy-*y,*rad); } void drawpolygon(int *x,int *y) { int i; for(i=0;i<n;i++) line(midx+*(x+i),midy-*(y+i),midx+*(x+((i+1)%n)),midy-*(y+((i+1)%n))); } 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 transformcirc() { xandy[0][0]=pox; xandy[1][0]=poy; xandy[2][0]=1.0; mulmat(); poxnew=(int)result[0][0]; poynew=(int)result[1][0]; } void transformpoly() { int i; for(i=0;i<n;i++) { xandy[0][0]=polx[i]; xandy[1][0]=poly[i]; xandy[2][0]=1; mulmat(); polxnew[i]=result[0][0]; polynew[i]=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,i,j,ch,dx,dy,rx,ry,chr; float sx,sy,ang; while(1) { makeIdentity(); printf("\nMENU\n-------\n1.CIRCLE\n2.POLYGON\n3.EXIT"); printf("\nEnter the choice "); scanf("%d",&choice); switch(choice) { case 1: readcircle(); initgraph(&gd,&gm,NULL); screendiv(); drawcircle(&pox,&poy,&r); break; case 2: printf("\nEnter the no of sides "); scanf("%d",&n); readpoly(); initgraph(&gd,&gm,NULL); screendiv(); drawpolygon(&polx[0],&poly[0]); break; case 3: exit(0); default:printf("\nInvalid choice "); } while (!kbhit()); closegraph(); printf("\nMENU\n-------\n1.TRANSLATE\n2.ROTATION\n3.SCALING\n4.REFLECTION\n5.SHEAR"); printf("\nEnter the choice "); scanf("%d",&ch); switch(ch) { case 1: printf("\ndx and dy"); scanf("%d%d",&dx,&dy); trans[0][2]=dx; trans[1][2]=dy; initgraph(&gd,&gm,NULL); if(choice==1) { transformcirc(); screendiv(); drawcircle(&pox,&poy,&r); drawcircle(&poxnew,&poynew,&r); } else { transformpoly(); screendiv(); drawpolygon(&polx[0],&poly[0]); drawpolygon(&polxnew[0],&polynew[0]); } break; case 2: printf("\nEnter the reference points "); scanf("%d%d",&rx,&ry); printf("\nEnter the angle "); scanf("%f",&ang); ang=(float)(ang*(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); initgraph(&gd,&gm,NULL); if(choice==1) { transformcirc(); screendiv(); drawcircle(&pox,&poy,&r); drawcircle(&poxnew,&poynew,&r); } else { transformpoly(); screendiv(); drawpolygon(&polx[0],&poly[0]); drawpolygon(&polxnew[0],&polynew[0]); } break; case 3: printf("\nEnter the reference points "); scanf("%d%d",&rx,&ry); if(choice==1) { printf("\nEnter the scaling factor "); scanf("%f",&sx); rnew=(float)r*sx; pox=rx; poy=ry; initgraph(&gd,&gm,NULL); screendiv(); drawcircle(&pox,&poy,&r); drawcircle(&poxnew,&poynew,&rnew); } else { printf("\nEnter the scaling factor,sx and sy"); scanf("%f%f",&sx,&sy); trans[0][0]=sx; trans[0][2]=(1-sx)*rx; trans[1][1]=sy; trans[1][2]=(1-sy)*ry; transformpoly(); initgraph(&gd,&gm,NULL); screendiv(); drawpolygon(&polx[0],&poly[0]); drawpolygon(&polxnew[0],&polynew[0]); } break; case 4: printf("\nMENU FOR REFLECTION\n-------\n1.X AXIS\n2.Y AXIS\n3.ORGIN\n4.Y=X\n5.Y=-X"); printf("\nEnter the choice "); scanf("%d",&chr); switch(chr) { case 1: trans[1][1]=-1; break; case 2: trans[0][0]=-1; break; case 3: trans[0][0]=-1; trans[1][1]=-1; break; case 4: trans[0][0]=0; trans[0][1]=1; trans[1][0]=1; trans[1][1]=0; break; case 5: trans[0][0]=0; trans[0][1]=-1; trans[1][0]=-1; trans[1][1]=0; break; default:printf("\nInvlalid choice"); break; } initgraph(&gd,&gm,NULL); if(choice==1) { transformcirc(); screendiv(); drawcircle(&pox,&poy,&r); drawcircle(&poxnew,&poynew,&r); } else { transformpoly(); screendiv(); drawpolygon(&polx[0],&poly[0]); drawpolygon(&polxnew[0],&polynew[0]); } break; case 5: if(choice==1) printf("\nShear not possible "); else { printf("\nMENU FOR SHEAR\n-------\n1.X DIRECTION SHEAR RELATIVE TO X AXIS\n2.X DIRECTION SHEAR RELATIVE TO REFERENCE\n3.Y DIRECTION SHEAR RELATIVE TO REFERENCE"); printf("\nEnter the choice "); scanf("%d",&chr); printf("\nEnter the shear factor "); scanf("%f",&sx); switch(chr) { case 1: trans[0][1]=sx; break; case 2: printf("\nEnter y reference point "); scanf("%d",&ry); trans[0][1]=sx; trans[0][2]=-1*sx*ry; break; case 3: printf("\nEnter x reference point "); scanf("%d",&rx); trans[1][0]=sx; trans[1][2]=-1*sx*rx; break; default:printf("\nInvlalid choice"); } initgraph(&gd,&gm,NULL); transformpoly(); screendiv(); drawpolygon(&polx[0],&poly[0]); drawpolygon(&polxnew[0],&polynew[0]); } break; default:printf("\nInvlalid choice"); } while (!kbhit()); closegraph(); } return 0; }

Sutherland-Hodgeman Polygon Clipping

#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 drawpolygon(); void showpoints(); void clippolygon() { float m=0.0; np=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; np++; j++; } else if(p2.x==wmin.x) { pout[j]=p2; np++; 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; np=np+2; j=j+2; } else { p2=ipt; pout[j]=ipt; np++; j++; } } } newpin(); outtextxy(220,10,"After left clip"); drawpolygon(); while (!kbhit()); closegraph(); initgraph(&gd,&gm,"NULL"); np=0; 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; np++; j++; } else if(p2.x==wmax.x) { pout[j]=p2; np++; 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; np=np+2; j=j+2; } else { p2=ipt; pout[j]=ipt; np++; j++; } } } newpin(); outtextxy(220,10,"After right clip"); drawpolygon(); while (!kbhit()); closegraph(); initgraph(&gd,&gm,"NULL"); np=0; 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; np++; j++; } else if(p2.y==wmin.y) { pout[j]=p2; np++; 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; np=np+2; j=j+2; } else { p2=ipt; pout[j]=ipt; np++; j++; } } } newpin(); outtextxy(220,10,"After bottom clip"); drawpolygon(); while (!kbhit()); closegraph(); initgraph(&gd,&gm,"NULL"); np=0; 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; np++; j++; } else if(p2.y==wmax.y) { pout[j]=p2; np++; 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; np=np+2; j=j+2; } else { p2=ipt; pout[j]=ipt; np++; j++; } } } newpin(); outtextxy(220,10,"clipped Polygon"); showpoints(); drawpolygon(); } void newpin() { int i,j; for(j=0,i=0;j<np;j++,i++) pin[i]=pout[j]; n=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); } int main() { int c,j; 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); } initgraph(&gd,&gm,"NULL"); outtextxy(220,10,"Orginal Polygon"); drawpolygon(); while (!kbhit()); closegraph(); initgraph(&gd,&gm,"NULL"); clippolygon(); while (!kbhit()); closegraph(); return 0; }

cohen-sutherland Line clipping

#include<stdio.h> #include<graphics.h> #define left 1 #define right 2 #define top 8 #define bottom 4 int xa,ya,xb,yb; int wxmax,wxmin,wymax,wymin; int regioncode(int x,int y) { int ret=0; if(x<wxmin) ret|=left; if(x>wxmax) ret|=right; if(y<wymin) ret|=bottom; if(y>wymax) ret|=top; return ret; } int clipline() { int rcode1=regioncode(xa,ya); int rcode2=regioncode(xb,yb); int x=0,y=0; float m=0.0; if(xb!=xa) m=(float)(yb-ya)/(float)(xb-xa); while(1) { if(rcode1==0&&rcode2==0) return 1; else if((rcode1&rcode2)!=0) return 0; else { int outcode=(rcode1==0)?rcode2:rcode1; if(outcode&top) { if(m!=0) x=xa+(float)(1/m)*(float)(wymax-ya); else x=xa; y=wymax; } else if(outcode&bottom) { if(m!=0) x=xa+(float)(1/m)*(float)(wymin-ya); else x=xa; y=wymin; } else if(outcode&left) { x=wxmin; y=ya+(float)(m)*(float)(wxmin-xa); } else if(outcode&right) { x=wxmax; y=ya+(float)(m)*(float)(wxmax-xa); } if(rcode1==outcode) { xa=x; ya=y; rcode1=regioncode(xa,ya); } else { xb=x; yb=y; rcode2=regioncode(xb,yb); } } } } int main() { int gd=DETECT,gm=VGAMAX; printf("\nEnter the diagonal of rectangle, wxmin,wymin,wxmax,wymax"); scanf("%d%d%d%d",&wxmin,&wymin,&wxmax,&wymax); printf("\nxmin ymin (%d,%d)\n xmax ymax (%d,%d)",wxmin,wymin,wxmax,wymax); printf("\nEnter the cordinates of line "); scanf("%d%d%d%d",&xa,&ya,&xb,&yb); initgraph(&gd,&gm,NULL); printf("Before Clipping "); rectangle(wxmin,wymin,wxmax,wymax); line(xa,ya,xb,yb); while(!kbhit()); closegraph(); initgraph(&gd,&gm,NULL); printf("After Clipping "); rectangle(wxmin,wymin,wxmax,wymax); if(clipline()) line(xa,ya,xb,yb); while(!kbhit()); closegraph(); return 0; }

Filling - recursive and non-recursive method

3. Write a program to implement following seed filling algorithms on polygon, circle and an ellipse.
    a) floodfill   b) boundary fill

Recursive and Non Recursive method
#include<stdio.h> #include<graphics.h> #define ROUND(a)((int)(a+0.5)) struct stack { int x; int y; struct stack *next; }*p,*temp,*head=NULL; typedef struct stack t1; void push(int,int); void pop(); int gd=DETECT,gm=VGAMAX; int xf,yf; void plotcircle(int,int,int,int); void plotellipse(int,int,int,int); void floodfill4(int,int,int,int); void floodfill4nonrec(int,int,int,int); void boundaryfill4(int,int,int,int); void boundaryfill4nonrec(int,int,int,int); void midpointcircle(int,int,int); void midpointellipse(int,int,int,int); void lineDDA(int,int,int,int); int main() { int xc,yc,rx,ry,r,n,choice,ch,newcolor; int oldcolor=getbkcolor(); printf("\nEnter the fill color "); scanf("%d",&newcolor); while(1) { printf("\nMENU\n-------\n1.FLOOD FILL\n2.BOUNDARY FILL\n3.FLOOD FILL NONRECURSIVE\n4.BOUNDARY FILL NONRECURSIVE\n5.EXIT"); printf("\nEnter the choice "); scanf("%d",&ch); if(ch>4) exit(0); printf("\nMENU\n-------\n1.CIRCLE\n2.ELLIPSE\n3.POLYGON\n4.EXIT"); printf("\nEnter the choice "); scanf("%d",&choice); switch(choice) { case 1: printf("\nEnter the center "); scanf("%d%d",&xc,&yc); printf("\nEnter the radius "); scanf("%d",&r); printf("\nEnter the inner point "); scanf("%d%d",&xf,&yf); initgraph(&gd,&gm,NULL); midpoint(xc,yc,r); break; case 2: printf("\nEnter the center "); scanf("%d%d",&xc,&yc); printf("\nEnter the minor and major axes "); scanf("%d%d",&rx,&ry); printf("\nEnter the inner point "); scanf("%d%d",&xf,&yf); initgraph(&gd,&gm,NULL); midpointellipse(xc,yc,rx,ry); break; case 3: printf("\nEnter the no of sides "); scanf("%d",&n); int i,ax[10],ay[10]; for(i=0;i<n;i++) { printf("\nEnter the %d co=ordinates ",i+1); scanf("%d%d",&ax[i],&ay[i]); } printf("\nEnter the inner point "); scanf("%d%d",&xf,&yf); initgraph(&gd,&gm,NULL); for(i=0;i<n;i++) lineDDA(ax[i],ay[i],ax[(i+1)%n],ay[(i+1)%n]); break; case 4: exit(0); break; default:printf("\nInvalid Choice"); } if(ch==1) floodfill4(xf,yf,newcolor,oldcolor); else if(ch==2) boundaryfill4(xf,yf,newcolor,10); else if(ch==2) floodfill4nonrec(xf,yf,newcolor,oldcolor); else boundaryfill4nonrec(xf,yf,newcolor,10); while (!kbhit()); closegraph(); } return 0; } 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 floodfill4(int x,int y,int new,int old) { if(getpixel(x,y)==old) { putpixel(x,y,new); floodfill4(x+1,y,new,old); floodfill4(x-1,y,new,old); floodfill4(x,y+1,new,old); floodfill4(x,y-1,new,old); } } 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); } } void floodfill4nonrec(int x,int y,int new,int old) { if(getpixel(x,y)==new) return; push(x,y); putpixel(x,y,new); while(head!=NULL) { pop(); x=p->x; y=p->y; if(getpixel(x+1,y)==old) { putpixel(x+1,y,new); push(x+1,y); } if(getpixel(x-1,y)==old) { putpixel(x-1,y,new); push(x-1,y); } if(getpixel(x,y+1)==old) { putpixel(x,y+1,new); push(x,y+1); } if(getpixel(x,y-1)==old) { putpixel(x,y-1,new); push(x,y-1); } } } void boundaryfill4nonrec(int x,int y,int new,int bdry) { if(getpixel(x,y)==bdry) return; push(x,y); putpixel(x,y,new); while(head!=NULL) { pop(); x=p->x; y=p->y; if((getpixel(x+1,y)!=bdry)&&(getpixel(x+1,y)!=new)) { putpixel(x+1,y,new); push(x+1,y); } if((getpixel(x-1,y)!=bdry)&&(getpixel(x-1,y)!=new)) { putpixel(x-1,y,new); push(x-1,y); } if((getpixel(x,y+1)!=bdry)&&(getpixel(x,y+1)!=new)) { putpixel(x,y+1,new); push(x,y+1); } if((getpixel(x,y-1)!=bdry)&&(getpixel(x,y-1)!=new)) { putpixel(x,y-1,new); push(x,y-1); } } } 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,10); putpixel(xc-x,yc+y,10); putpixel(xc+x,yc-y,10); putpixel(xc-x,yc-y,10); } void lineDDA(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),10); for(k=0;k<steps;k++) { x+=xinc; y+=yinc; putpixel(ROUND(x),ROUND(y),10); } } void push(int x1,int y1) { temp=(t1 *)malloc(sizeof(t1)); temp->x=x1; temp->y=y1; if(head==NULL) { head=temp; temp->next=NULL; } else { temp->next=head; head=temp; } } void pop() { p=head; head=head->next; }

Circle Drawing

2. Write a program to draw a circle using parametric and midpoint algorithm.
#include<stdio.h> #include<graphics.h> #include<math.h> #define ROUND(a)((int)(a+0.5)) void parametric(int,int,int); void midpoint(int,int,int); void plot(int,int,int,int); int main() { int xc,yc,r,choice,gd=DETECT,gm=VGAMAX; printf("\nEnter the center "); scanf("%d%d",&xc,&yc); printf("\nEnter the radius "); scanf("%d",&r); while(1) { printf("\nMENU\n-------\n1.PARAMETRIC\n2.MIDPOINT\n3.EXIT"); printf("\nEnter the choice "); scanf("%d",&choice); initgraph(&gd,&gm,NULL); switch(choice) { case 1: parametric(xc,yc,r); break; case 2: midpoint(xc,yc,r); break; case 3: exit(0); break; default:printf("\nInvalid Choice"); } while (!kbhit()); closegraph(); } return 0; } void parametric(int xc,int yc,int r) { int i,x,y; for(i=0;i<=90;i++) { x=r*cos(3.1412*i/180); y=r*sin(3.1412*i/180); plot(xc,yc,x,y); } } 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); }

Line Drawing

1. Write a program to draw line using DDA algorithm and Bresenham's algorithm.
#include<stdio.h> #include<graphics.h> #include<math.h> #define ROUND(a)((int)(a+0.5)) void DDA(int,int,int,int); void Bresenham(int,int,int,int); int main() { int xa,ya,xb,yb,choice,gd=DETECT,gm=VGAMAX; printf("\nEnter first cordinates "); scanf("%d%d",&xa,&ya); printf("\nEnter second cordinates "); scanf("%d%d",&xb,&yb); while(1) { printf("\nMENU\n-------\n1.DDA\n2.BRESENTHAM\n3.EXIT"); printf("\nEnter the choice "); scanf("%d",&choice); initgraph(&gd,&gm,NULL); switch(choice) { case 1: DDA(xa,ya,xb,yb); break; case 2: Bresenham(xa,ya,xb,yb); break; case 3: exit(0); break; default:printf("\nInvalid Choice"); } 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),10); for(k=0;k<steps;k++) { x+=xinc; y+=yinc; putpixel(ROUND(x),ROUND(y),10); } } void Bresenham(int xa,int ya,int xb,int yb) { int dx=abs(xa-xb),dy=abs(ya-yb),p,twody,twodydx,x,y,xEnd; p=2*dy-dx; twody=2*dy; twodydx=2*(dy-dx); if(xa>xb) { x=xb; y=yb; xEnd=xa; } else { x=xa; y=ya; xEnd=xb; } putpixel(x,y,10); while(x<xEnd) { x++; if(p<0) p+=twody; else { y++; p+=twodydx; } putpixel(x,y,10); } }