ppm1024.cpp

#include <iostream>
#include <cmath>
#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 short)(pow((x),1.0/3.0))  // cube root

unsigned short red_fn(int i,int j){
    // YOUR CODE HERE
}
unsigned short green_fn(int i,int j){
    // YOUR CODE HERE
}
unsigned short blue_fn(int i,int j){
    // YOUR CODE HERE
}

void pixel_write(int,int);
FILE *fp;
int main(){
    fp = fopen("MathPic","wb");
    fprintf(fp, "P6\n%d %d\n1023\n", DIM, DIM);
    for(int j=0;j<DIM;j++)
        for(int i=0;i<DIM;i++)
            pixel_write(i,j);
    fclose(fp);
    return 0;
}
void pixel_write(int i, int j){
    static unsigned short color[3];
    color[0] =   red_fn(i,j)&DM1;
    color[1] = green_fn(i,j)&DM1;
    color[2] =  blue_fn(i,j)&DM1;
    fwrite(color, 2, 3, fp);
}