((Subscript[Ϛ, a,k])≈N(0,(Subscript[Σ, Ϛ a,k]^2)),
(Subscript[Σ, Ϛ a,k]^2)=[((Subscript[σ, Ϛ ax]^2)000(Subscript[σ, Ϛ ay])000(Subscript[σ, Ϛ az]))]+((μ^2)/3)((\[LeftBracketingBar]Δ (Subscript[z, a,k])\[RightBracketingBar])^2)(Subscript[I, 3×3])(37))