b = 0 (slider from -5 to 5, step 0.1)
a = 0 (slider from -5 to 5, step 0.1)
d: (3 * (a - b))
c: a + b
cor(t) = ((5 * ((c * t) + d)) * (t^(2) - 1))
x0 = -1
B0 = Point(Line[(x0, 0), (x0, 1)])
x1 = -0.8
B1 = Point(Line[(x1, 0), (x1, 1)])
x2 = -0.6
B2 = Point(Line[(x2, 0), (x2, 1)])
x3 = -0.4 (slider from -5 to 5, step 1)
B3 = Point(Line[(x3, 0), (x3, 1)])
x4 = -0.2 (slider from -5 to 5, step 1)
B4 = Point(Line[(x4, 0), (x4, 1)])
x5 = 0
B5 = Point(Line[(x5, 0), (x5, 1)])
x6 = 0.2
B6 = Point(Line[(x6, 0), (x6, 1)])
x7 = 0.4
B7 = Point(Line[(x7, 0), (x7, 1)])
x8 = 0.6
B8 = Point(Line[(x8, 0), (x8, 1)])
x9 = 0.8
B9 = Point(Line[(x9, 0), (x9, 1)])
x10 = 1
B10 = Point(Line[(x10, 0), (x10, 1)])
M = 1000
ψ = Spline({(x(B0), (y(B0) + cor(x(B0))) / M), (x(B1), (y(B1) + cor(x(B1))) / M), (x(B2), (y(B2) + cor(x(B2))) / M), (x(B3), (y(B3) + cor(x(B3))) / M), (x(B4), (y(B4) + cor(x(B4))) / M), (x(B5), (y(B5) + cor(x(B5))) / M), (x(B6), (y(B6) + cor(x(B6))) / M), (x(B7), (y(B7) + cor(x(B7))) / M), (x(B8), (y(B8) + cor(x(B8))) / M), (x(B9), (y(B9) + cor(x(B9))) / M), (x(B10), (y(B10) + cor(x(B10))) / M)}, 3)
f = CurveCartesian(x(ψ(t)), (M * y(ψ(t))) - cor(x(ψ(t))), t, 0, 1)