d = 1
c = -1.5
b = 0
a = 0.125
f(x) = a * x^(3.0) + b * x^(2.0) + c * x + d
g = Derivative(f(x))
A, B = Intersect(g, xAxis)
e = OrthogonalLine(A, xAxis)
h = OrthogonalLine(B, xAxis)
E_1 = Intersect(f, e)
E_2 = Intersect(f, h)
k = Derivative(g(x))
C = Intersect(k, xAxis)
i = OrthogonalLine(C, xAxis)
W = Intersect(f, i)
N_1, N_2, N_3 = Intersect(f, xAxis)
text1: "f(x) = " + f