[SciPy-User] Help - scipy.integrate
Waléria Antunes David
waleriantunes@gmail....
Wed Nov 10 08:25:45 CST 2010
Yeah, i'm considering omega_k = 0 ..yeah i implemented equation (12) :
http://pastebin.com/SHTKcMNU ... line 47
Thanks,
On Wed, Nov 10, 2010 at 12:07 PM, Warren Weckesser <
warren.weckesser@enthought.com> wrote:
>
>
> 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
>
>
> _______________________________________________
> SciPy-User mailing list
> SciPy-User@scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.scipy.org/pipermail/scipy-user/attachments/20101110/f7d9f40b/attachment-0001.html
More information about the SciPy-User
mailing list