Skip to main content
8 events
when toggle format what by license comment
Nov 16, 2017 at 4:29 history edited Carl Woll CC BY-SA 3.0
grammar
Nov 16, 2017 at 2:28 comment added bbgodfrey @CarlWoll What version of Mathematica are you using? When I run your code with 11.2 on Windows 10, the blue curve does not appear at all, although the InterpolationFunction for the blue curve (in the second quadrant) certainly exists.
Nov 16, 2017 at 2:25 comment added bbgodfrey @EmilioPisanty The first InterpolationFunction must be extrapolating in some way to obtain the blue curve in the third and fourth quadrants, because (as the error message notes) s stops at 0.9165575855524158.
Nov 15, 2017 at 11:32 comment added Emilio Pisanty @J.M. Frankly, I don't understand why it's encountering any hiccups at that point - the RHS maps to 1, and that's smack in the safest zone of the principal branch of Power[]. The kink in Carl's plot seems to be extraneous, though - if I run that NDSolve, I get the same mxst error, and looking at the internals of the first InterpolatingFunction, it doesn't have any data at all in either the third or the fourth quadrants; the solution just stops at $s\approx 0.916\leftrightarrow y=0$.
Nov 15, 2017 at 10:51 comment added J. M.'s missing motivation @Emilio, that might be related to the need to take the principal branch of Power[]. I'm guessing the solution you actually want corresponds to a different branch.
Nov 15, 2017 at 7:30 comment added Emilio Pisanty This has an equivalent problem with the kink when the curve reaches the negative real axis. This one should be smooth and stay in the third quadrant.
Nov 15, 2017 at 1:54 comment added J. M.'s missing motivation Was going to suggest this myself... (+1)
Nov 15, 2017 at 1:45 history answered Carl Woll CC BY-SA 3.0