3. Approximation of an integral in C. Consider the integral from 0 to xmax of the shifted Pareto density, f(x) = (p-1) cp-1 (x+c)-p, for x 0, and f(x) = 0 otherwise, where c>0 and p>1 are parameters. Let c = 5 and p = 3. Write a C function called paretoint(xmax,c,p) that approximates this integral over a grid of 1 million values ranging from x = 0 to xmax. Note that technically paretoint() is not only going to be a function of xmax, c, and p, but will also have another input variable which will store the result. Call your C function