Function satisfying the relation $f(x+y)=f(x)+f(y)-(e^{-x}-1)(e^{-y}-1)+1$ Let f be the differentiable function satisfying the relation  $f\left( {x + y} \right) = f\left( x \right) + f\left( y \right) - \left( {{e^{ - x}} - 1} \right)\left( {{e^{ - y}} - 1} \right) + 1$; $\forall x,y \in R$ and $\mathop {\lim }\limits_{h \to 0} \frac{{f'\left( {1 + h} \right) + f\left( h \right) - {e^{ - 1}}}}{h}$ exist. The value of $\int\limits_0^1 {f\left( x \right)dx}  = \_\_\_\_\_\_\_$.
My approach is as follow
$f\left( {x + y} \right) = f\left( x \right) + f\left( y \right) - \left( {{e^{ - x}} - 1} \right)\left( {{e^{ - y}} - 1} \right) + 1$
$\Rightarrow f\left( {x + y} \right) - f\left( x \right) = f\left( y \right) - \left( {{e^{ - x}} - 1} \right)\left( {{e^{ - y}} - 1} \right) + 1$
$ \Rightarrow \mathop {\lim }\limits_{h \to 0} \frac{{f\left( {x + y} \right) - f\left( x \right)}}{h} = \mathop {\lim }\limits_{h \to 0} \frac{{f\left( y \right) - \left( {{e^{ - x}} - 1} \right)\left( {{e^{ - y}} - 1} \right) + 1}}{h}$
$ \Rightarrow \mathop {\lim }\limits_{h \to 0} \frac{{f\left( {x + y} \right) - f\left( x \right)}}{h} = \mathop {\lim }\limits_{h \to 0} \frac{{f\left( y \right) - \left( {{e^{ - x}} - 1} \right)\left( {{e^{ - y}} - 1} \right) + 1}}{h}$
y=h
$ \Rightarrow f'\left( x \right) = \mathop {\lim }\limits_{h \to 0} \frac{{f\left( {x + h} \right) - f\left( x \right)}}{h} = \mathop {\lim }\limits_{h \to 0} \frac{{f\left( h \right) - \left( {{e^{ - x}} - 1} \right)\left( {{e^{ - h}} - 1} \right) + 1}}{h}$
How will I proceed from here

If you define $g(x) = f(x) + e^{-x}$, then the equation becomes $g(x + y) = g(x) + g(y)$.
It is well-known that all continuous solutions to that equation are linear, i.e. there exists $c$ such that $f(x) = cx$.
The rest should be easy.
