A = (1, 1)
B = (5, 2)
C = (10, 5)
xA: x(A)
xB: x(B)
xC: x(C)
yA: y(A)
yB: y(B)
yC: y(C)
a: (xA * (yC - yB) - xB * yC + xC * yB + (xB - xC) * yA) / (xA * (xC^(2.0) - xB^(2.0)) - xB * xC^(2.0) + xB^(2.0) * xC + xA^(2.0) * (xB - xC))
b: -1.0 * (xA^(2.0) * (yC - yB) - xB^(2.0) * yC + xC^(2.0) * yB + (xB^(2.0) - xC^(2.0)) * yA) / (xA * (xC^(2.0) - xB^(2.0)) - xB * xC^(2.0) + xB^(2.0) * xC + xA^(2.0) * (xB - xC))
c: (xA * (xC^(2.0) * yB - xB^(2.0) * yC) + xA^(2.0) * (xB * yC - xC * yB) + (xB^(2.0) * xC - xB * xC^(2.0)) * yA) / (xA * (xC^(2.0) - xB^(2.0)) - xB * xC^(2.0) + xB^(2.0) * xC + xA^(2.0) * (xB - xC))
f(x) = a * x^(2.0) + b * x + c
g = Function(f(x), 0.0, x(A))
h = Function(f(x), x(A), x(C))
p = Function(f(x), x(C), 100.0)
Text1: "f(x) = " + (Simplify[f(x)])