Numerical integration/Gauss-Legendre Quadrature: Difference between revisions
Numerical integration/Gauss-Legendre Quadrature (view source)
Revision as of 22:10, 2 August 2021
, 2 years ago→{{header|Mathematica}}
(→{{header|J}}: Use more descriptive names and tacit) |
|||
Line 1,116:
{{out}}<pre>20.035577718386</pre>
=={{header|Mathematica}}/{{header|Wolfram Language}}==
code assumes function to be integrated has attribute Listable which is true of most built in Mathematica functions
<lang Mathematica>gaussLegendreQuadrature[func_, {a_, b_}, degree_: 5] :=
Line 1,123:
weights = 2 (1 - nodes^2)/(degree LegendreP[degree - 1, nodes])^2;
(b - a)/2 weights.func[(b - a)/2 nodes + (b + a)/2]]
gaussLegendreQuadrature[Exp, {-3, 3}]</lang>
{{out}}<pre>20.0356</pre>
|