Изменить работу с массивами (коэффициент корреляции, фазовые портреты)

import sys
import math
import matplotlib.pyplot as plt
import array
import numpy as np
 
N=26000
dt = 0.01
al=7 
bet=9
gam=0
 
def h(x):
    m0 = -(8.0/7.0)
    m1 = -(5.0/7.0)
    return m1*x + 0.5*((m0-m1)*(abs(x+1) - abs(x-1))) 
 
x = [0]*N
y = [0]*N
z = [0]*N
t = [0]*N
xd = [0]*N
yd = [0]*N
zd = [0]*N
x2 = [0]*N
y2 = [0]*N
z2 = [0]*N
x3 = [0]*N
y3 = [0]*N
z3 = [0]*N
x4 = [0]*N
y4 = [0]*N
z4 = [0]*N
masx = []
x[0] = 0.01
y[0] = 0.02
z[0] = 0.15
 
for i in range(0,len(t)):
        t[i]=i*dt
 
for i in range(1,len(x)):
    x[i] = x[i-1] + al*dt*(y[i-1]-x[i-1]-h(x[i-1]))
    y[i] = y[i-1] + dt*(x[i-1]-y[i-1]+z[i-1])
        z[i] = z[i-1] + dt*(gam*z[i-1]-bet*y[i-1])
 
xd[0] = x[0]+0.005 
yd[0] = y[0]+0.005
zd[0] = z[0]+0.005
masx.append(xd[0])
 
for id in range(1,len(xd)):
    xd[id] = xd[id-1] + al*dt*(yd[id-1]-xd[id-1]-h(xd[id-1]))
    yd[id] = yd[id-1] + dt*(xd[id-1]-yd[id-1]+zd[id-1])
        zd[id] = zd[id-1] + dt*(gam*zd[id-1]-bet*yd[id-1])
 
x2[0] = xd[0]+0.005  
y2[0] = yd[0]+0.005
z2[0] = zd[0]+0.005
masx.append(x2[0])
 
for i2 in range(1,len(x2)):
        x2[i2] = x2[i2-1] + al*dt*(y2[i2-1]-x2[i2-1]-h(x2[i2-1]))
        y2[i2] = y2[i2-1] + dt*(x2[i2-1]-y2[i2-1]+z2[i2-1])
        z2[i2] = z2[i2-1] + dt*(gam*z2[i2-1]-bet*y2[i2-1])
 
x3[0] = x2[0]+0.005
y3[0] = y2[0]+0.005
z3[0] = z2[0]+0.005
masx.append(x3[0])
 
for i3 in range(1,len(x3)):
        x3[i3] = x3[i3-1] + al*dt*(y3[i3-1]-x3[i3-1]-h(x3[i3-1]))
        y3[i3] = y3[i3-1] + dt*(x3[i3-1]-y3[i3-1]+z3[i3-1])
        z3[i3] = z3[i3-1] + dt*(gam*z3[i3-1]-bet*y3[i3-1])
 
x4[0] = x3[0]+0.005
y4[0] = y3[0]+0.005
z4[0] = z3[0]+0.005
masx.append(x4[0])
 
for i4 in range(1,len(x4)):
        x4[i4] = x4[i4-1] + al*dt*(y4[i4-1]-x4[i4-1]-h(x4[i4-1]))
        y4[i4] = y4[i4-1] + dt*(x4[i4-1]-y4[i4-1]+z4[i4-1])
        z4[i4] = z4[i4-1] + dt*(gam*z4[i4-1]-bet*y4[i4-1])
 
#plt.xlabel(r'$t$') 
#plt.ylabel(r'$x,y,z$') 
#plt.title(r'$x(t),y(t),z(t)$')
#plt.grid(True) 
#plt.plot(t,x,t,y,t,z)
#plt.plot(t,xd,t,yd,t,zd)
#plt.plot(x,y)
#plt.plot(t,x,t,xd)
#plt.show()
c = np.corrcoef(x,xd)
r=[]
r.append(round(c[0,1],5))
c1 = np.corrcoef(x,x2)
r.append(round(c1[0,1],5))
c2 = np.corrcoef(x,x3)
r.append(round(c2[0,1],5))
c3 = np.corrcoef(x,x4)
r.append(round(c3[0,1],5))
print r
print masx
plt.plot(masx,r)
plt.show()