Ошибка типа: ошибка только для массивов размера 1 в алгоритме Ньютона-Рафсона

Question

Я работаю над научным расчетом в Python и использую алгоритм Ньютона-Рафсона, но столкнулся с проблемой и, к сожалению, не смог ее решить.

Я выполняю этот код в Spyder (Python3.6), и моя машина - Windows 10.

N=361

t3 = mt.pi/4   #Initial Guess for Newton-Raphson
t4 = mt.pi/2   #Initial Guess for Newton-Raphson

for i in range(1,N):

    theta2[i-1] = (i-1)*(2*mt.pi)/(N-1)

    #Newton-Raphson Calculation
    for j in range(1,6):
        phi11 = a*mt.cos((i-1)*(2*mt.pi)/(N-1)) + b*mt.cos(t3) - c*mt.cos(t4) - d
        phi12 = a*mt.sin((i-1)*(2*mt.pi)/(N-1)) + b*mt.sin(t3) - c*mt.sin(t4)

        phi = np.array([[phi11], [phi12]])

После этого вычисления я переопределяю массив theta2 с помощью цикла for.