0

c 値の範囲をループしたいですか? それ、どうやったら出来るの？外側に for ループまたは while ループを置いてみました:

注: 値が 120 の 4 から 15 の間の c が必要です

for c in range(4,15):
    i=0
    while i< nt-1: #if i<nt-1 is true, the body below will be excuted
        #Update the acceleration using the following equation
        x[i+1] = (-y[i]-z[i])*deltat+x[i]
        #Update the velocity using the follwing equation
        y[i+1] =(x[i]+a*y[i])*deltat+y[i]
        #Update the displacement, using follwing equation
        z[i+1] = (b+z[i]*(x[i]-c))*deltat+z[i]
        #go to next time step
        i=i+1

また

c=arange(4,15,1)
for p in c:
   while i< nt-1: #if i<nt-1 is true, the body below will be excuted
        #Update the acceleration using the following equation
        x[i+1] = (-y[i]-z[i])*deltat+x[i]
        #Update the velocity using the follwing equation
        y[i+1] =(x[i]+a*y[i])*deltat+y[i]
        #Update the displacement, using follwing equation
        z[i+1] = (b+z[i]*(x[i]-p))*deltat+z[i]
        #go to next time step
        i=i+1
4

1 に答える 1

2

すでに NumPy を使用しているので、おそらく次のようにするのが最も自然です。

for c in np.linspace(4, 15, 120):
   ...

このlinspace()呼び出しは 120 個の値を生成します。

In [33]: np.linspace(4, 15, 120)
Out[33]: 
array([  4.        ,   4.09243697,   4.18487395,   4.27731092,
         4.3697479 ,   4.46218487,   4.55462185,   4.64705882,
         4.7394958 ,   4.83193277,   4.92436975,   5.01680672,
        ...
        14.35294118,  14.44537815,  14.53781513,  14.6302521 ,
        14.72268908,  14.81512605,  14.90756303,  15.        ])
于 2012-12-01T17:31:59.270 に答える