I'm trying to solve the brachistochrone problem by numerical optimization. I started with the minimal surface problem to get a foot in the door: I discretize the integral as a sum of constributions that depend on the function values at specific points, then use autograd libraries for optimizing a non-linear scalar loss function