Given $a,b\in\mathbb{R}$ with $a < b$ and defining $F(z):=\int_0^z f(s) \, ds$ with $z \in \mathbb{R}$, how can one establish that $$F(a+b)=F(a)+f(a)b+ b^2\int_0^1 (1-s)f'(a+sb) \, ds,$$ which is part of a larger proof on page 508 of PDE Evans (2nd edition)?
Should I perform a change of variable to the second term on the RHS?
First, $$ F(a+b)-F(a)=\int_a^{a+b}f(t)\,dt. $$ Then, integrate by parts, $$ F(a+b)-F(a)=[(t-a-b)f(t)]_{a}^{a+b}-\int_a^{a+b}(t-a-b) f'(t)\,dt. $$ Now, let $t=a+sb$, $$ \begin{aligned} F(a+b)-F(a)&=[(t-a-b)f(t)]_{a}^{a+b}+b^2\int_0^1 (1-s)f'(a+sb)\,ds\\ &=f(a)b+b^2\int_0^1 (1-s)f'(a+sb)\,ds. \end{aligned} $$