自动生成三角形和四面体的积分点和权重
通常人们会参考一本书或一本书来找到单位三角形和四面体的积分点和权重。我正在寻找一种自动计算这些点和权重的方法。以下Mathematica代码示例计算单位线(四/六面体)元素的积分权重和点: unitGaussianQuadraturePoints[points_] := Sort[x /. Solve[Evaluate[LegendreP[points, x] == 0], {x}], ! OrderedQ[N[{#1, #2}]] &]; unitGaussianQuadratureWeights[points_] := Module[{gps, f, int, integr, vars, eqns}, gps = unitGaussianQuadraturePoints[points]; f[0, 0] := 1; f[0., 0] := 1.; f[x_, n_] := x^n; int = Integrate[f[x, #], x] & /@ Range[0, points - 1]; integr = …