精彩作图

原文： 有没有一段代码，让你觉得人类的智慧也可以璀璨无比？

原文是C++代码，且产生的图形是PPM格式，下面程序是C语言版本的，输出为BMP格式的图片。

#include <stdio.h>
#include <math.h> 
#include <stdlib.h>
#define DIM 1024
#define DM1 (DIM-1)
#define _sq(x) ((x)*(x)) // square
#define _cb(x) abs((x)*(x)*(x)) // absolute value of cube
#define _cr(x) (unsigned char)(pow((x),1.0/3.0)) // cube root
 
unsigned char GR(int,int);
unsigned char BL(int,int);


unsigned char RD(int i,int j){
// YOUR CODE HERE
}
unsigned char GR(int i,int j){
// YOUR CODE HERE
}
unsigned char BL(int i,int j){
// YOUR CODE HERE
}




FILE *fp;
void pixel_write(int i, int j)
{
static unsigned char color[3];
color[0] = RD(i,j)&255;
color[1] = GR(i,j)&255;
color[2] = BL(i,j)&255;
fwrite(color, 1, 3, fp);
}


int main()
{
char bmphead[] = {0X42 ,0X4D ,0X36 ,0X00 ,0X30 ,0X00 ,0X00 ,0X00 ,0X00 ,0X00 ,0X36 ,0X00 ,0X00 ,0X00 ,
             0X28 ,0X00 ,0X00 ,0X00 ,0X00 ,0X04 ,0X00 ,0X00 ,0X00 ,0X04 ,0X00 ,0X00 ,0X01 ,0X00 ,
 0X18 ,0X00 ,0X00 ,0X00 ,0X00 ,0X00 ,0X00 ,0X00 ,0X30 ,0X00 ,0X00 ,0X00 ,0X00 ,0X00 ,
 0X00 ,0X00 ,0X00 ,0X00 ,0X00 ,0X00 ,0X00 ,0X00 ,0X00 ,0X00 ,0X00 ,0X00} ;

fp = fopen("MathPic.bmp","wb");
fwrite(bmphead, 1, 54, fp);

for(int j=0;j<DIM;j++)
for(int i=0;i<DIM;i++)
pixel_write(i,j);

fclose(fp);
return 0;
}

更多精彩参见 出处

示例

unsigned char RD(int i,int j){
   double a=0,b=0,c,d,n=0;
   while((c=a*a)+(d=b*b)<4&&n++<880)
   {b=2*a*b+j*8e-9-.645411;a=c-d+i*8e-9+.356888;}
   return 255*pow((n-80)/800,3.);
}
unsigned char GR(int i,int j){
   double a=0,b=0,c,d,n=0;
   while((c=a*a)+(d=b*b)<4&&n++<880)
   {b=2*a*b+j*8e-9-.645411;a=c-d+i*8e-9+.356888;}
   return 255*pow((n-80)/800,.7);
}
unsigned char BL(int i,int j){
   double a=0,b=0,c,d,n=0;
   while((c=a*a)+(d=b*b)<4&&n++<880)
   {b=2*a*b+j*8e-9-.645411;a=c-d+i*8e-9+.356888;}
   return 255*pow((n-80)/800,.5);
}

----------------------------------------------------------------------------------------------------------------------------

unsigned char RD(int i,int j){
    float s=3./(j+99);
    float y=(j+sin((i*i+_sq(j-700)*5)/100./DIM)*35)*s;
    return (int((i+DIM)*s+y)%2+int((DIM*2-i)*s+y)%2)*127;
}
unsigned char GR(int i,int j){
    float s=3./(j+99);
    float y=(j+sin((i*i+_sq(j-700)*5)/100./DIM)*35)*s;
    return (int(5*((i+DIM)*s+y))%2+int(5*((DIM*2-i)*s+y))%2)*127;
}
unsigned char BL(int i,int j){
    float s=3./(j+99);
    float y=(j+sin((i*i+_sq(j-700)*5)/100./DIM)*35)*s;
    return (int(29*((i+DIM)*s+y))%2+int(29*((DIM*2-i)*s+y))%2)*127;
}