sier_gasket.c

包含了多个实现Sierpinski分形（三角形、正方形、六边形）的Matlab程序和C程序。

#include "graphics.h"
void Sier_Gasket(n,x1,y1,x2,y2,x3,y3)
int n;
int x1,y1,x2,y2,x3,y3;
{
 int tri[8];
 if(n==0)return;
 tri[0]=(x1+x2)/2;
 tri[1]=(y1+y2)/2;
 tri[2]=(x2+x3)/2;
 tri[3]=(y2+y3)/2;
 tri[4]=(x1+x3)/2;
 tri[5]=(y1+y3)/2;
 tri[6]=tri[0];
 tri[7]=tri[1];
 moveto(tri[0],tri[1]);
 lineto(tri[2],tri[3]);
 lineto(tri[4],tri[5]);
 lineto(tri[6],tri[7]);
 floodfill((x1+x2+x3)/3,(y1+y2+y3)/3,RED);
 Sier_Gasket(n-1,x1,y1,(x1+x2)/2,(y1+y2)/2,(x1+x3)/2,(y1+y3)/2);
 Sier_Gasket(n-1,(x1+x2)/2,(y1+y2)/2,x2,y2,(x2+x3)/2,(y2+y3)/2);
 Sier_Gasket(n-1,(x1+x3)/2,(y1+y3)/2,(x2+x3)/2,(y2+y3)/2,x3,y3);
}

main()
{
int gdriver=DETECT,gmode;
initgraph(&gdriver,&gmode,"");
cleardevice();
setfillstyle(1,2);
setcolor(RED);
moveto(10,420);
lineto(10,420);
lineto(320,20);
lineto(640,420);
lineto(10,420);
setfillstyle(1,3);
Sier_Gasket(1,10,420,320,20,640,420);
getch();
Sier_Gasket(2,10,420,320,20,640,420);
getch();
Sier_Gasket(6,10,420,320,20,640,420);
getch();
closegraph();
}