第C語言圣誕樹的實(shí)現(xiàn)示例效果如下：

#define_CRT_SECURE_NO_WARNINGS1

#includemath.h

#includestdio.h

#includestdlib.h

#includestring.h

#definePI3

floatsx,sy;

typedeffloatMat[4][4];

typedeffloatVec[4];

voidscale(Mat*m,floats){

Mattemp={{s,0,0,0},{0,s,0,0},{0,0,s,0},{0,0,0,1}};

memcpy(m,temp,sizeof(Mat));

voidrotateY(Mat*m,floatt){

floatc=cosf(t),s=sinf(t);

Mattemp={{c,0,s,0},{0,1,0,0},{-s,0,c,0},{0,0,0,1}};

memcpy(m,temp,sizeof(Mat));

voidrotateZ(Mat*m,floatt){

floatc=cosf(t),s=sinf(t);

Mattemp={{c,-s,0,0},{s,c,0,0},{0,0,1,0},{0,0,0,1}};

memcpy(m,temp,sizeof(Mat));

voidtranslate(Mat*m,floatx,floaty,floatz){

Mattemp={{1,0,0,x},{0,1,0,y},{0,0,1,z},{0,0,0,1}};

memcpy(m,temp,sizeof(Mat));

voidmul(Mat*m,Mata,Matb){

Mattemp;

for(intj=0;jj++)

for(inti=0;ii++){

temp[j][i]=0.0f;

for(intk=0;kk++)

temp[j][i]+=a[j][k]*b[k][i];

memcpy(m,temp,sizeof(Mat));

voidtransformPosition(Vec*r,Matm,Vecv){

Vectemp={0,0,0,0};

for(intj=0;jj++)

for(inti=0;ii++)

temp[j]+=m[j][i]*v[i];

memcpy(r,temp,sizeof(Vec));

floattransformLength(Matm,floatr){

returnsqrtf(m[0][0]*m[0][0]+m[0][1]*m[0][1]+m[0][2]*m[0][2])*r;

floatsphere(Vecc,floatr){

floatdx=c[0]-sx,dy=c[1]-sy;

floata=dx*dx+dy*dy;

returnar*rsqrtf(r*r-a)+c[2]:-1.0f;

floatopUnion(floatz1,floatz2){

returnz1z2z1:z2;

floatf(Matm,intn){

floatz=-1.0f;

for(floatr=0.0f;r0.8f;r+=0.02f){

Vecv={0.0f,r,0.0f,1.0f};

transformPosition(v,m,v);

z=opUnion(z,sphere(v,transformLength(m,0.05f*(0.95f-r))));

if(n0){

Matry,rz,s,t,m2,m3;

rotateZ(rz,1.8f);

for(intp=0;pp++){

rotateY(ry,p*(2*PI/6));

mul(m2,ry,rz);

floatss=0.45f;

for(floatr=0.2f;r0.8f;r+=0.1f){

scale(s,ss);

translate(t,0.0f,r,0.0f);

mul(m3,s,m2);

mul(m3,t,m3);

mul(m3,m,m3);

z=opUnion(z,f(m3,n-1));

ss*=0.8f;

returnz;

floatf0(floatx,floaty,intn){

sx=x;

sy=y;

Matm;

scale(m,1.0f);

returnf(m,n);

intmain(intargc,char*argv[]){

//這個(gè)里可以加打印信息喲，為他添加專屬的printf

intn=argc1atoi(argv[1]):3;

floatzoom=argc2atof(argv[2]):1.0f;

for(floaty=0.8f;y-0.0f;y-=0.02f/zoom,putchar('\n'))

for(floatx=-0.35f;x0.35f;x+=0.01f/zoom){

floatz=f0(x,y,n);

if(z-1.0f){

floatnz=0.001f;

floatnx=f0(x+nz,y,n)-z;

floatny=f0(x,y+nz,n)-z;

floatnd=sqrtf(nx*nx+ny*ny+nz*nz);

floatd=(nx-