For the 3D case, the variable Eddington-Tensor method is used which is accurate
to the order O(u/v) and we neglect acceleration terms
[1]}, [2]}, [3]}, [4]}, [5]}. The system of
equations is given by
\( {1\over c} {\delta F_\nu \over \delta t}+{1\over c} \nabla (u F_\nu + c \nabla (T_\nu E_\nu ))=-k_\nu F_\nu \) 
\(\rho {D(E_\nu /\rho ) \over D t} + E_\nu T_\nu \nabla u - \nabla [ {c T_\nu \over \kappa _nu \rho q_\nu } \nabla (q_\nu E_\nu ) ] = 4 \pi j_\nu - \kappa _\nu \rho c E_\nu \) 

Для случая 3D используется метод переменного тензора Эддингтона, который точен по порядку O(u/v), и мы пренебрегаем ускорениями [1], [2], [3], [4], [5]. Система уравнений задается следующим образом:



\( {1\over c} {\delta F_\nu \over \delta t}+{1\over c} \nabla (u F_\nu + c \nabla (T_\nu E_\nu ))=-k_\nu F_\nu \) 

\(\rho {D(E_\nu /\rho ) \over D t} + E_\nu T_\nu \nabla u - \nabla [ {c T_\nu \over \kappa _nu \rho q_\nu } \nabla (q_\nu E_\nu ) ] = 4 \pi j_\nu - \kappa _\nu \rho c E_\nu \)