... θ
r
r
H
ϕ
−
∂
∂ϕ
r
r
H
θ
iωεs
r
r
r
E
θ
=
1
r sin θ
∂
∂ϕ
(s
r
H
r
) −
1
r
∂
r
r
r
r
H
ϕ
(3.3.2)
iωεs
r
r
r
E
ϕ
=
1
r
∂
r
r
r
r
H
θ
−
∂
∂θ
(s
r
H
r
)
22
where
s
r
= 1 +
σ
iω
σ
∗
=
1
r
r
σdr
r
r
= ... Consequently, the permittivity ε
and the permeability µ are 3 × 3 tensors rather then scalars. Therefore, nine scalar
numb...
... Kronecker-Weierstrass form
{diag(I
r
, O
n r
), diag(D, I
n r
)} with a symmetric and positive define matrix D ∈ R
r r
. Here I
r
and
O
n r
stand for the identity and zero matrices of appropriate ... 05.10.
References
[1] W. Lucht, K. Strehmel, C. Eichler-Liebenow, Linear Partial Differential Algebraic Equation, Report No .18 (1997) 430.
[2] W. Marszalek, Z. Trzaska, A Bounda...
... of predictor–corrector methods 349
44 Order and Stability Barriers 352
440 Survey of barrier results 352
441 Maximum order for a convergent k-step method 353
442 Orderstarsforlinearmultistepmethods ... compensated summation
z = zeros(size(y));
fori=1:n
term = h*f(y) + z;
newy = y + term;
z = term - (newy - y);
y = newy;
end
truncation error, the accumulated rounding error wo...