Let $\Phi$ and $\varphi$ denote the cumulative density function and the density function of a standard normal random variable. I'd like to calculate
$$I=\int_{-\infty}^{\infty}e^{cx}\varphi(x)\Phi(a+bx)dx.$$
This is similar to Expected Value of Normal Random Variable times its CDF where the term $e^{ax}$ is replaced with a linear term. Tried to apply similar tricks here but didn't work out.