/* nag_specfun_bessel_i1_scaled (s18cfc) Example Program.
*
* Copyright 2024 Numerical Algorithms Group.
*
* Mark 30.2, 2024.
*/
#include <nag.h>
#include <stdio.h>
int main(void) {
Integer exit_status = 0;
double x, y;
/* Skip heading in data file */
scanf("%*[^\n]");
printf("nag_specfun_bessel_i1_scaled (s18cfc) Example Program Results\n");
printf(" x y\n");
while (scanf("%lf", &x) != EOF)
{
/* nag_specfun_bessel_i1_scaled (s18cfc).
* Scaled modified Bessel function exp(-|x|) I_1(x)
*/
y = nag_specfun_bessel_i1_scaled(x);
printf("%12.3e%12.3e\n", x, y);
}
return exit_status;
}