I attempt to solve the brachistochrone problem numerically. I am using a direct method which considers the curve ##y(x)## as a Lagrange polynomial evaluated at fixed nodes ##x_i##, and the time functional as a multivariate function of the ##y_i##. The classical statement of the problem requires...
Whenever I try to reduce, invert, or factorize matrices of size 10000x10000 my whole system suddenly restarts without any error message. This happens on both julia and MATLAB when I run the command A=rand(10000,10000)^-1;
I suspect the problem is due to the linear algebra libraries (I have...