2017-06-30 17 views
0
from numpy import zeros,linspace 

N = 100 
points = N**2 
x1,x2 = -2,2 
y1,y2 = -2,2 

m_array = zeros([N,N],float) 

i,j = -1,0 

for x in linspace(x1,x2,N): 
    i += 1 
    for y in linspace(y1,y2,N): 
     if x == 0 and y == 0: 
      continue 
     else: 
      c = complex(x,y)   
      z = 0 
      if abs(c)<2: 
       for k in range(0,101): 
        zprime = abs(z) + c 
        z = zprime 
       if abs(z) < 2: 
        m_array[i,j] = 1 
        j += 1 
        continue 
       else: 
        j += 1 
        continue 
      else: 
       j += 1 
       continue 

N.B. 'k'の範囲は任意ですが、比較的大きくする必要があります。Python 3エラー:「IndexError:インデックス140は、サイズ1の軸1の範囲外です」

私はウェブサイト上の以前の質問のいくつかを読んだことがありますが、私自身のコードで問題を見つけることはできません。これは密度プロットを使用してマンデルブロ集合をプロットする私の試みです。

答えて

0

ループはfor y in linspace(y1,y2,N):ループに続行します。ループはありません。ループはfor x in linspace(x1,x2,N):です。

if abs(z) < 2: 
    m_array[i,j] = 1 
    j += 1 
    continue 

したがって、jは、迅速な応答を超える99

+0

感謝することができます。あなたの助けを借りて、ちょっとしたやり直しの後、私はそれを働かせました。 –

+0

@DominicWilliamson問題ありません。あなたも私にピースアップアップしてもらえますか? – lf215

+0

こんにちは、私はしますが、私はまだそうするには十分な評判がありません:( –

関連する問題