In a recent post, I asked for help deriving a computable formula for $\text{Var}[\max(V-K,0)]$ based on the approach on p. 262 of ths book. $V$ is a lognormally distributed random variable and $K$ is a constant. And help did I get.
Now I need to extend this to the more general covariance case. I.e. I need a computable expression for
$$\text{Cov}[\max(V_1-K_1,0), \max(V_2-K_2, 0)]$$
where as before $V_1, V_2$ are lognormally distributed while $K_1, K_2$ are constants.
As before, I know the key is to find an expression for
$$E[\max(V_1-K_1,0) \max(V_2-K_2, 0)]$$
But I am stuck.