Now both analytical and numerical ways work pretty fine. Next  step is to find out about infinity issue in \int\limits_{0}^{infty} H(s, b) b J_{0} (bq) db