best practices

• Functionalize, deduplicate
• use standard separate .f files from Netlib, especially for Bessel, rather than jamming everything into a huge file
• add jinc(), J'1()
• use associate to make code neater