What about BesselK function for non-integer order? It is a very important function which SAS does not have. Any idea how can we implement it?