[SciPy-User] Help - scipy.integrate
Wed Nov 10 08:07:49 CST 2010
On Wed, Nov 10, 2010 at 7:40 AM, Warren Weckesser <
> On Wed, Nov 10, 2010 at 6:02 AM, Waléria Antunes David <
> firstname.lastname@example.org> wrote:
>> Hi all,
>> I have this equation: http://img52.imageshack.us/i/equao.jpg/
>> I need to accomplish the integration of it and then plot on a graph ... I
>> managed to do the following: http://pastebin.com/kkS3EX0m .....but i
>> don't know if i'm using the function integrate correctly.
>> See the line 22 and def func (line 33).
> In your equation it appears that the function S is applied to
> sqrt(abs(omega_k))*(the integral), but I don't see S in your code. What is
Partially answering my own question... google found what I think is your
reference. It looks like you are computing equation (11) from this paper:
http://arxiv.org/abs/astro-ph/0402512. So S is "sinn", where "sinn” is sinh
for omega_k > 0 and sin for omega_k < 0. The formula simplifies to your
code if omega_k = 0. Are you only considering omega_k = 0?
Also, the plot that you showed is mu vs. z, not d_L vs. z. Have you also
implemented equation (12) somewhere?
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the SciPy-User