Pasted

#include <stdio.h>
 
#define ANCHURA 512
#define ALTURA  512
 
void guardaPGMd(char* nombre, int anchura, int altura, double *pixels,
                double pixel_min, double pixel_max)
                {
                    int i, j, ij, p;
                    FILE* imagen;
 
                    imagen=fopen(nombre,"wb");
                    fprintf(imagen,"P2");
                    fprintf(imagen,"#guardaPGMd %s\n", nombre);
                    fprintf(imagen,"%d %d\n",anchura,altura);
                    fprintf(imagen,"255\n");
 
                    ij=0;
                    for(i=0;i<altura; i++)
                    {
                        for(j=0;j<anchura;j++)
                            {
                                p=(int)(255*(pixels[ij]-pixel_min))/(pixel_max-pixel_min);
                                if(p<0) p=0;
                                if(p>255) p=255;
 
                                fprintf(imagen," %d",p);
 
                                ij++;
                            }
                            fprintf(imagen,"\n");
                    }
                    fclose(imagen);
 
                    return;
                }
int main(int argc, char **argv)
{
    int i, j;
    double x,y,z;
    double f[ALTURA][ANCHURA];
 
    for(i=0;i<ALTURA;i++)
    {
        for(j=0;j<ANCHURA;j++)
        {
            x=2*(double)j/(ANCHURA-1)-1.0;
            y=2*(double)i/(ALTURA-1)-1.0;
            z=1.0/(1+(2*x*x+3*y*y));
 
            f[i][j]=z;
        }
 
    }
 
    guardaPGMd("prueba.pgm",ANCHURA,ALTURA,(double*)f,0,1.0);
 
    return 0;
}