2014-09-11から1日間の記事一覧

Haskellで拡散方程式を書いてみた。その一

拡散方程式をHaskellで書いてみた。拡散方程式の陽解法にあるCのプログラムをHaskellに書き直した。diffusion.hs nmax :: Double nmax = 20 n :: Integer n = 20000 dt = 1.0 / 10000.0 dx = 1.0 / nmax a = dt / dx / dx -- initial data u_init = [2.0 * (…