Anonymous user
Numerical integration/Gauss-Legendre Quadrature: Difference between revisions
Numerical integration/Gauss-Legendre Quadrature (view source)
Revision as of 02:21, 23 August 2018
, 5 years agoAdjust epsilon
(x != x1 is definitely a no-no and not well behaved on some systems) |
(Adjust epsilon) |
||
Line 123:
x1 = x;
x -= lege_eval(N, x) / lege_diff(N, x);
} while ( fdim( x, x1) >
/* fdim( ) was introduced in C99, if it isn't available
* on your system, try fabs( ) */
| |||