open Fmat.Numerics open Fmat.Numerics.MatrixFunctions open Fmat.Numerics.BasicStat let calcPi n = let x = rand [1;n] let y = rand [1;n] let d = x .* x + y .* y let circ = new Matrix(d .< 1.0) let m = sum(circ,1) float(m)/(float)n*4.0 let pi = calcPi 4000000