Вместо доступа к E_R
напрямую из функции попробуйте следующее:
def function():
global y1
global y2
y1 = list(np.array([E_R, E_theta, E_z]))
y2 = list(np.array([B_R, B_theta, B_z]))
return y1, y2
print(function())
print(y1[1])
Это приведет к доступу к E_R
. Просто получите доступ к любому элементу по индексу, в котором он появляется. Возможно, это не самое красивое решение и не самый Pythoni c, но оно работает.