-- https://mathworld.wolfram.com/BrachistochroneProblem.html x = theta - sin(theta) y = 1 - cos(theta) f = (d(x,theta)^2 + d(y,theta)^2) / (2 g y) f = simplify(f) f thetaA = -pi thetaB = 0 A = eval((x,y), theta, thetaA) A B = eval((x,y), theta, thetaB) B "Time to go from A to B" t = defint(sqrt(f), theta, thetaA, thetaB) t Run