El truco aquí es darse cuenta de que xi y x∗i son en este caso constantes.
cov(^yi,^y∗i)=cov(^β0+^β1xi,^β0+^β1x∗i)=cov(^β0,^β0)+cov(^β0,^β1x∗i)+cov(^β1xi,^β0)+cov(^β1xi,^β1x∗i)=var(^β0)+x∗icov(^β0,^β1)+xicov(^β1,^β0)+xix∗ivar(^β1)=var(^β0)+cov(^β0,^β1)(xi+x∗i)+xix∗ivar(^β1) .
Ahora el resto debería ser fácil ya que el var(^β0),var(^β1),cov(^β0,^β1) son bien conocidos. Vea aquí los valores exactos de varianza/covarianza: Cómo derivar la matriz de varianza-covarianza de los coeficientes en la regresión lineal
Editar:
En respuesta a tu comentario, realmente es una cuestión de gusto hasta dónde lleves esto a partir de ahora. Yo procedería de la siguiente manera:
var(^β0)+cov(^β0,^β1)(xi+x∗i)+xix∗ivar(^β1)=σ2∑x2in∑(xi−¯X)2+(xi+x∗i)−¯Xσ2∑(xi−¯X)2+xix∗iσ2∑(xi−¯X)2=σ2∑(xi−¯X)2(∑x2in−¯X(xi+x∗i)+xix∗i)