/*Translated by FOR_C, v3.4.2 (-), on 07/09/115 at 08:33:14 */
/*FOR_C Options SET: ftn=u io=c no=p op=aimnv pf=,p_sbi0k0 s=dbov str=l x=f - prototypes */
#include <math.h>
#include "fcrt.h"
#include <stdio.h>
#include <stdlib.h>
#include "p_sbi0k0.h"
/*     program DRSBI0K0
 *>> 1996-05-28 DRSBI0K0 Krogh Moved formats up.
 *>> 1994-10-19 DRSBI0K0 Krogh  Changes to use M77CON
 *>> 1992-03-18 DRSBI0K0 CLL Added "c" to "program" line above.
 *>> 1990-11-21 WVS (edited by CLL)
 *     Demonstration driver for single precision hyperbolic Bessel
 *     function subprograms.
 *
 *     Compute the Wronskian relation
 *
 *     z = x * (I1(x)*K0(x) + I0(x)*K1(x)) - 1.0
 *
 *     z should be approximately zero.
 *
 *     ------------------------------------------------------------------
 *--S replaces "?": DR?BI0K0, ?BI0K0, ?BI1K1
 *     ------------------------------------------------------------------ */
 
int main( )
{
	long int info, ix;
	float bi0, bi1, bk0, bk1, x, z;
 
 
 
	printf("   X       I0(X)        I1(X)        K0(X)        K1(X)        Z\n   -       -----        -----        -----        -----        -\n");
 
	x = 0.0e0;
	sbi0k0( x, &bi0, &bk0, 1, &info );
	sbi1k1( x, &bi1, &bk1, 1, &info );
	printf(" %4.1f %12.8f %12.8f     INFINITY     INFINITY\n", x, bi0, bi1);
	for (ix = 5; ix <= 50; ix += 5)
	{
		x = ix/10.0e0;
		sbi0k0( x, &bi0, &bk0, -3, &info );
		sbi1k1( x, &bi1, &bk1, -3, &info );
		z = x*(bi1*bk0 + bi0*bk1) - 1.0e0;
		printf(" %4.1f %12.8f %12.8f %12.8f %12.8f  %9.2g\n", x, bi0, bi1, bk0, bk1, z);
	}
	exit(0);
} /* end of function */