Mandelbrot

mandelbrot.c

#include <stdio.h>

#define R 50
#define C 150
#define K 50

int main()
{
    int r, c;
    for (r=0; r<R; r++)
    {
        for (c=0; c<C; c++)
        {
            double x = 0;
            double y = 0;

            double x0 = -2.5 + 4.0 * ((double)c)/C;
            double y0 = -1.5 + 3.0 * ((double)r)/R;

            int k;
            for (k=0; k<K; k++)
            {
                if (x*x + y*y > 2.0)
                {
                    printf(" ");
                    break;
                }

                double xnew = x*x - y*y + x0;
                double ynew = 2*x*y + y0;
                x = xnew;
                y = ynew;

            }
            if (k == K)
                printf("*");
        }
        printf("\n");
    }

    return 0;
}