| 424 | } |
| 425 | |
| 426 | int nsl_diff_sixth_deriv(const double* x, double* y, const size_t n, int order) { |
| 427 | switch (order) { |
| 428 | case 1: |
| 429 | return nsl_diff_sixth_deriv_first_order(x, y, n); |
| 430 | /*TODO: higher order */ |
| 431 | default: |
| 432 | printf("nsl_diff_sixth_deriv() unsupported order %d\n", order); |
| 433 | return -1; |
| 434 | } |
| 435 | } |
| 436 | |
| 437 | int nsl_diff_sixth_deriv_first_order(const double* x, double* y, const size_t n) { |
| 438 | if (n < 7) |