subroutine single_point_mandelbrot(z,c,zout)
! *************************************************
! *      compute single point mandlebrot          *
! *************************************************
implicit none
integer :: i 
complex*16 :: z, c, zout
do i = 1,100 
  z = z*z + c
  if ((real(z)**2 + aimag(z)**2 ) .gt. 1000.0) then
    exit
  endif
enddo
zout = z 
end subroutine