Integration with a variable integration bound

Hi,

i am using a 2D axisymmetric model. For my simulation i need to evaluate for each point (r,z) the integral of f(r,z')dz' from 0 to z keeping r constant.
Can you tell me how to implement this in comsol?