Mandelbrot set (fractal) in python

Iseng ngerjain exercise, he..

Mandelbrot set itu apa? Set dari bilangan kompleks c yang kalau diurutkan menjadi (c, c^2 + c, (c^2+c)^2 + c, ((c^2 + c)^2 + c)^2 + c, \ldots ) tidak akan mencapai infinit. :D

Dalam perhitungan dapat di tulis: z_{n} = z_{n}^2 + c .

import numpy as np
import matplotlib.pyplot as plt

def mandelbrot(N_max, threshold, nx, ny):
    # a grid of c-values
    x = np.linspace(-4, 4, nx)
    y = np.linspace(-4, 4, ny)

    c = x[:,np.newaxis] + 1j*y[np.newaxis,:] # Initial complex set

    # iteration
    z = c
    for j in xrange(N_max):
        z = z**2 + c

    mandelbrot_set = (abs(z) < threshold)

    return mandelbrot_set

# let's iterate over N_max
for i in xrange(30):
    mandelbrot_set = mandelbrot(i, 50., 401, 401)
    plt.imshow(mandelbrot_set.T, extent=[-4, 4, -4, 4])
    plt.gray()
    plt.savefig('mandelbrot_'+str(i)+'.png')

N_max = 6
mandelbrot_6

N_max = 25
mandelbrot_25

Ayo ubah jadi gif, hehe

 convert -delay 10 $(ls -v mandelbrot_*.png) output.gif 

output

Python for Scientist (mathematics, science, and engineering)

Learning Resources

Tutorial

Book

  • Python Scripting for Computational Science – Hans Petter Langtangen

Ipython notebook

of course by default: Google and Stackoverflow

Important Packages

Intergrated python packages

Specific Fields
Astronomy & Astrophysics

  • Astropy –  a community python library for astronomy
  • Ureka – a collection of useful astronomy software that is generally centered around Python and IRAF

 Please add more in comment!

Tikz by Example

Dalam beberapa bulan terakhir saya belajar menggambar di latex (tikz atau pgf), untuk keperluan penulisan buku dan kepuasan batin.. hahah..

Saya kumpulkan beberapa hasil karya tersebut yang kebanyakan code-nya bersumber dari contekan di stackoverflow.. silahkan jika perlu mencontek dan/atau menambahkan hasil karyanya ke sini:

https://github.com/ridlo/tikz_by_example

(termasuk saya sendiri.. hahah :D ).

Contoh:
eksentrisAnomali

hiperbola

solar_wind

sunset

kepler1

Solusi OSK Astronomi 2015

Kemarin makan malam, eh di akhir malah ngebahas soal (toweng), karena banyak anak2 haus ilmu yang bertanya di WA dll :v,  jadi akhirnya kami buat solusi OSK Astronomi 2015… LOL.

Jika ingin sekaligus menjajal kemampuan lagi, bisa coba versi online berikut (numpang di salah satu lembaga bimbel olim, hehe)

Kalau mau langsung download, silahkan: link download
Terdapat update yang dilakukan, terakhir tanggal 14 Februari pukul 00:30,
jika ada yang pernah mendownload sebelumnya sebaiknya download kembali versi terbaru.

Kalau dirasa masih ada yang salah bisa dimasukkan sebagai komentar disini! :D