/* PROGRAM BY KIRAN */


	#include<stdio.h>
	#include<conio.h>
	#include<math.h>
	#include<graphics.h>
	#define pi 3.14

	float xa,ya,xb,yb,xc,yc;
	float xr,yr;
	int i,j,k,z,t;
	float a[3][3],c[3][3],o,r;
	float sx,sy,sum;
	int s,e;
	float identity[3][3];
	struct GOOD
		    {
		     float x;
		     float y;
		    }pts[50];
	main()
      {
	int gm=DETECT,gd=DETECT;
	initgraph(&gm,&gd,"c:\\turboc\\bgi");
	clrscr();
	cleardevice();
	printf("\nEnter the points for rotation : ");
	scanf("%f%f",&xr,&yr);
	printf("\nEnter the angle for rotation: ");
	scanf("%f",&r);
	printf("\nEnter the sx and sy for scaling :");
	scanf("%f%f",&sx,&sy);
	printf("\nEnter the sides of figure :");
	scanf("%d",&s);


	    for(e=1;e<=s;e++)
	   {
	   printf("\nEnter the point of sides one by one :");
	   scanf("%f%f",&pts[e].x,&pts[e].y);
	   }
	   clrscr();
	   cleardevice();

	   for(e=1;e<s;e++)
	  {

	   line(pts[e].x,pts[e].y,pts[e+1].x,pts[e+1].y);
	   }

	   line(pts[e].x,pts[e].y,pts[1].x,pts[1].y);

	getche();

	ident();
       for(e=1;e<=s;e++)
       {
	identity[0][0]=pts[e].x;
	identity[1][0]=pts[e].y;
	identity[2][0]=1;
	call();
	}

	   for(e=1;e<s;e++)
	  {

	   line(pts[e].x,pts[e].y,pts[e+1].x,pts[e+1].y);
	   }

	   line(pts[e].x,pts[e].y,pts[1].x,pts[1].y);
	  getche();
	  closegraph();
	 }



	call()
       {
	o=(r*pi)/180;
	trs();
	for(i=0;i<3;i++)
	{
	 for(j=0;j<3;j++)
	 {
	  sum=0;
	  for(k=0;k<3;k++)
	  {
	   sum=sum+a[i][k]*identity[k][j];
	   c[i][j]=sum;
	  }
	 }
	}

	 pts[e].x=c[0][0];
	 pts[e].y=c[1][0];

	return;
       }

	ident()
       {
	identity[0][0]=1;
	identity[0][1]=0;
	identity[0][2]=0;
	identity[1][0]=0;
	identity[1][1]=1;
	identity[1][2]=0;
	identity[2][0]=0;
	identity[2][1]=0;
	identity[2][2]=1;
	return;
       }

	trs()
       {
	a[0][0]=sx*cos(o);
	a[0][1]=-sy*sin(o);
	a[0][2]=-xr*sx*cos(o)+yr*sy*sin(o)+xr;
	a[1][0]=sx*sin(o);
	a[1][1]=sy*cos(o);
	a[1][2]=-xr*sx*sin(o)-yr*sy*cos(o)+yr;
	a[2][0]=0;
	a[2][1]=0;
	a[2][2]=1;
	return;
       }