PedagoSpline

Morspline(<Point A>,<Point B>,<Poignee A>,<Poignee B>)

f(x) = If[x(A) ≤ x ≤ x(B), (x^(3) * ((y(A2) - y(A)) / (x(A2) - x(A)) / ((x(A))^(2) + (x(B))^(2) - ((2 * x(A)) * x(B))) + (y(B2) - y(B)) / (x(B2) - x(B)) / ((x(A))^(2) + (x(B))^(2) - ((2 * x(A)) * x(B))) - (2 * y(A) / ((x(A))^(3) - (x(B))^(3) + ((3 * x(A)) * (x(B))^(2)) - ((3 * (x(A))^(2)) * x(B)))) + (2 * y(B) / ((x(A))^(3) - (x(B))^(3) + ((3 * x(A)) * (x(B))^(2)) - ((3 * (x(A))^(2)) * x(B)))))) + (x^(2) * (((y(A2) - y(A)) / (x(A2) - x(A)) * ((-x(A)) - (2 * x(B))) / ((x(A))^(2) + (x(B))^(2) - ((2 * x(A)) * x(B)))) + ((y(B2) - y(B)) / (x(B2) - x(B)) * ((-2 * x(A)) - x(B)) / ((x(A))^(2) + (x(B))^(2) - ((2 * x(A)) * x(B)))) + (y(A) * ((3 * x(A)) + (3 * x(B))) / ((x(A))^(3) - (x(B))^(3) + ((3 * x(A)) * (x(B))^(2)) - ((3 * (x(A))^(2)) * x(B)))) + (y(B) * ((-3 * x(A)) - (3 * x(B))) / ((x(A))^(3) - (x(B))^(3) + ((3 * x(A)) * (x(B))^(2)) - ((3 * (x(A))^(2)) * x(B)))))) + (x * (((y(A2) - y(A)) / (x(A2) - x(A)) * ((x(B))^(2) + ((2 * x(A)) * x(B))) / ((x(A))^(2) + (x(B))^(2) - ((2 * x(A)) * x(B)))) + ((y(B2) - y(B)) / (x(B2) - x(B)) * ((x(A))^(2) + ((2 * x(A)) * x(B))) / ((x(A))^(2) + (x(B))^(2) - ((2 * x(A)) * x(B)))) - (((6 * x(A)) * x(B)) * y(A) / ((x(A))^(3) - (x(B))^(3) + ((3 * x(A)) * (x(B))^(2)) - ((3 * (x(A))^(2)) * x(B)))) + (((6 * x(A)) * x(B)) * y(B) / ((x(A))^(3) - (x(B))^(3) + ((3 * x(A)) * (x(B))^(2)) - ((3 * (x(A))^(2)) * x(B)))))) + ((((-(y(A2) - y(A)))) / (x(A2) - x(A)) * x(A)) * (x(B))^(2) / ((x(A))^(2) + (x(B))^(2) - ((2 * x(A)) * x(B)))) - (((y(B2) - y(B)) / (x(B2) - x(B)) * (x(A))^(2)) * x(B) / ((x(A))^(2) + (x(B))^(2) - ((2 * x(A)) * x(B)))) + (y(A) * ((-(x(B))^(3)) + ((3 * x(A)) * (x(B))^(2))) / ((x(A))^(3) - (x(B))^(3) + ((3 * x(A)) * (x(B))^(2)) - ((3 * (x(A))^(2)) * x(B)))) + (y(B) * ((x(A))^(3) - ((3 * (x(A))^(2)) * x(B))) / ((x(A))^(3) - (x(B))^(3) + ((3 * x(A)) * (x(B))^(2)) - ((3 * (x(A))^(2)) * x(B))))]