Stack.

 It's basically an unused cymbal thrown on top of another broken one at the cymbal stand, :P

 #edisiError

Lorenz Attractor in Python with Visual Module.

 So much faster than matplotlib 3d projection

from visual import *
from random import uniform
display(center=(0,0,0),         #pusat display
        background=(1,1,1), 
        #autoscale=False,        #agar display tidak otomatis mengikuti obyek 
        width=600, 
        height=600,
        forward=(-0.4,-0.3,-1)) #arah kamera
x   = 1.
y   = 1.
z   = 1.

dt  = 1./64.
s   = 10.
b   = 8./3.
r   = 28
w   = 0.
dw  = 0.01
    
n   = 0
while 1:
    rate(1)
    while n<3000:
        n   += 1

        xdot    = s * (y-x)
        ydot    = x*r -x*z -y
        zdot    = x*y -b*z

        x       = x+xdot*dt
        y       = y+ydot*dt
        z       = z+zdot*dt

        w       +=dw
        if w>1 or w<0:
            dw   = -dw

        sphere(pos=(x,y,z),radius=.7,color=(0.,w,0))

Timbal balik

 Alfa suka roti bakar "tesoberry" buatan ibuknya atau buatanku (yang hampir pasti selalu gosong, :) ).

 Kadang dia ingin rasa lain.

 Beberapa hari sekali minta roti bakar keju dan "yang coklat untuk adek". Beta memang suka  coklat.

 Jadilah, jika di pagi hari dia pesan sewaktu saya akan berangkat ngampus (saat dia kuantar ke rumah budhe), maka pulang harus bawa roti bakar bandung rasa coklat dan keju terpisah buat mereka.

 Awalnya saya langganan beli  roti bakar di pertigaan yang lumayan dekat dengan rumah.  Langganan, yeah, karena di situ selalu sepi, jarang sekali ngantri.

 Ternyata juga jarang sekali penjualnya ada meski dagangannya ada, entah kenapa.

 Ehm, sepi, ternyata itu. Tetap saja kadang saya rela menunggu lima belas menit untuk orangnya datang. Pulang ke rumah dengan tangan kosong bukan pilihan yang bagus, meski Alfa Beta terkadang juga tidak protes.

 Mengajari mereka menepati janji dengan contoh nyata, :)




 Suatu malam, pulang ngampus, dapat "order" untuk bawa roti-bakar. Ndilalah tukang roti bakar di pertigaan situ tak ada, berikut jualannya. Tutup.

 Terpaksa balik ke arah kampus. Targetnya ke perempatan Galunggung. Dulu sepertinya ada.

 Eh, ternyata selepas perempatan Dieng ternyata ada. Coba ah.

 Kebetulan, hanya ada satu pembeli, mbak cantik yang sedang duduk di sadel bebek matic putih, pesanannya hampir selesai, relatif sepi.

 Sehabis turun dari si ducati langsung disapa

 "Rasa apa mas?"

 "Coklat keju, mas"

 "Dicampur atau pisah?"

 "Pisah"

 "Tunggu sebentar ya"

 "nggih"

 Pesanan sebelumku telah jadi, terbungkus kertas minyak coklat, dan sudah masuk kantong plastik hitam. Entah rasa apa.

 "sudah mbak"

 Si mbak cantik, mengeluarkan uang seratus ribu.

 Tanpa diduga, mas penjual roti minta ijin ke warung tenda sebelah  untuk tukar recehan.

 Wow, kejadian langka. Harusnya, eh, biasanya, penjual akan sedikit protes atau setidaknya tanya "apakah ada uang kecil?'.

 Kesan pertama.

 ...

 Roti panggangku sudah separuh jalan. Tinggal dihias cokelat dan keju.

 Si mas memberi taburan cokelat meses dengan royal, sampai-sampai ketika roti ditutup jadi langsung berubah wujud, obesitas.

 Hal yang sama juga terjadi pada parutan keju. Hampir setengah batang dihbaiskan hanya untuk rotiku.

 Panik, meroh]goh dompet, kulihat daftar harga. Sama, tujuhbelas ribu. What the...

 "Ini mas, " tiba-tiba  kantong plastik berisi roti yang terbungkus kertas minyak cokelat disodorkan ke tanganku.

 Kubayar dengan tiga lembar lima ribuan dan selembar dua ribu, yang diterima dengan senang hati, ada bahan untuk uang kembalian.

 Dan malam itu, di rumah, kami berempat seperti pesta roti mewah, :)

 ...

 Hari-hari berikutnya selalu pesan ke tempat itu.

 Selalu bayar dengan uang pas.

 Selalu di jam yang hampir sama.

 Selalu dengan jaket gunung kuning menyala.

 Selalu bayar dengan uang pas.

 ...

 Suatu malam, pembeli sedang banyak-banyaknya. OK, tak apa, nunggu sambil main-main handphone.

 Satu demi satu pembeli terlayani hingga tinggal dua orang. Sepasang cowok-cewek yang naik motor touring dan aku. Rotinya sudah selesai. Ada satu roti yang juga siap dibungkus

 Ups, lupa pesan.

"Eh mas, coklat keju"

 "Sudah tahu kok, ni tinggal mbungkus "

 Eh.

 Si mas penjual sudah hafal. Mungkin.

...

Beberapa hari setelah itu. Datang ke sana, sengaja diam untuk mengetes.

 Dan..., coklat-keju. Otomatis. Valid, si mas penjual roti telah hafal, hehehe.

 Bakal ke situ terus.

 Sepertinya dihafal oleh seseorang memang bukan apa-apa, tak ada artinya, bagi beberapa orang.

 Tetapi itu juga merupakan kemewahan, meskipun kecil.

 Dan itu juga merupakan nilai lebih, meskipun kecil.

 Juga merupakan nilai tambah, nilai jual, meningkatkan daya saing, meskipun kecil.

 Itu yang membuat saya ke situ terus.

The Cookie Monster dan Sang Pengunyah

Forced to get Force using Center of Mass

 I know. Using solely center of mass as main contributor in force calculation is bad idea in nearly homogenous n-body system.

 Well, I'm curious about that. So, just bite it.

 And, yeah, bad idea.

 The system is fine, generally.

 Each body move as usual, at a glimpse.

 But if we go detail about that. There's strange behavior here and there.


 On previous code with brute force calculation, we could saw, binary system consist of two body orbiting each other. But, with this center-of-mass-code, we see no one. A body that have very close neighbor don't make a twin system. It continues orbiting center of mass.

 As we have system with -body with same mass, it should be natural if two close-body will have greater gravity force than with any other body. It should have "priority"; orbiting each other.

 Ok, this "forced"-center-mass-code is weird, but it'll useful in other system. :)

from visual import *
from random import uniform,random

l       = 17.
dl      = .01
display(center=(0,0,0),background=(1,1,1),autoscale=False, width=600, height=600,forward=(-0.7,-0.7,-1))
distant_light(direction=(1,1,1), color=color.red)
#local_light(pos=(0,0,0), color=(0,0,1))
#sphere(pos=(0,0,0), color=(0,0,0), material=materials.emissive, opacity=0.9)

box(color=color.white, pos=(0,l/2,0),length=l,height=dl, width=l, opacity=0.3)
box(color=color.white, pos=(0,-l/2,0),length=l,height=dl, width=l, opacity=0.3)
box(color=color.white, pos=(l/2,0,0),length=dl,height=l, width=l, opacity=0.3)
box(color=color.white, pos=(-l/2,0,0),length=dl,height=l, width=l, opacity=0.3)
box(color=color.white, pos=(0,0,l/2),length=l,height=l, width=dl, opacity=0.3)

bola    = [] 
n       = 11
dv      = .1        #kecepatan setelah menabrak dinding
G       = 31.
for i in range(n):
    ball            = sphere (pos=(uniform(-7,7),uniform(-7,7),uniform(-7,7)), radius=.2, color=(random(),random(),random()))
    ball.v          = vector(uniform(-1,1),uniform(-1,1),uniform(-1,1))
    
    bola.append(ball)


dt = 1./32.
bola[0].pos=(1,1,1)


def hitungGaya(i,rcm):
    r           = bola[i].pos-rcm
    jarak       = mag(r)
    arah        = -norm(r)
    gaya        = G*1./(pow(jarak,2)+.1)*arah*n
    return gaya
def hitungPusatMassa():
    xcm = 0.
    ycm = 0.
    zcm = 0.
    for i in arange(n):
        xcm += bola[i].x
        ycm += bola[i].y
        zcm += bola[i].z
    xcm /= n
    ycm /= n
    zcm /= n
    return vector(xcm,ycm,zcm)
        
def proses():
    global bola
    rcm = hitungPusatMassa()
    for i in arange(n):
        gaya        = hitungGaya(i,rcm)
        r           = bola[i].pos
        #tambah gaya dari pusat
        '''
        jarak       = mag(bola[i].pos)
        arah        = -norm(bola[i].pos)
        gaya        += 10*G*1./(pow(jarak,2)+.05)*arah
        '''
        a           = gaya
        v           = bola[i].v
        v           += a*dt
        r           += v*dt
        #cek pantul ke tembok
        if bola[i].x > l/2.:
            bola[i].x   = l/2.
            v.x         *= -dv
        elif bola[i].x < -l/2.:
            bola[i].x   = -l/2.
            v.x         *= -dv
        elif bola[i].y > l/2.:
            bola[i].y   = l/2.
            v.y         *= -dv
        if bola[i].y < -l/2.:
            bola[i].y   = -l/2.
            v.y         *= -dv
        if bola[i].z > l/2.:
            bola[i].z   = l/2.
            v.z         *= -dv
        if bola[i].z < -l/2.:
            bola[i].z   = -l/2.
            v.z         *= -dv

        bola[i].pos = r
        
while 1:
    rate (10)
    proses()

.

Menjinakkan Gravitasi.

 Cari cara agar mereka jadi enak untuk dilihat.

 Di kode sebelumnya, bola sudah dicegah agar tidak melarikan diri dengan menambahkan kotak, mereka akan memantul jika menabrak kotak.

 Beres, sekarang tambahkan interaksi antar bola, interaksi sederhana, gravitasi.

 Ok, pake metode brutal, untuk tiap-tiap bola, hitung satu-satu interaksi gravitasi dengan seluruh bola, kemudian jumlahkan.





(ide, cari pusat masa seluruh bola, perlakukan interaksi gaya seperti kode sebelumnya)

 (saat nulis tadi, sambil mikir, ternyata itu ide buruk, tapi boleh juga, lihat saja nanti, bagus atau tidak)

 (tapi tidak sekarang, tidak masuk kode yang sekarang)

 (kode di posting ini tetap pake metode brutal)

 Hasilnya adalah, ..., ok, mereka tetap memantul saat menabrak tembok, tetapi di dalam kotak mereka bergerak dengan sangat cepat, bola liar.

 (yep, metode brutal menghasilkan hasil brutal, tetapi tidak selalu...)

 Ide licik, tambahkan koefisien restitusi di dinding, tumbukan dengan dinding tidak lenting sempurna, kecepatannya tinggal 10% setelah memantul

 Dengan demikian mereka lebih jinak di dalam kotak, :)

from visual import *
from random import uniform,random

l       = 17.
dl      = .01
display(center=(0,0,0),background=(1,1,1),autoscale=False, width=600, height=600,forward=(-0.7,-0.7,-1))
distant_light(direction=(1,1,1), color=color.red)
#local_light(pos=(0,0,0), color=(0,0,1))
#sphere(pos=(0,0,0), color=(0,0,0), material=materials.emissive, opacity=0.9)

box(color=color.white, pos=(0,l/2,0),length=l,height=dl, width=l, opacity=0.3)
box(color=color.white, pos=(0,-l/2,0),length=l,height=dl, width=l, opacity=0.3)
box(color=color.white, pos=(l/2,0,0),length=dl,height=l, width=l, opacity=0.3)
box(color=color.white, pos=(-l/2,0,0),length=dl,height=l, width=l, opacity=0.3)
box(color=color.white, pos=(0,0,l/2),length=l,height=l, width=dl, opacity=0.3)

bola    = [] 
n       = 11
dv      = .1        #kecepatan setelah menabrak dinding
G       = 31.
for i in range(n):
    ball            = sphere (pos=(uniform(-7,7),uniform(-7,7),uniform(-7,7)), radius=.3, color=(random(),random(),random()))
    ball.v          = vector(uniform(-1,1),uniform(-1,1),uniform(-1,1))
    
    bola.append(ball)


dt = 1./32.
bola[0].pos=(1,1,1)


def hitungGaya(i):
    gaya    = vector(0,0,0)
    for j in arange(n):
        if i!=j:  
            r           = bola[i].pos-bola[j].pos
            jarak       = mag(r)
            arah        = -norm(r)
            gaya        += G*1./(pow(jarak,2)+.1)*arah
    return gaya
def proses():
    global bola
    for i in arange(n):
        gaya        = hitungGaya(i)
        #tambah gaya dari pusat
        r           = bola[i].pos
        jarak       = mag(bola[i].pos)
        arah        = -norm(bola[i].pos)
        gaya        += 10*G*1./(pow(jarak,2)+.05)*arah

        a           = gaya
        v           = bola[i].v
        v           += a*dt
        r           += v*dt
        #cek pantul ke tembok
        if bola[i].x > l/2.:
            bola[i].x   = l/2.
            v.x         *= -dv
        elif bola[i].x < -l/2.:
            bola[i].x   = -l/2.
            v.x         *= -dv
        elif bola[i].y > l/2.:
            bola[i].y   = l/2.
            v.y         *= -dv
        if bola[i].y < -l/2.:
            bola[i].y   = -l/2.
            v.y         *= -dv
        if bola[i].z > l/2.:
            bola[i].z   = l/2.
            v.z         *= -dv
        if bola[i].z < -l/2.:
            bola[i].z   = -l/2.
            v.z         *= -dv

        bola[i].pos = r
        
while 1:
    rate (10)
    proses()

.

Bermain dengan Gaya Sentral di Python.

 Merangkak dari kode sebelumnya, dengan tambahan kotak yang mengungkung bola-bola yang mengitari titik asal dengan berbagai kecepatan. Bola-bola tersebut tak bisa menembus kotak, jika menabrak salah satu dinding kotak, maka dia akan memantul

from visual import *
from random import uniform,random

l       = 17.
dl      = .01
display(center=(0,0,0),background=(1,1,1),autoscale=False, width=600, height=600,forward=(-0.7,-0.7,-1))
distant_light(direction=(1,1,1), color=color.red)

box(color=color.white, pos=(0,l/2,0),length=l,height=dl, width=l, opacity=0.3)
box(color=color.white, pos=(0,-l/2,0),length=l,height=dl, width=l, opacity=0.3)
box(color=color.white, pos=(l/2,0,0),length=dl,height=l, width=l, opacity=0.3)
box(color=color.white, pos=(-l/2,0,0),length=dl,height=l, width=l, opacity=0.3)
box(color=color.white, pos=(0,0,l/2),length=l,height=l, width=dl, opacity=0.3)

bola    = [] 
n       =   11

G       = 31.
for i in range(n):
    ball            = sphere (pos=(uniform(1,7),0,uniform(-7,7)), radius=.3, color=(random(),random(),random()))
    ball.v          = vector(0,uniform(1,3),0)
    
    bola.append(ball)


dt = 1./32.
bola[0].pos=(1,1,1)
while 1:
    rate (100)
    for i in arange(n):
        r           = bola[i].pos
        v           = bola[i].v
        jarak       = mag(bola[i].pos)
        arah        = -norm(bola[i].pos)
        gaya        = G*1./(pow(jarak,2)+.1)*arah
        a           = gaya
        v           += a*dt
        r           += v*dt
        #cek pantul ke tembok
        if bola[i].x > l/2.:
            bola[i].x   = l/2.
            v.x         *= -1
        elif bola[i].x < -l/2.:
            bola[i].x   = -l/2.
            v.x         *= -1
        elif bola[i].y > l/2.:
            bola[i].y   = l/2.
            v.y         *= -1
        if bola[i].y < -l/2.:
            bola[i].y   = -l/2.
            v.y         *= -1
        if bola[i].z > l/2.:
            bola[i].z   = l/2.
            v.z         *= -1
        if bola[i].z < -l/2.:
            bola[i].z   = -l/2.
            v.z         *= -1

        bola[i].pos = r
        

.