Numerical integration/Gauss-Legendre Quadrature: Difference between revisions
Content added Content deleted
(x != x1 is definitely a no-no and not well behaved on some systems) |
(Adjust epsilon) |
||
Line 123: | Line 123: | ||
x1 = x; |
x1 = x; |
||
x -= lege_eval(N, x) / lege_diff(N, x); |
x -= lege_eval(N, x) / lege_diff(N, x); |
||
} while ( fdim( x, x1) > |
} while ( fdim( x, x1) > 2e-16 ); |
||
/* fdim( ) was introduced in C99, if it isn't available |
/* fdim( ) was introduced in C99, if it isn't available |
||
* on your system, try fabs( ) */ |
* on your system, try fabs( ) */ |