[SciPy-User] Help - scipy.integrate

[Waléria Antunes David]

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 at enthought.com> wrote:

>
>
> 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
>
>
> _______________________________________________
> SciPy-User mailing list
> SciPy-User at 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.html>