Throughout, ∥⋅∥2 and ⟨⋅,⋅⟩2 are the mean-square norm and inner product of Square-Integrable Random Variables and the Mean-Square Inner Product.
Claim 1. By the definitions of the dot product and the matrix-vector product, applied pointwise on Ω,
Y⋅(MY)=i=1∑pj=1∑pMijYiYj.
For each pair (i,j) the product YiYj is integrable with E[∣YiYj∣]≤∥Yi∥2∥Yj∥2, by Cauchy-Schwarz and Triangle Inequalities for the Mean-Square Norm applied to ∣Yi∣ and ∣Yj∣, which are square-integrable by the closure properties of Square-Integrable Random Variables and the Mean-Square Inner Product (their squares agree with those of Yi,Yj). Hence the finite sum is integrable and, by linearity of the integral (Linearity and Monotonicity of the Lebesgue Integral),
E[Y⋅(MY)]=i,j∑MijE[YiYj].
By the definition of the covariance, E[YiYj]=Cov(Yi,Yj)+E[Yi]E[Yj]=Cij+μiμj. Therefore
E[Y⋅(MY)]=i,j∑MijCij+i,j∑Mijμiμj=tr(M⊤C)+μ⋅(Mμ),
using claim 4 of Basic Properties of the Trace for the first sum and the definitions of dot product and matrix-vector product for the second. Finally, C is symmetric (Cov(Yi,Yj)=Cov(Yj,Yi) by the symmetry of its defining formula), so claims 2-3 of Basic Properties of the Trace give tr(M⊤C)=tr((M⊤C)⊤)=tr(C⊤M)=tr(CM)=tr(MC), with (M⊤C)⊤=C⊤M by claim 3 of Componentwise Estimates, Transpose Identities, and Indefinite Riemann Integrals.
Claim 2. As in claim 1, for each t,
E[Yt⋅(M(t)Zt)]=i=1∑pj=1∑qMij(t)E[YtiZtj],
all terms defined and finite. It therefore suffices, by the continuity of sums and products of continuous real functions (Sums and Products of Continuous Real-Valued Functions), to show that each function t↦E[YtiZtj]=⟨Yti,Ztj⟩2 is continuous on [a,b]. Fix s,t∈[a,b]. Adding and subtracting E[YsiZtj] and applying Cauchy-Schwarz and Triangle Inequalities for the Mean-Square Norm twice,
E[YtiZtj]−E[YsiZsj]≤∥Yti−Ysi∥2∥Ztj∥2+∥Ysi∥2∥Ztj−Zsj∥2.
By claim 4 of Basic Properties of the Mean-Square Riemann Integral, the functions t↦∥Yti∥2 and t↦∥Ztj∥2 are continuous on the compact interval [a,b], hence bounded there by Extreme Value Theorem on a Compact Interval; and ∥Yti−Ysi∥2→0, ∥Ztj−Zsj∥2→0 as t→s by mean-square continuity. Hence the right-hand side tends to 0 as t→s, which is the asserted continuity at s. □