a = 0.125
b = 0.1
c = -1
d = 0
e = 1
f(x) = a * x^(4.0) + b * x^(3.0) + c * x^(2.0) + d * x + e
g = Derivative(f(x))
h = Derivative(g(x))
N_1, N_2, N_3, N_4 = Intersect(f, xAxis)
A, B, C = Intersect(g, xAxis)
i = OrthogonalLine(A, xAxis)
j = OrthogonalLine(B, xAxis)
E_1 = Intersect(f, i)
E = Intersect(j, yAxis)
F = Intersect(j, yAxis)
E_2 = Intersect(f, j)
k = OrthogonalLine(C, xAxis)
E_3 = Intersect(f, k)
D, G = Intersect(h, xAxis)
l = OrthogonalLine(D, xAxis)
m = OrthogonalLine(G, xAxis)
W_1 = Intersect(f, l)
I = Intersect(h, m)
W_2 = Intersect(f, m)
text1: "f(x) = " + f