python シミュレーション　渋滞学

前提

プログラミング初心者です。知識の乏しさゆえ変なことを聞いていたらすいません。。

自分の環境は、macbook pro(2019年購入)、コードはVSCで書いてローカル環境でpython書いてます。

実現したいこと

下記サイトにある動画のようなシミュレーション動画を作りたいです。
https://qiita.com/jabberwocky0139/items/e2526fc5ee3b0dbf144b
下記サイトの数理モデルの派生なので、このサイトを参考にしていますがなかなかうまく作れません。

サイトのコードを参考に作ったのですが以下のようなエラーが出ました。
初心者なのでどういった意味のエラーがわからず、困っています、

発生している問題・エラーメッセージ

Traceback (most recent call last):
 line 42, in <module>
    var = odeint(trafficJam, np.append(x, v), t, args=(K, L), full_output=False)
          ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 line 242, in odeint
    output = _odepack.odeint(func, y0, t, args, Dfun, col_deriv, ml, mu,
             ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 line 21, in trafficJam
    v_dot = K * (V(s) - v)
                 ^^^^
line 13, in V
    return V1 + V2*np.tanh(C1(s-Ic)-C2)
                           ^^^^^^^^
TypeError: 'float' object is not callable

該当のソースコード

python
1import numpy as np
2from scipy.integrate import odeint
3import matplotlib.pyplot as plt
4
5def trafficJam(f, t, K, L):
6    """
7    微分方程式の定義(odeint)
8    f = [x, v]
9    """
10    x, v = f[:N], f[N:]
11
12    def V(s):
13        return V1 + V2*np.tanh(C1(s-Ic)-C2)
14
15
16    x_dot = v
17    x = np.append(x, [x[0]])
18    s = np.diff(x)
19    s = np.array([i if i > 0 else L + i for i in s])
20
21    v_dot = K * (V(s) - v)
22
23    return np.append(x_dot, v_dot)
24
25
26# 各定数（台数、、）
27N, K, L = 30, 0.85, 60
28
29#定数定義
30K, V1, V2, C1, C2, Ic = 0.85, 6.75, 7.91, 0.13, 157, 5.0
31
32# 初期位置
33x = np.arange(N) * L / N
34
35# 初期速度
36v = 6.75 # 最初の1台がブレーキをかける
37
38# 時刻
39t = np.arange(0, 200, 0.1)
40
41# 微分方程式を解く 
42var = odeint(trafficJam, np.append(x, v), t, args=(K, L), full_output=False)
43x_arr, v_arr = var[:, :N], var[:, N:]
44
45
46
47
48
49import matplotlib.animation as animation
50
51r = L / (2 * np.pi)
52fig = plt.figure(figsize=(6, 6))
53
54ims = []
55for index, x in enumerate(x_arr):
56    theta = x / r
57    im = plt.plot(r * np.cos(theta), r * np.sin(theta), '.k', markersize=12)
58    ims.append(im)
59
60Writer = animation.writers['ffmpeg']
61# writer = Writer(fps=15, bitrate=1800)
62writer = Writer(fps=60)
63ani = animation.ArtistAnimation(fig, ims, interval=10)
64
65plt.legend()
66ani.save("output.mp4", writer=writer)
67plt.show()