/*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_sbesj0 s=dbov str=l x=f - prototypes */
#include <math.h>
#include "fcrt.h"
#include <stdio.h>
#include <stdlib.h>
#include "p_sbesj0.h"
/*     DRSBESJ0
 *>> 2009-10-26 DRSBESJ0 Krogh Moved function calls from print statement.
 *>> 1996-05-30 DRSBESJ0 Krogh Added external statement.
 *>> 1994-10-19 DRSBESJ0 Krogh  Changes to use M77CON
 *>> 1994-09-13 DRSBESJ0 CLL Typed all variables.
 *>> 1994-09-01 DRSBESJ0 WVS Moved formats to top for C conversion
 *>> 1994-08-09 DRSBESJ0 WVS set up for CHGTYP
 *>> 1992-04-29 DRSBESJ0 CAO Replaced '1' in format.
 *>> 1987-12-09 DRSBESJ0 Lawson  Initial Code.
 *--S replaces "?": ?BESJ0, ?BESJ1, ?BESY0, ?BESY1, DR?BESJ0
 *     Demonstration driver for Bessel functions.
 *
 *     Z = (PI / 2) * X * ( J1(X)*Y0(X) - J0(X)*Y1(X) ) - 1.0
 * */
 
int main( )
{
	long int ix;
	float bej0, bej1, bey0, bey1, x, z;
	static float pi2 = 1.5707963267948966192313216e0;
 
 
 
 
	printf("   X       J0(X)        J1(X)        Y0(X)        Y1(X)        Z\n   -       -----        -----        -----        -----        -\n");
 
	x = 0.e0;
	bej0 = sbesj0( x );
	bej1 = sbesj1( x );
	printf(" %4.1f %12.7f %12.7f    -INFINITY    -INFINITY\n", x, bej0, bej1);
 
	for (ix = 5; ix <= 50; ix += 5)
	{
		x = ix/10.e0;
		bej0 = sbesj0( x );
		bej1 = sbesj1( x );
		bey0 = sbesy0( x );
		bey1 = sbesy1( x );
		z = pi2*x*(bej1*bey0 - bej0*bey1) - 1.0e0;
		printf(" %4.1f %12.7f %12.7f %12.7f %12.7f  %9.2g\n", x, bej0, bej1, bey0, bey1,
		   z);
	}
	exit(0);
 
} /* end of function */