# Mathematica: Problems with integrating and plotting a function

1. Feb 11, 2009

### Bkkkk

Hey, I have been trying to make Mathematica integrate the following function:

$$r_{e}=\frac{c}{H_{0}}\int_{0}^{z}\frac{dz}{\left[\Omega_{NR}(1+z)^{3}+\Omega_n(1+z)^{3-n}+\Omega_{\Lambda}\right]^{\frac{1}{2}}}$$

There are two cases for which this equation needs to be integrated, seperately;

(i) When $$\Omega_{\Lambda}=0$$
(ii) $$\Omega_n=0$$

n is real, so is z and the sum of the density parameters is unity, thus if we take case (ii) $$\Omega_{\Lambda} = 1 - \Omega_{NR}$$

However Mathematica outputs the formatted version of the input line. Next I thought I would just define a function which includes the above integral, and then plot it for varying values of $$\Omega_{NR}$$ and n, however the problem is that often times the Kernel crashes or the curves plotted have gaps in them and/or end unexpectedly.

I would of course prefer to have the integrated function as I presume it would save on plotting time as I need make plots for various values of n for case (i).

Thanks

Last edited: Feb 11, 2009