分形--细胞自动装置2

#include "glib.h"
int A=1,B=0,C=1,MAX_COL=4,Y=256;
void main()
{
  int x,y;
  static int f[640],g[640];
  ginit();
  set0(0,72);
  for (x=0;x<640;x++)
  {
 if (320==x)
 {
   f[x]=1;
   putpixel(x,0,f[x]);
 }
 else
   f[x]=0;
  }

  for (y=1;y<Y;y++)
  {
 for (x=0;x<640;x++)
 {
   if (!x || 639==x)
  if (!x)
    g[x]=B*f[x]+C*f[x+1];
  else
    g[x]=A*f[x-1]+B*f[x];
   else
  g[x]=A*f[x-1]+B*f[x]+C*f[x+1];
   if (g[x])
   {
  g[x]=(g[x]-1)%MAX_COL+1;
  putpixel(320+(x-320)/sqrt(3.0),y,g[x]);
   }
 }
 for (x=0;x<640;x++)
   f[x]=g[x];
  }
  gend();
}