I have this function:

f(r,z) = r*(1 + g(r,z))

The function g(r,z) is a modification of the student t-distribution in z, where the degrees of freedom depend continuously on r.

I would like to integrate this function f(r,z) from 0 to R with respect to r.

Unfortunately Maple could not find the integral, so I would like to do this numerically in MATLAB.

I have no clue where to start though. I've integrated a function of one variable in MATLAB before, but only with constant limits.

Any help would be greatly appreciated!