Given a probability distribution with joint CDF:
$F(x,y)= 1 - e^{-x}-e^{-y}+e^{-(x+y + \alpha x y)} \quad \text{for} \quad x >0 ,\quad y>0, \quad \alpha \in [0,1]$
We can find the corresponding PDF through the differentiation of both variables and after simplifying we obtain:
$f(x,y) = -\alpha e^{-x-y-\alpha x y} + e^{-x-y-\alpha xy}(-1-\alpha x)(-1-\alpha y)$
The next step would be to obtain the marginal PDFs from $f(x,y)$ , however as the function is clearly non-separable into $f(x,y) = g(x)h(y)$ along with the variables clearly being dependent, this seems to be a quite nontrivial. So far I have concluded that given the form of the PDF, the marginals must be distributed exponentially for some parameter $\lambda(\alpha)$, however I am unsure of how to proceed to show this to be the case without having to perform excessive computations.