Fortran - Conjunto de Mandelbrot

Pido ayuda a algun programador que no le importase echarme una manilla con un programa que se me esta atravesando un poco. Lo que me gustaría conseguir es un programa en fortran que sea capaz de calcularme el conjunto de Mandelbrot en el plano complejo considerando que pertenece al conjunto si se llega a un recuento de 100. Hay que empezar por el punto (-2,-1) hasta el (1,1) y llevar la información a un fichero para luego de ahí pasarlo a gnuplot. Les agradezco la ayuda de veras pues no tengo ni la menor idea de como plantear el algoritmo. Muchísimas gracias!!