Skip to main content
11 events
when toggle format what by license comment
Aug 24, 2023 at 0:36 vote accept physicsu83
Sep 13, 2019 at 5:20 history edited Αλέξανδρος Ζεγγ CC BY-SA 4.0
added 1 character in body
Sep 13, 2019 at 5:15 answer added Bob Hanlon timeline score: 2
Sep 12, 2019 at 18:27 comment added bbgodfrey You are changing the question. Please open a new question that includes both the code above and the boundary condition at large r that you are trying to match. My guess is that you are trying to match a separatrix at large `r', which typically is difficult but possible numerically.
Sep 12, 2019 at 17:42 comment added physicsu83 Actually, my differential equation is a boundary value problem. With x[0.0001]=Pi and x[infinity]=0. So I used the shooting method and converted it into an initial value problem with x[0.0001]=Pi and x'[0.0001]=-(some guessed value). But my gussed value gives oscillation at large distance. But I need x[infinity]=0 not oscillation.
Sep 12, 2019 at 17:35 comment added march I'm sorry but I don't understand your question. The solution is the solution. Are you talking about modifying your model so that it does that or defining a new function that is the solution above until $r=20$ and then zero afterward? (And what do you mean by zero? Do you mean it's at the origin after that or that the vertical component is zero?)
Sep 12, 2019 at 17:20 history edited physicsu83 CC BY-SA 4.0
added 7 characters in body
Sep 12, 2019 at 17:12 comment added physicsu83 Thank you very much. It worked. Do you have any idea that how I can get rid of oscillation at r>20? I need solution to be zero as like it is before r=20.
Sep 12, 2019 at 17:03 comment added march ParametricPlot[{x[r] /. First@sols4, D[x[r] /. First@sols4, r]} // Evaluate, {r, 0.00001, 100}]
Sep 12, 2019 at 16:45 review First posts
Sep 12, 2019 at 18:03
Sep 12, 2019 at 16:44 history asked physicsu83 CC BY-SA 4.0