Esta solución no es ni mucho menos la más sencilla; utilizaremos curvas de Bézier (véase wiki ).
Supongamos que $(d-c)>\gamma (b-a)$ .
Construiremos una curva cúbica de Bézier con $P_0= (a,c)$ - tu punto de partida, $P_3 = (b,d)$ - que terminas el punto, y $P_1= (b,c+\gamma(b-a))$ - la dirección de la derivada en el punto inicial, de forma similar, $P_2=(a, d-\gamma(b-a))$ - la dirección en el punto de llegada.
La curva de Bézier se escribe entonces en el plano $\Bbb R^2$ :
$$x(t) =a\cdot (1-t)^3 + b\cdot 3t(1-t)^2+ a\cdot 3t^2(1-t)+b\cdot t^3,$$ $$y(t) =c\cdot (1-t)^3 + (c+\gamma(b-a))\cdot 3t(1-t)^2+ (d-\gamma(b-a))\cdot 3t^2(1-t)+d\cdot t^3.$$
Es fácil demostrar que $t\to x(t)$ y $t\to t$ son estrictamente monótonas (basta con estudiar la derivada), por lo tanto, la aplicación inversa $x\to t(x)$ está bien definida (incluso más - es $C^1$ ).
Veamos la derivada de $y$ cuando $x=a$ (es decir $t=0$ ):
$$\frac{dy}{dt}\big|_{t=0} = \frac{dy}{dx}\big|_{x=a}\cdot\frac{dx}{dt}\big|_{t=0}. $$ Obtenemos $$ -3c+3(c+\gamma(b-a)) = \frac{dy}{dx}\big|_{x=a} \cdot 3(b-a), $$ que da $$\frac{dy}{dx}\big|_{x=a}=\gamma.$$ La prueba para $\frac{dy}{dx}\big|_{x=b}=\gamma $ se hace de forma similar.
Por monotonía estricta de $x $ y $y$ concluimos que la aplicación $y(x)$ es de hecho $[a,b]\to[c,d]$ y satisface las condiciones que usted impuso.
El caso $(d-c)<\gamma (b-a)$ requiere vectores ligeramente modificados $P_1$ y $P_2$ :
$$P_1=(a+(d-c)/\gamma,d),\quad P_2=(b-(d-c)/\gamma,c).$$
El caso $\gamma = \frac{d-c}{b-a}$ es trivial.
La intuición detrás de la elección de $P_1$ : queremos una derivada $\gamma$ Así pues, partimos de $(a,c)$ y vaya con la misma derivada constante $\gamma$ hasta llegar al límite de nuestra caja $[a,b]\times [c,d]$ . Del mismo modo, para $P_2$ simplemente vamos en otra dirección.