用vc++实现IFS分形算法画一棵树

void IFS()
{
	double a[][8]={
		0,0.195,-0.488 ,0.344 ,0.433 ,0.4431 ,0.2452 ,0.25 ,\
		0,0.462 ,0.414 ,-0.252 ,0.361 ,0.2511 ,0.5692 ,0.25,\
		0,-0.058 ,-0.07 ,0.453 ,-0.111 ,0.5976 ,0.0969 ,0.25 ,\
		0,-0.035 ,0.07 ,-0.469 ,-0.022 ,0.4884 ,0.5069 ,0.2,\
		0,-0.637 ,0 ,0 ,0.501 ,0.8562 ,0.2513 ,0.05};
	double aa,bb;
	aa = a[0][1];
	double x = 900;

	srand( (unsigned)time( NULL ) );

	double x0,y0;
	x0=0;
	y0=0;
	initgraph(1000,1000);
	//setorigin(getwidth() / 2, getheight() / 2);         // 设置坐标原点为窗口中心
	for (int k=0;k<10000;k++){
		int r = rand()%100;
		double x1,y1;

		if(r<=25){
			x1 = a[0][1]*x0+a[0][2]*y0+a[0][5];
			y1 = a[0][3]*x0+a[0][4]*y0+a[0][6];
		}else if(r<=50){
			x1 = a[1][1]*x0+a[1][2]*y0+a[1][5];
			y1 = a[1][3]*x0+a[1][4]*y0+a[1][6];
		}else if (r<=75){
			x1 = a[2][1]*x0+a[2][2]*y0+a[2][5];
			y1 = a[2][3]*x0+a[2][4]*y0+a[2][6];
		}else if (r<=95){
			x1 = a[3][1]*x0+a[3][2]*y0+a[3][5];
			y1 = a[3][3]*x0+a[3][4]*y0+a[3][6];
		}else{
			x1 = a[4][1]*x0+a[4][2]*y0+a[4][5];
			y1 = a[4][3]*x0+a[4][4]*y0+a[4][6];
		}
		x0 = x1;
		y0 = y1;
		putpixel(x1*x,1000-y1*x,HSLtoRGB(200,r*6,200));


		}
	getch();
	closegraph();
}

用到了 www.easyx.cn