RuntimeWarning:除算で無効な値が検出されました

| | | | | |

「春のボール」モデルにオイラーの方法を使用してプログラムを作成する必要があります

  from pylab import * from math import * m = 0.1 Lo = 1 tt = 30 k = 200 t = 20 g = 9.81 dt = 0.01 n = int((ceil(t / dt)))km = k / m r0 = [-5,5 * sqrt(3)] v0 = [-5,5 * sqrt(3)] a = zeros((n、2))r = zeros((n、2))v = zeros((n、2))t = zeros((n、2))r [1、 :] = r0 v [1、:] = v0 for i in range(n-1):rr = dot(r [i、:]、r [i、:])** 0.5 a = -g + km * cos(tt)*(rr-L0)* r [i、:] / rr v [i + 1、:] = v [i、:] + a * dt r [i + 1、:] = r [i 、:] + v [i + 1、:] * dt t [i + 1] = t [i] + dt #print norm(r [i、:])plot(r [:、0]、r [: 、1])xlim(-100,100)ylim(-100,100)xlabel( "x [m]")ylabel( "y [m]")show() 

これを取得し続けますエラー:

  a = -g + km * cos(tt)*(rr-L0)* r [i、:] / rr RuntimeWarning:除算で無効な値が見つかりました 

わかりません。コードの何が問題になっていますか?