zzu数学 实验二圆周率pi的计算

时间:2021-11-11 19:08:29

zzu数学 实验二圆周率pi的计算

示例代码:

Select[Eigenvalues[RandomReal[1, {5, 5}]], Abs[#] < 1 &]
(*挑选特征值小于1的矩阵特征值*)


app[n_] :=
4 N[Length[Select[RandomReal[{-1, 1}, {n, 2}], Dot[#, #] <= 1 &]]/n]
(*dot这里表示的就是点的长度其实因为list是n行两列的*)
TableForm[
Table[n = 10^k;
p = app[n]; {n, p, \[Pi] - p}, {k, 1, 6}],(*给出表格形式,大括号表示要列出的元素*)
TableHeadings -> {{}, {"n", "approximation", "error"}}(*给出行列指标*)]


n = 1000000
points = RandomReal[1, {n, 2}];
4*N[Length[Select[points, Dot[#, #] <= 1 &]]/n]



list = RandomInteger[9, {10, 2}](*产生9内10行2列的向量*)
f = Total[#] < 10 &;(*定义f为字符串*)
Select[list, f](*挑选出和小于10的对子*)
Cases[list, x_ /; f[x]](*挑选将元素代入f中满足表达式的值*)