I am trying to integrate the following expression in Maple. It produces the correct result, but not in the format that I want.

BesselProblem := BesselJ(0, BesselJZeros(0, 1)*r)*BesselJ(2, BesselJZeros(0, 1)*r)*r^3;

int(BesselProblem, r = 0 .. 1);

returns the expression