Delphi ist eigentlich kein Taschenrechner, der sowas einfach löst.
Wenn du willst, das Delphi so etwas für dich löst, müsstest du dafür ein eigenes Programm schreiben (oder finden).
Bei so einem linearen Gleichungssystem sollte dir das Gaußsche Eliminationsverfahren helfen.
Die Gleichungen könntest du dann in eine Matrix umschreiben und dann mit dem Gauß auflösen lassen.
Code:
a² a q -f(a)
xm² xm q -f(xm)
b² b q -f(b)