Я пытаюсь построить график увеличения температуры в зависимости от угла на приборе с блокировкой приливов anet. Я думал, что это будет просто для l oop, но Python дает мне эту ошибку: TypeError: только скалярные целочисленные массивы могут быть преобразованы в скалярный индекс.
Вот код, который у меня есть.
import matplotlib.pylab as plt
import numpy as np
theta = np.arange(0,(np.pi)/2,100)
L_s = 2.7880945368086e+23
sigma = 5.67e-8
D = 3.76987e+9
A2 = 0.2
X = 8.0000004300E-03
f1 = 1
for i in range(theta):
T_eq=(((1-A2)*(f1*X+1)*L_s*np.cos(theta[i]))/(4*np.pi*sigma*(D**2)))**(1./4.)
print (T_eq)
fig = plt.figure()
plt.plot(theta,T_eq)
plt.grid() #show grid
plt.title('Range of temperature as a function of angle, tidal locked body, Mars analog')
plt.legend(bbox_to_anchor=(1.7,1),loc='upper right',borderaxespad=0.)
plt.xlabel('x') #label axis
plt.ylabel('y') #label axis
plt.show() #show figure
Я не понимаю, в чем проблема. Я повторяю уравнение в диапазоне тета. Как обстоят дела с этими скалярными массивами? Любая помощь будет оценена. Спасибо.