search

vavachan / brownian_motion

0 stars 0 forks source link

fluctuations don't die out. #3

Closed vavachan closed 9 years ago

vavachan commented 9 years ago

Omega=1 gamma=4 temp=1 using euler maruyama omega_1_gamm_4_em Omega=1 gamma=4 temp=1 using Runge kutta omega_1_gamm_4_rk Omega=2 gamma=1 temp=1 using euler maruyama omega_2_gamm_1_em Omega=1 gamma=4 temp=1 using RK omega_2_gamma_1_rk Omega=2 gamma=1 temp=2 using RK omega_2_gamma_1_temp_2_rk In all cases I used dt =.005 , The saturation value is (Kb_T)/(m_omega^2) Since m=1 , kb=1 , for T=1 and omega = 2 , the saturation value is 1/4 but the simulation give 1/2 , but it scales according to the formula.