Ein Weg geht mit dem erweiterten Mittelwertsatz (siehe Abschnitt "Erweiterter Mittelwertsatz ...").
Ich setze mal aus Bequemlichkeit
x1=a,x2=b
Die linke Seite deiner Gleichung lässt sich so umformen, dass sie zum erweiterten MWS passt:
a−b1(af(b)−bf(a))=aba−bbf(b)−af(a)=...
...=b1−a1bf(b)−af(a)(1)
Jetzt wenden wir den verallgemeinerten MWS an, mit den Funktionen:
g(x)=xf(x) und h(x)=x1
Dazu brauchen wir
g′(x)=x2f′(x)x−f(x)
h′(x)=−x21
Damit haben wir
(1)=h′(ξ)g′(ξ)=−ξ21ξ2f′(ξ)ξ−f(ξ)=f(ξ)−f′(ξ)ξ
p.s.: Wegen ab>0 enthält das betrachtete Intervall nicht die Null.