[SciPy-User] Help - scipy.integrate

[Warren Weckesser]

On Wed, Nov 10, 2010 at 7:40 AM, Warren Weckesser <
warren.weckesser at enthought.com> wrote:

>
>
> On Wed, Nov 10, 2010 at 6:02 AM, Waléria Antunes David <
> waleriantunes at gmail.com> 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
> S?
>
>

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?

Warren
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20101110/df05bb39/attachment.html>