modified lommel function