Как изменить значения вектора?

Question

[Игорь]

Объявил переменную и задал для нее значения:

dphi = x[0:2] * 0.1
      
p     = np.zeros((4)) 
dp    = np.zeros((4))
      
dp[0] = (-dphi[0] * p[1] - dphi[1] * p[2] - dphi[2] * p[3])/2
dp[1] = ( dphi[0] * p[0] + dphi[2] * p[2] - dphi[1] * p[3])/2
dp[2] = ( dphi[1] * p[0] - dphi[2] * p[1] + dphi[0] * p[3])/2
dp[3] = ( dphi[2] * p[0] + dphi[1] * p[1] - dphi[0] * p[1])/2

p = p + dp

При запуске программа выдает ошибку:
dp[0] = (-dphi[0] * p[1] - dphi[1] * p[2] - dphi[2] * p[3])/2
IndexError: index 2 is out of bounds for axis 0 with size 2
Что не так и как это исправить?

Comments

[Игорь]

Переменная х
J = np.array([1400, 830, 730])
My = np.array([0, 0, 0])
Mv = np.array([0, 0, 0])
M = My + Mv
x = np.zeros((3))
x[0] = (M[0] - (J[2]-J[1])*x[1]*x[2])/J[0]
x[1] = (M[1] - (J[0]-J[2])*x[0]*x[2])/J[1]
x[2] = (M[2] - (J[1]-J[0])*x[0]*x[1])/J[2]

Answer 1

[Andrey Druzhaev]

Вы пытаетесь использовать 3 элемент массива dphi в выражении dphi[2] * p[3]. А создали массив только с 2-мя элементами: dphi = x[0:2] * 0.1.

Судя по всему надо сделать так: dphi = x[: 3] * 0.1.