Declare Function y0 CDECL Lib "msvcrt.dll" Alias "_y0" (ByVal x As Double) As Double
Declare Function y1 CDECL Lib "msvcrt.dll" Alias "_y1" (ByVal x As Double) As Double
Declare Function yn CDECL Lib "msvcrt.dll" Alias "_yn" (ByVal n As Integer, ByVal x As Double) As Double
Declare Function j0 CDECL Lib "msvcrt.dll" Alias "_j0" (ByVal x As Double) As Double
Declare Function j1 CDECL Lib "msvcrt.dll" Alias "_j1" (ByVal x As Double) As Double
Declare Function jn CDECL Lib "msvcrt.dll" Alias "_jn" (ByVal n As Integer, ByVal x As Double) As Double
Uses "console"
Dim c, n As Long
Dim x As Ext
x = 2.387
n = 3
Console_WriteLine "Bessel functions for x = "&Format$(x, 16 )
Console_WriteLine " Kind   Order  Function     Result" 
Console_WriteLine " First  0      j0( x )     "&Format$(j0( x ), 16 )
Console_WriteLine " First  1      j1( x )     "&Format$(j1( x ), 16 )
For c = 2 To 4
  Console_WriteLine " First  "&LTrim$(Str$(c))&"      jn("&Str$(c)&", x )  "&Format$(jn( c, x ), 16 )
Next
Console_WriteLine " Second 0      y0( x )     "&Format$(y0( x ), 16 )
Console_WriteLine " Second 1      y1( x )     "&Format$(y1( x ), 16 )
For c = 2 To 4
  Console_WriteLine" Second "&LTrim$(Str$(c))&"      yn( "&Str$(c)&", x )  "&Format$(yn( c, x ), 16 )
Next
Console_WriteLine "All done. Press any key to finish"
Console_WaitKey