/*  modified bessel example */


dlibrary ribesl; 
/* 
**  The dlibrary statement can be put into your startup
**   file, and the functions in the .dll can be used just like
**   GAUSS intrinsic functions.  If you have more than one
**   dll in your dll library, list them all in the dlibrary statement.
**   
**   When you put the dlibrary statement into your startup file,
**   you should delete the dlibrary statement from this file.
*/



x = { .1, .2, .3, .4 };

print;
print "modified Bessels of zero and first order";
print mbesseli(x,2,0);

print;
print "modified Bessels of zero order";
print mbesseli0(x);

print;
print "modified Bessels of first order";
print mbesseli1(x);

print;
print "modified exponentially scaled Bessels of zero and first order";
print mebesseli(x,2,0);

print;
print "modified exponentially scaled Bessels of zero order";
print mebesseli0(x);

print;
print "modified exponentially scaled Bessels of first order";
print mebesseli1(x);