How can evaluate the integral (r*bessilei(1,k*r)) dr, k is a constant?

How can I evaluate the integral?

fun = @(r) r.*besseli(1,5.*r);

Answer 1

[Walter Roberson]

syms k r
solution = int(besseli(1, k*r), r)

solution =