python print()の際の見栄えについて

Question

import numpy as np
import matplotlib.pyplot as plt

def f(x,t):
    return  x*(5-x)

a = 0.0
b = 2.0
N = 0.01
h = (b-a)/N
i = 0

t = np.arange(a,b,h)
x_0 = 1

x = np.empty(N)
x[0] = x_0

while i < N-1:
    k_1 = h * f(x[i],t[i])
    k_2 = h * f(x[i] + k_1 /2 , t[i] + h/2 )
    k_3 = h * f(x[i] + k_2 /2 , t[i] + h/2 )
    k_4 = h * f(x[i] + k_3 , t[i] + h )      
    x[i+1] = x[i] + 1/6 * (k_1 + 2*k_2 + 2*k_3 + k_4 ) 
    
    if i == N-1:
      break

    i += 1

plt.plot (t, x)
plt.xlabel("t")
plt.ylabel("x(t)")
plt.show()
```ルンゲクッタ法を用いたコードなのですが、実行すると、
- float' object cannot be interpreted as an integer
  File "C:\Users\hot40\OneDrive\Desktop\runngekyttapra3.py", line 16, in <module>
    x = np.empty(N)

とエラーが出てしまいます。何がいけないのかと、対応策を教えていただきたいです。

Comments

[melian]

N = 0.01 となっていますが、N = 100 の書き間違いでしょうか？

[8524ba23]

ソース中の変数a, b, N, h, i, t, xについて何を表しているのか？という説明コメントを記載ください。
（このあたりに混乱が見られるようです）

[Ackngawe.-]

Nは分割数なので全くその通りです。笑　ありがとうございます

Answer 1

x = np.empty(N)はx = np.empty(0.01)となり、0.01個の未初期化の要素を生成するというコードになってしまっています。
そこは分割区間の総数を与えるべきかと思います。

Comments

[Ackngawe.-]

ありがとうございます！ミスに気が付きました