Skip to main content
Article
Licensed
Unlicensed Requires Authentication

High order weak approximation for irregular functionals of time-inhomogeneous SDEs

  • EMAIL logo
Published/Copyright: February 20, 2021

Abstract

This paper shows a general weak approximation method for time-inhomogeneous stochastic differential equations (SDEs) using Malliavin weights. A unified approach is introduced to construct a higher order discretization scheme for expectations of non-smooth functionals of solutions of time-inhomogeneous SDEs. Numerical experiments show the validity of the method.

Funding statement: This work is supported by JSPS KAKENHI (grant number 19K13736), MEXT, Japan, and JST PRESTO (grant number JPMJPR2029), Japan.

A Proof of Lemma 3.7

For k and for multi-indices αe{0,1,,d}re, re, e=1,,k such that αe2, e=1,,k, we have

e=1k𝔹αe,t,sStrat=σSi=1kri,σ(1)<<σ(r1),σ(r1+1)<<σ(r1+r2),,σ(i=1k-2ri+1)<<σ(i=1k-1ri),σ(i=1k-1ri+1)<<σ(i=1kri)𝔹σ-1(α1αk),t,sStrat.

By Proposition 3.5, it holds that

δy(Bt,s),e=1k𝔹αe,t,sStrat𝔻𝔻-
=β=σ-1(α1αk)pγ,p,σSi=1kri,σ(1)<<σ(r1),σ(r1+1)<<σ(r1+r2),σ(i=1k-2ri+1)<<σ(i=1k-1ri),σ(i=1k-1ri+1)<<σ(i=1kri)12η(β,γ)δy(t,s),𝔹γ,t,sItô𝔻𝔻-.

By integration by parts, we have

δy(Bt,s),𝔹γ,t,sItô𝔻𝔻-=1|γ|!(s-t)|γ|γδy(Bt,s),1𝔻𝔻-
=1|γ|!(s-t)|γ|-|γ|δy(Bt,s),𝔹γ,t,sSkor𝔻𝔻-.

Then we obtain the assertion.

B Proof of Lemma 3.8

By integration by parts, we have

IS(X¯st,x),𝔹γ,t,sSkor𝔻𝔻-=S(X¯st,x),HI(X¯st,x,𝔹γ,t,sSkor)𝔻𝔻-.

The computation of HI(X¯st,x,𝔹γ,t,sSkor) is reduced to the computation of Skorohod integrals due to the following:

H(i)(X¯st,x,𝔹γ,t,sSkor)=(s-t)-1j{1,,d}M(i)(j)(t,x){(Bsj-Btj)𝔹γ,t,sSkor-tsDj,u𝔹γ,t,sSkor𝑑u}
=(s-t)-1j{1,,d}M(i)(j)(t,x)𝔹(j)γ,t,sSkor.

Then we obtain the assertion.

C Proof of Lemma 3.9

Instead of the standard integration by parts on the Wiener space, we apply a Bismut type formula when N=d:

i=1diS(X¯st,x),Gi𝔻𝔻-=i,j=1dS(X¯st,x),1s-t[V-1]ji(t,x){Gj(Bsi-Bti)-tsDi,uGjdu}𝔻𝔻-

for G=(G1,,Gd)(𝔻)d. Iteratively using this formula, we obtain the assertion.

D Proof of Lemma 4.1

In the proof, we use the method of integration by parts with uniform non-degeneracy of the Malliavin covariance matrix in the sense of Watanabe. We introduce the scaling SDE with a parameter ε(0,1]:

Xst,x,(ε)=x+ε2tsV0(u,Xut,x,(ε))𝑑u+εi=1dtsVi(u,Xut,x,(ε))𝑑Bui,

with Xtt,x,(ε)=xN. We define

Fε(t,s,x)=(Xst,x,(ε)-x-ε2V0(t,x)(s-t))/ε,

which is expanded by using the Stratonovich–Taylor expansion (Kloeden and Platen (1999)) as

Fε(t,s,x)=i=1dVi(t,x)Bt,si+k=2e|α|=kεα-1V^α1V^αr-1Vαr(t,x)𝔹α,t,sStrat+re+1ε(t,s,x),

where re+1ε(t,s,x) is the residual. Here, a family of Wiener functionals Fε(t,s,x), ε(0,1], is uniformly non-degenerate and it holds that

E[φ(Xst,x)]=Nφ(x+ε2V0(t,x)(s-t)+εy)δy(Fε(t,s,x)),1𝔻𝔻-dy|ε=1.

We put F0(t,s,x):=i=1dVi(t,x)Bt,si. By Taylor expansion and integration by parts, we have

δy(Fε(t,s,x)),1𝔻𝔻-
=δy(F0(t,s,x),1𝔻+j=12m+1εjk=1ji=1kbi=j+k,bi2,i=1,,kI=(I1,,Ik){1,,N}k1k!𝔻-
   ×δy(F0(t,s,x),HI(F0(t,s,x),ν=1kαν{0,1,,d}rν,αν=bν,rνV^α1νV^αrν-1νVαrννIν(t,x)𝔹αν,t,sStrat)𝔻𝔻-
   +ε2m+201(1-λ)2m+1(2m+1)!δy(Fελ(t,s,x)),Rmελ(t,s,x)𝔻𝔻-dλ,

where

Rmη(t,s,x)=(2m+2)k=12m+2i=1kbi=2m+2+k,bi2,i=1,,kI=(I1,,Ik){1,,N}k1k!
×HI(Fη(t,s,x),i=1kEbi,Iiη(t,s,x)),

with Wiener functionals Ebi,Iiη(t,s,x)𝔻, i=1,,k, satisfying that for , p>1,

supη[0,1],xNEbi,Iiη(t,s,x),p=O((s-t)bi/2).

Then it holds that

E[φ(Xst,x)]=E[φ(X¯st,x)]+j=12m+1k=1ji=1kbi=j+k,bi2,i=1,,kI=(I1,,Ik){1,,N}k1k!
×E[φ(X¯st,x)HI(X¯st,x,ν=1kαν{0,1,,d}rν,αν=bν,rνV^α1νV^αrν-1νVαrννIν(t,x)𝔹αν,t,sStrat)]
+01(1-η)2m+1(2m+1)!E[φ(X^sη,t,x)Rmη(t,s,x)]𝑑η,

where

X^sλ,t,x=x+V0(t,x)(s-t)+ηFη(t,s,x),η[0,1].

The estimate of the weight Rmη(t,s,x) in the error term is reduced to the following: for I{1,,N}k, b1,,bk2 such that i=1kbi=2m+2+k, for p>1, there is C>0 such that for xN, t<s and η[0,1],

HI(Fη(t,s,x),i=1kEbi,Iiη(t,s,x))pC(s-t)(2N-1/2)|I|-Ndet(σFη(t,s,x))-1(8N+4)|I|(|I|+1)p(|I|+1)(2|I|-1)
×DFη(t,s,x)|I|,(16N+8)N|I|2,H2N|I|(2|I|-1)i=1kEbi,Iiη(t,s,x))|I|,(4N+2)|I|p

by [9, Corollary 3.7]. Note that DFη(t,s,x)=DXsη,t,x/η. The factor DFη(t,s,x),q,H is bounded, for and q>1, by

supη[0,1],xNDFη(t,s,x),q,H=O((s-t)1/2)

by [9, Theorem 2.19]. Also, by [9, Theorem 3.5], for q>1 we have

supη[0,1],xNdet(σFη(t,s,x))-1q=O((s-t)-N).

Since for and q>1,

supη[0,1],xNi=1kEbi,Iiη(t,s,x)),q=O((s-t)(2m+2+k)/2),

we have that for p>1 there exists C>0 such that

supη[0,1],xNHI(Fη(t,s,x),i=1kEbi,Iiη(t,s,x))pC(s-t)-|I|/2(s-t)(2m+2+|I|)/2=C(s-t)m+1

for t<s. Therefore, there is C>0 such that

supxN|E[φ(Xst,x)]-{E[φ(X¯st,x)]+j=12m+1k=1j{bi}i=1k,i=1kbi=j+k,bi2,i=1,,kI=(I1,,Ik){1,,N}k1k!
   ×E[φ(X¯st,x)HI(X¯st,x,ν=1kαν{0,1,,d}rν,αν=bν,rνV^α1νV^αrν-1νVαrννIν(t,x)𝔹αν,t,sStrat)]}|
Cφ(s-t)m+1

for all φCb(N) and t<s.

E Proof of Lemma 4.2

By applying Lemma 3.7 and Lemma 3.8 with (2.3) to the expansion in Lemma 4.1, we have

E[φ(Xst,x)]=𝒯t,sφ(x)+t,s1φ(x),

where

𝒯t,sφ(x)=E[φ(X¯st,x)]+j=12m+1k=1jb={bi}i=1k,i=1kbi=j+k,bi2,i=1,,kI=(I1,,Ik){1,,N}k1k!αe=(α1e,,αree){0,1,,d}re,αe=be,re,e=1,,k
×e=1kV^α1eV^αre-1eVαreeIe(t,x)β=σ-1(α1αk)pγ,p,σSi=1kri,σ(1)<<σ(r1),σ(r1+1)<<σ(r1+r2),σ(r1++rk-1+1)<<σ(r1++rk)12η(γ,β)1|γ|!(s-t)|γ|-|γ|
×κ{1,,d}|I|(s-t)-|I|MIκ(t,x)E[φ(X¯st,x)𝔹κγ,t,sSkor]

and t,s1φ satisfies that there exists C>0 such that

t,s2φCφ(s-t)m+1for t<s.

If we use the smoothness of φ, the expectations of type E[φ(X¯st,x)𝔹κγ,t,sSkor] with length |γ|m+1 in 𝒯t,sφ(x) will be of O((s-t)m+1) due to the following identity: for |γ|m+1,

(s-t)|γ|-K(γ)1|γ|!κ{1,,d}|I|(s-t)-|I|MIκ(t,x)E[φ(X¯st,x)𝔹κγ,t,sSkor]
=(s-t)|γ|-K(γ)1|γ|!E[Iφ(X¯st,x)𝔹γ,t,sSkor]
(E.1)=(s-t)|γ|1|γ|!α{1,,N}|γ|k=1|γ|Vγkαk(t,x)E[Iαφ(X¯st,x)],

where the calculation in the proof of Lemma 3.7 and the duality formula (2.1) are applied. Then we define

t,s2φ(x)=𝒯t,sφ(x)-E[φ(X¯st,x))𝒲(m)(t,s,x,Bt,s)],

i.e.

E[φ(Xst,x)]-E[φ(X¯st,x)𝒲(m)(t,s,x,Bt,s)]=t,s1φ(x)+t,s2φ(x).

Here, in E[φ(X¯st,x)𝒲(m)(t,s,x,Bt,s)], expectations of type E[φ(X¯st,x)𝔹κγ,t,sSkor] with length |γ|m+1 are all removed.

By (E.1), we see that t,s2φ(x) has the form

t,s2φ(x)=(s-t)m+1νC()(t,x)E[β()φ(X¯st,x)]

for some ν=ν(m), multi-indices β() and bounded smooth functions C():[0,T]×N with bounded derivatives, =1,,ν. Note that (E.1) is equivalent to the following:

(s-t)|γ|-K(γ)1|γ|!κ{1,,d}|I|(s-t)-|I|MIκ(t,x)E[φ(X¯st,x)𝔹κγ,t,sSkor]
=(s-t)|γ|-K(γ)1|γ|!Nφ(y)Iδy(X¯st,x),𝔹γ,t,sSkor𝔻𝔻-dy
=(s-t)|γ|1|γ|!α{1,,N}|γ|k=1|γ|Vγkαk(t,x)Nφ(y)Iαδy(X¯st,x),1𝔻𝔻-dy.

Then we have the assertion.

F Proof of Lemma 4.3

Since 𝒲(m)(t,s,x,Bs-Bt) is given by the sum of the form

v(t,x)(s-t)q1𝒫q2(Bt,s),q1+q22,q1,q2,

where vCb([0,T]×N) and 𝒫r(ξ) is polynomial of ξ=(ξ1,,ξd) of degree r, there is C>0 such that

|Qt,s(m)f(x)|f𝒲(m)(t,s,x,Bt,s)2f(1+C(s-t))

for all t<s and xN, where we used the following property of the Brownian motion: for k=1,,d,

(F.1)E[(Bt,sk)r]={0,r is odd,r!2r/2(r/2)!(s-t)r/2,r is even.

G Proof of Lemma 4.4

By Lemma 4.2 and Lemma 4.3, we immediately have:

(Q0,T/n(m)Q(k-1)T/n,kT/n(m))kT/n,(k+1)T/n1P(k+1)T/n,Tf
C0kT/n,(k+1)T/n1P(k+1)T/n,TfC11nm+1f,

where C0,C1>0 are some constants independent of f:N, n1 and kn-1.

H Proof of Lemma 4.7

First, note that for r1, p=2e, e and ν (=ν(m)) there exists C(T)>0 (independent of n) such that

(H.1)ζ,kxr,pC(T)

for all xN and kn, which can be checked in the following way. Let

ΞT/nk:=i=1k𝒲(m)((i-1)T/n,iT/n,X¯(T/n)0,x((i-1)T/n),B(i-1)T/n,iT/n).

Note that

𝒲(m)((i-1)T/n,iT/n,X¯(T/n)0,x((i-1)T/n),B(i-1)T/n,iT/n)

is given by the sum of the form

v((i-1)T/n,X¯(T/n)0,x((i-1)T/n))(T/n)q1𝒫q2(B(i-1)T/n,iT/n),q1+q22,q1,q2,

where 𝒫r(B(i-1)T/n,iT/n) is polynomial of degree r. Let {t}t0 be the filtration generated by the Brownian motion. It holds that

E[|ΞT/nk|2e]=E[|ΞT/nk-1|2e|𝒲(m)((k-1)T/n,kT/n,X¯0,x((k-1)T/n),B(k-1)T/n,kT/n)|2e]
=E[|ΞT/nk-1|2eE[|𝒲(m)((k-1)T/n,kT/n,X¯0,x((k-1)T/n),B(k-1)T/n,kT/n)|2e|(k-1)T/n]]
=E[|ΞT/nk-1|2eE[|𝒲(m)((k-1)T/n,kT/n,,B(k-1)T/n,kT/n)|2e]|=X¯0,x((k-1)T/n)].

Also, there is C>0 such that

E[|𝒲(m)((k-1)T/n,kT/n,y,B(k-1)T/n,kT/n)|2]1+CT/n

for all yN, n1 and kn. Then, for all e,

(H.2)E[|ΞT/nk|2e]=(1+CT/n)k

for all kn.

Next, we show the upper bound of

E[DΞT/nkHp]=E[(0T|DsΞT/nk|2𝑑s)p/2].

By Hölder’s inequality, we see that

(0T|DsΞT/nk|2𝑑s)p/2c0T|DsΞT/nk|2e𝑑s

for some c>0. Then we will prove

E[0T(DsΞT/nk)2e𝑑s](1+CT/n)k

by induction. For k=1 it is obvious. We assume that there exists C>0 (independent of n) such that

(H.3)E[0T(DsΞT/nk-1)2e𝑑s](1+CT/n)k-1.

The chain rule of the Malliavin derivative gives

DsΞT/nk=DsΞT/nk-1𝟏s(k-1)T/n𝒲(m)((k-1)T/n,kT/n,X¯0,x((k-1)T/n),B(k-1)T/n,kT/n)
+ΞT/nk-1Ds𝒲(m)((k-1)T/n,kT/n,X¯0,x((k-1)T/n),B(k-1)T/n,kT/n)𝟏(k-1)T/nskT/n

and

D,s𝒲(m)((k-1)T/n,kT/n,X¯0,x((k-1)T/n),B(k-1)T/n,kT/n)
=r=1Nrv((k-1)T/n,X¯0,x((k-1)T/n)D,sX¯0,x,r((k-1)T/n)𝟏s(k-1)T/n(T/n)q1𝒫q2(B(k-1)T/n,kT/n)
   +v((k-1)T/n,X¯0,x((k-1)T/n)(T/n)q1l𝒫q2(B(k-1)T/n,kT/n)𝟏(k-1)T/nskT/n.

Here, l𝒫q2(B(k-1)T/n,kT/n) is a polynomial of order q2-1. By expanding (DsΞT/nk)2e and using the property of the polynomial of the Brownian motion (F.1), we have

E[0T(DsΞT/nk)2eds]=E[E[0T(DsΞT/nk)2eds|(k-1)T/n]]
E[G0,(0,(k-1)T/n)](1+O(T/n))+1lξ(e)E[Gl,(0,(k-1)T/n)]O(T/n)
(1+O(T/n))k

for some (k-1)T/n-measurable functions Gl,(0,(k-1)T/n), l=0,1,,ξ(e) with ξ(e), where the assumption (H.3) and the estimates

E[|ΞT/nk-1|2e](1+c1T/n)k-1

(by H.2) and

E[(0(k-1)T/n|DsX¯0,x((k-1)T/n)|2𝑑s)p]c2

(by [1, 2]) are applied (with Hölder’s inequality), where c1,c2>0 are constants which do not depend on n1. Therefore, we have

E[|ΞT/nk|p]+E[(0T|DsΞT/nk|2𝑑s)p/2]K1(1+K2T/n)kK1eK2T

for some K1,K2>0 independent of n1. We can proceed in the same manner to get (H.1).

We now focus on deriving formulas (4.1), (4.3) and the estimates (4.2), (4.4). In the following, we use the properties of the process X¯(T/n)0,x. Note that X¯(T/n)0,x can be written in the Itô process form as follows:

X¯(T/n)0,x(s)=x+0sV0(φ(u),X¯(T/n)0,x(φ(u)))𝑑u+i=1d0sVi(φ(u),X¯(T/n)0,x(φ(u)))𝑑Bui

for 0sT, where φ(s)=sup{kT/n:kT/ns}. In particular, under the uniform elliptic condition, X¯(T/n)0,x becomes an elliptic Itô process.

  1. For 0k[n2]-1, by the definition of Qt,s(m), we have

    (Q0,T/n(m)Q(k-1)T/n,kT/n(m))kT/n,(k+1)T/n2P(k+1)T/n,Tf(x)
    =E[kT/n,(k+1)T/n2P(k+1)T/n,Tf(X¯(T/n)0,x(kT/n))ΞT/nk]
    =(Tn)m+1νE[β()P(k+1)T/n,Tf(X¯(T/n)0,x((k+1)T/n)ζ,kx].

    In this case, we apply Proposition 4.5 with (H.1) to get

    (Q0,T/n(m)Q(k-1)T/n,kT/n(m))kT/n,(k+1)T/n2P(k+1)T/n,Tf
    C0(Tn)m+1ν|β()|P(k+1)T/n,TfC11nm+1f,

    where C0,C1>0 are constants independent of f:N, kn and n1.

  2. For [n2]kn-1, by applying the integration by parts of Kusuoka and Stroock [9], one has

    (Q0,T/n(m)Q(k-1)T/n,kT/n(m))kT/n,(k+1)T/n2P(k+1)T/n,Tf(x)
    =(Tn)m+1νE[P(k+1)T/n,Tf(X¯(T/n)0,x((k+1)T/n)Hβ()(X¯(T/n)0,x((k+1)T/n),ζ,kx)].

    We give the estimate for the weight

    Hβ()(X¯(T/n)0,x((k+1)T/n),ζ,kx).

    Let X¯(T/n)0,x(t), T2tT and G𝔻. Then, for p>1, there exist cp,α1,α2,q>0, k1,k21 such that

    Hβ()(X¯(T/n)0,x(t),G)pcpdet(σX¯(T/n)0,x(t))-1α1k1DX¯(T/n)0,x(t)|β()|,α2,Hk2G|β()|,q.

    We use the estimate in Proposition 4.6: for q>1,

    supT/2tT,xNdet(σX¯(T/n)0,x(t))-1qK1(T),

    And we use the following bound: for r, q>1,

    suptT,xNDX¯(T/n)0,x(t)r,q,HK2(T)

    by [1, 2] with (H.1). Therefore, we have that for p>1,

    sup[n/2]kn-1,xNHβ()(X¯(T/n)0,x((k+1)T/n),ζ,kx)pK3(T),

    where the constants Ki(T)>0, i=1,2,3, can be chosen uniformly in n1. Therefore, we have

    (Q0,T/n(m)Q(k-1)T/n,kT/n(m))kT/n,(k+1)T/n2P(k+1)T/n,TfC21nm+1P(k+1)T/n,TfC21nm+1f

    for some C2>0 independent of f:N, kn and n1.

References

[1] V. Bally and D. Talay, The law of the Euler scheme for stochastic differential equations. I. Convergence rate of the distribution function, Probab. Theory Related Fields 104 (1996), no. 1, 43–60. 10.1007/BF01303802Search in Google Scholar

[2] V. Bally and D. Talay, The law of the Euler scheme for stochastic differential equations. II. Convergence rate of the density, Monte Carlo Methods Appl. 2 (1996), no. 2, 93–128. 10.1515/mcma.1996.2.2.93Search in Google Scholar

[3] J.-M. Bismut, Large Deviations and the Malliavin Calculus, Progr. Math. 45, Birkhäuser, Boston, 1984. Search in Google Scholar

[4] L. G. Gyurkó and T. Lyons, Rough paths based numerical algorithms in computational finance, Mathematics in Finance, Contemp. Math. 515, American Mathematical Society, Providence (2010), 17–46. 10.1090/conm/515/10120Search in Google Scholar

[5] Y. Iguchi and T. Yamada, Operator splitting around Euler–Maruyama scheme and high order discretization of heat kernels, ESAIM Math. Model. Numer. Anal. (2020), 10.2139/ssrn.3510133. 10.2139/ssrn.3510133Search in Google Scholar

[6] N. Ikeda and S. Watanabe, Stochastic Differential Equations and Diffusion Processes, 2bd ed., North-Holland Math. Libr. 24, North-Holland Publishing, Amsterdam, 1981. Search in Google Scholar

[7] P. E. Kloeden and E. Platen, Numerical Solution of Stochastic Differential Equations, Appl. Math. (New York) 23, Springer, Berlin, 1992. 10.1007/978-3-662-12616-5Search in Google Scholar

[8] S. Kusuoka, Approximation of expectation of diffusion process and mathematical finance, Taniguchi Conference on Mathematics Nara ’98, Adv. Stud. Pure Math. 31, Mathematical Society of Japan, Tokyo (2001), 147–165. 10.2969/aspm/03110147Search in Google Scholar

[9] S. Kusuoka and D. Stroock, Applications of the Malliavin calculus. I, Stochastic Analysis (Katata/Kyoto 1982), North-Holland Math. Library 32, North-Holland, Amsterdam (1984), 271–306. 10.1016/S0924-6509(08)70397-0Search in Google Scholar

[10] T. Lyons and N. Victoir, Cubature on Wiener space, Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 460 (2004), 169–198. 10.1098/rspa.2003.1239Search in Google Scholar

[11] P. Malliavin, Stochastic Analysis, Grundlehren Math. Wiss. 313, Springer, Berlin, 1997. 10.1007/978-3-642-15074-6Search in Google Scholar

[12] P. Malliavin and A. Thalmaier, Stochastic Calculus of Variations in Mathematical Finance, Springer Finance, Springer, Berlin, 2006. Search in Google Scholar

[13] R. Naito and T. Yamada, A third-order weak approximation of multidimensional Itô stochastic differential equations, Monte Carlo Methods Appl. 25 (2019), no. 2, 97–120. 10.1515/mcma-2019-2036Search in Google Scholar

[14] D. Nualart, The Malliavin Calculus and Related Topics, 2nd ed., Probab. Appl., Springer, Berlin, 2006. Search in Google Scholar

[15] A. Takahashi, Asymptotic expansion approach in finance, Large Deviations and Asymptotic Methods in Finance, Springer Proc. Math. Stat. 110, Springer, Cham (2015), 345–411. 10.1007/978-3-319-11605-1_13Search in Google Scholar

[16] A. Takahashi and T. Yamada, An asymptotic expansion with push-down of Malliavin weights, SIAM J. Financial Math. 3 (2012), no. 1, 95–136. 10.1137/100807624Search in Google Scholar

[17] A. Takahashi and T. Yamada, A weak approximation with asymptotic expansion and multidimensional Malliavin weights, Ann. Appl. Probab. 26 (2016), no. 2, 818–856. 10.1214/15-AAP1105Search in Google Scholar

[18] S. Watanabe, Analysis of Wiener functionals (Malliavin calculus) and its applications to heat kernels, Ann. Probab. 15 (1987), no. 1, 1–39. 10.1214/aop/1176992255Search in Google Scholar

[19] T. Yamada, An arbitrary high order weak approximation of SDE and Malliavin Monte Carlo: Analysis of probability distribution functions, SIAM J. Numer. Anal. 57 (2019), no. 2, 563–591. 10.1137/17M114412XSearch in Google Scholar

Received: 2020-08-10
Revised: 2021-01-20
Accepted: 2021-02-02
Published Online: 2021-02-20
Published in Print: 2021-06-01

© 2021 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 27.4.2026 from https://www.degruyterbrill.com/document/doi/10.1515/mcma-2021-2085/html?lang=en
Scroll to top button