拡散方程式をHaskellで書いてみた。拡散方程式の陽解法にあるCのプログラムをHaskellに書き直した。diffusion.hs import qualified Data.Vector.Unboxed as D nmax :: Double nmax = 200 dt :: Double dt = 1.0 / 100000.0 dx :: Double dx = 1.0 / nmax a :…

拡散方程式を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 * (…