Timeline for Overflow with NDSolve
Current License: CC BY-SA 4.0
3 events
| when toggle format | what | by | license | comment | |
|---|---|---|---|---|---|
| Jun 7, 2020 at 12:36 | comment | added | Michael E2 |
@xzczd For reasons I can't fathom, the shooting method is varying x near x == 0 instead of integrating to x == rin. I'd think any programmer would understand that stepping outside the specified domain is inviting trouble. You can fix it by using Sqrt[x + $MinMachineNumber] instead of Sqrt[x], which won't make a difference if x > 10^-292, since x + $MinMachineNumber rounds to exact x in machine precision. It's a pretty ugly fix though.
|
|
| Jun 7, 2020 at 4:16 | comment | added | xzczd♦ |
It's a bit surprising that WorkingPrecision -> 20 seems to be necessary to make Shooting method work.
|
|
| Jun 7, 2020 at 3:59 | history | answered | Michael E2 | CC BY-SA 4.0 |