[SciPy-User] Sequence not work

[Robert Kern]

On Wed, Mar 27, 2013 at 3:11 AM, El suisse <elsuizo37 at gmail.com> wrote:

> Hi
> I have to represent the following sequence:
>
> [image: z[2]=2]
> [image: z[n+1]=2^{n-1/2}\sqrt{1-\sqrt{1-4^{1-n}z[n]^{2}}]  [image:
> n=2,3,4....]
>
>
> and my code is as follows:
>
>
> #!/usr/bin/env python
>
>
> import matplotlib.pyplot as plt
>
> import numpy as np
>
>
> num_muestras = 100
>
> z = np.zeros(num_muestras)
>
>
>
> z[2] = 2.0  #En Python se cuenta desde el cero
>
> for n in range(2,num_muestras-1):
>
>
>     z[n+1] = np.power(2,n-(1/2)) * (np.sqrt(1- np.sqrt(1-np.power(4,1-n) *
> np.power(z[n],2))) )
>

(1/2) should be 0.5 . (1/2) gives 0, not 0.5. The division of integers in
Python 2.x defaults to returning an integer, not a floating point number
(for a variety of reasons that you can Google for if you want an
explanation).

-- 
Robert Kern
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20130327/e0939d80/attachment.html>