Go to the documentation of this file.
9 #include "../../Exceptions/error.h"
16 void ratint(
double *xa,
double *ya,
int n,
double x,
double *y,
double *dy)
19 double w,t,hh,h,dd,*c,*d;
24 for (i=0;i<=n-1;i++) {
47 throw error_msg(
"RATINT_ERROR",
"Error in routine ratint");
53 *y += (*dy=(2*ns < (n-
m) ? c[ns+1] : d[ns--]));