Question: BesselJ Integration

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

(1/48)*BesselJZeros(0, 1)^2*hypergeom([3/2, 2], [1, 3, 4], -BesselJZeros(0, 1)^2)

I know that an expression without the hypergeom function exists. I can calculate it in Mathematic and get

mathematica:=(-BesselJ(1,BesselJZeros(0,1))^2*(-8+BesselJZeros(0,1)^2))/(6*BesselJZeros(0,1)^2);

Is there any way to get Maple to give me the answer without the hypergeom functions? It would be great if there was a way to expand it into the format without the hypergeom.

Thank you.