[SciPy-User] Help - scipy.integrate
Warren Weckesser
warren.weckesser@enthought....
Wed Nov 10 08:07:49 CST 2010
On Wed, Nov 10, 2010 at 7:40 AM, Warren Weckesser <
warren.weckesser@enthought.com> wrote:
>
>
> On Wed, Nov 10, 2010 at 6:02 AM, Waléria Antunes David <
> waleriantunes@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
More information about the SciPy-User
mailing list