How to plot an unstable attractor?Why does DSolve return two solutions for my ODE?Second Order ODE (Bessel Function) with a dependent variable in the BCDSolve doesn't give all the solutionsUsing DSolve for a coupled differential equationDSolve not satisfying initial conditionshow to To specify initial conditions for a system of ode?DAE with NDSolve -monitor numerical noiseDSolve - Unable to obtain plot of solution - 2nd order ODEPartial differential equation heat/diffusion equation 3dCannot solve ODE question with Initial Value
Is this story about US tax office reasonable?
Infinite Sequence based on Simple Rule
How many chess players are over 2500 Elo?
Why do Russians call their women expensive ("дорогая")?
Mother abusing my finances
What is the object moving across the ceiling in this stock footage?
Placing bypass capacitors after VCC reaches the IC
Can you heal a summoned creature?
Under what law can the U.S. arrest International Criminal Court (ICC) judges over war crimes probe?
1960s sci-fi novella with a character who is treated as invisible by being ignored
Canon 70D often overexposing or underexposing shots
Are there situations when self-assignment is useful?
What is the largest (size) solid object ever dropped from an airplane to impact the ground in freefall?
Windows 10 Programs start without visual Interface
What are the benefits of cryosleep?
At what point in European history could a government build a printing press given a basic description?
Is there any use case for the bottom type as a function parameter type?
What is the most important source of natural gas? coal, oil or other?
Ticket sales for Queen at the Live Aid
Integrating an absolute function using Mathematica
Were pen cap holes designed to prevent death by suffocation if swallowed?
How bitcoin nodes update UTXO set when their latests blocks are replaced?
Is this resistor leaking? If so, is it a concern?
Looking for a soft substance that doesn't dissolve underwater
How to plot an unstable attractor?
Why does DSolve return two solutions for my ODE?Second Order ODE (Bessel Function) with a dependent variable in the BCDSolve doesn't give all the solutionsUsing DSolve for a coupled differential equationDSolve not satisfying initial conditionshow to To specify initial conditions for a system of ode?DAE with NDSolve -monitor numerical noiseDSolve - Unable to obtain plot of solution - 2nd order ODEPartial differential equation heat/diffusion equation 3dCannot solve ODE question with Initial Value
$begingroup$
I'm trying to solve and plot the following in Mathematica:
eqns = x'[t] ==
x[t] - y[t] -
x[t] (x[t]^2 + y[t]^2) + (x[t] y[t])/Sqrt[x[t]^2 + y[t]^2],
y'[t] ==
x[t] + y[t] - y[t] (x[t]^2 + y[t]^2) -
x[t]^2/Sqrt[x[t]^2 + y[t]^2];
DSolve[eqns, x, y, t]
This is supposed to be an example of unstable attractor ODE. However, execution never ends and I don't manage to see the solution of the equation.
plotting differential-equations
edited 1 hour ago
user64494
3,88711323
asked 8 hours ago
JavierJavier
1305
$endgroup$
add a comment |
$begingroup$
I'm trying to solve and plot the following in Mathematica:
eqns = x'[t] ==
x[t] - y[t] -
x[t] (x[t]^2 + y[t]^2) + (x[t] y[t])/Sqrt[x[t]^2 + y[t]^2],
y'[t] ==
x[t] + y[t] - y[t] (x[t]^2 + y[t]^2) -
x[t]^2/Sqrt[x[t]^2 + y[t]^2];
DSolve[eqns, x, y, t]
This is supposed to be an example of unstable attractor ODE. However, execution never ends and I don't manage to see the solution of the equation.
plotting differential-equations
edited 1 hour ago
user64494
3,88711323
asked 8 hours ago
JavierJavier
1305
$endgroup$
$begingroup$
Try usingNDSolveinstead
$endgroup$
– b3m2a1
8 hours ago
add a comment |
$begingroup$
I'm trying to solve and plot the following in Mathematica:
eqns = x'[t] ==
x[t] - y[t] -
x[t] (x[t]^2 + y[t]^2) + (x[t] y[t])/Sqrt[x[t]^2 + y[t]^2],
y'[t] ==
x[t] + y[t] - y[t] (x[t]^2 + y[t]^2) -
x[t]^2/Sqrt[x[t]^2 + y[t]^2];
DSolve[eqns, x, y, t]
This is supposed to be an example of unstable attractor ODE. However, execution never ends and I don't manage to see the solution of the equation.
plotting differential-equations
edited 1 hour ago
user64494
3,88711323
asked 8 hours ago
JavierJavier
1305
$endgroup$
I'm trying to solve and plot the following in Mathematica:
eqns = x'[t] ==
x[t] - y[t] -
x[t] (x[t]^2 + y[t]^2) + (x[t] y[t])/Sqrt[x[t]^2 + y[t]^2],
y'[t] ==
x[t] + y[t] - y[t] (x[t]^2 + y[t]^2) -
x[t]^2/Sqrt[x[t]^2 + y[t]^2];
DSolve[eqns, x, y, t]
This is supposed to be an example of unstable attractor ODE. However, execution never ends and I don't manage to see the solution of the equation.
plotting differential-equations
plotting differential-equations
edited 1 hour ago
user64494
3,88711323
asked 8 hours ago
JavierJavier
1305
edited 1 hour ago
user64494
3,88711323
asked 8 hours ago
JavierJavier
1305
edited 1 hour ago
user64494
3,88711323
edited 1 hour ago
user64494
3,88711323
edited 1 hour ago
user64494
3,88711323
3,88711323
asked 8 hours ago
JavierJavier
1305
asked 8 hours ago
JavierJavier
1305
asked 8 hours ago
JavierJavier
1305
1305
$begingroup$
Try usingNDSolveinstead
$endgroup$
– b3m2a1
8 hours ago
add a comment |
$begingroup$
Try usingNDSolveinstead
$endgroup$
– b3m2a1
8 hours ago
$begingroup$
Try using
NDSolve instead$endgroup$
– b3m2a1
8 hours ago
$begingroup$
Try using
NDSolve instead$endgroup$
– b3m2a1
8 hours ago
add a comment |
2 Answers
2
active
oldest
votes
$begingroup$
To visualize a 2D system, I would start with StreamPlot:
vf = x', y' /. First@Solve[eqns /. f_[t] :> f, x', y']; (* strip the args *)
StreamPlot[vf, x, -2, 2, y, -2, 2]
You can use StreamPoints to highlight the structure and Epilog to mark the attractor at $(1,0)$:
ics = Cos[1/5], Sin[1/5], Red,
0.5, 0, Magenta, 1.5, 0., Magenta;
StreamPlot[vf, x, -2, 2, y, -2, 2,
StreamPoints -> Append[ics, Automatic],
Epilog -> White, EdgeForm[Black], Disk[1, 0, 0.03]]
answered 7 hours ago
Michael E2Michael E2
153k12208493
$endgroup$
$begingroup$
this was what I was looking for
$endgroup$
– Javier
7 hours ago
add a comment |
$begingroup$
eqns = x'[t] ==
x[t] - y[t] -
x[t] (x[t]^2 + y[t]^2) + (x[t] y[t])/Sqrt[x[t]^2 + y[t]^2],
y'[t] ==
x[t] + y[t] - y[t] (x[t]^2 + y[t]^2) -
x[t]^2/Sqrt[x[t]^2 + y[t]^2];
sol = NDSolve[Join[x[0]==1.5, y[0]==1.5, eqns], x, y, t, 0, 50];
ParametricPlot[x[t], y[t]/.sol//Evaluate, t, 0, 50, PlotRange->All]
answered 8 hours ago
b3m2a1b3m2a1
29.9k360176
$endgroup$
add a comment |
Your Answer
StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "387"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);
else
createEditor();
);
function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);
);
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f199145%2fhow-to-plot-an-unstable-attractor%23new-answer', 'question_page');
);
Post as a guest
Required, but never shown
2 Answers
2
active
oldest
votes
2 Answers
2
active
oldest
votes
active
oldest
votes
active
oldest
votes
$begingroup$
To visualize a 2D system, I would start with StreamPlot:
vf = x', y' /. First@Solve[eqns /. f_[t] :> f, x', y']; (* strip the args *)
StreamPlot[vf, x, -2, 2, y, -2, 2]
You can use StreamPoints to highlight the structure and Epilog to mark the attractor at $(1,0)$:
ics = Cos[1/5], Sin[1/5], Red,
0.5, 0, Magenta, 1.5, 0., Magenta;
StreamPlot[vf, x, -2, 2, y, -2, 2,
StreamPoints -> Append[ics, Automatic],
Epilog -> White, EdgeForm[Black], Disk[1, 0, 0.03]]
answered 7 hours ago
Michael E2Michael E2
153k12208493
$endgroup$
$begingroup$
this was what I was looking for
$endgroup$
– Javier
7 hours ago
add a comment |
$begingroup$
To visualize a 2D system, I would start with StreamPlot:
vf = x', y' /. First@Solve[eqns /. f_[t] :> f, x', y']; (* strip the args *)
StreamPlot[vf, x, -2, 2, y, -2, 2]
You can use StreamPoints to highlight the structure and Epilog to mark the attractor at $(1,0)$:
ics = Cos[1/5], Sin[1/5], Red,
0.5, 0, Magenta, 1.5, 0., Magenta;
StreamPlot[vf, x, -2, 2, y, -2, 2,
StreamPoints -> Append[ics, Automatic],
Epilog -> White, EdgeForm[Black], Disk[1, 0, 0.03]]
answered 7 hours ago
Michael E2Michael E2
153k12208493
$endgroup$
$begingroup$
this was what I was looking for
$endgroup$
– Javier
7 hours ago
add a comment |
$begingroup$
To visualize a 2D system, I would start with StreamPlot:
vf = x', y' /. First@Solve[eqns /. f_[t] :> f, x', y']; (* strip the args *)
StreamPlot[vf, x, -2, 2, y, -2, 2]
You can use StreamPoints to highlight the structure and Epilog to mark the attractor at $(1,0)$:
ics = Cos[1/5], Sin[1/5], Red,
0.5, 0, Magenta, 1.5, 0., Magenta;
StreamPlot[vf, x, -2, 2, y, -2, 2,
StreamPoints -> Append[ics, Automatic],
Epilog -> White, EdgeForm[Black], Disk[1, 0, 0.03]]
answered 7 hours ago
Michael E2Michael E2
153k12208493
$endgroup$
To visualize a 2D system, I would start with StreamPlot:
vf = x', y' /. First@Solve[eqns /. f_[t] :> f, x', y']; (* strip the args *)
StreamPlot[vf, x, -2, 2, y, -2, 2]
You can use StreamPoints to highlight the structure and Epilog to mark the attractor at $(1,0)$:
ics = Cos[1/5], Sin[1/5], Red,
0.5, 0, Magenta, 1.5, 0., Magenta;
StreamPlot[vf, x, -2, 2, y, -2, 2,
StreamPoints -> Append[ics, Automatic],
Epilog -> White, EdgeForm[Black], Disk[1, 0, 0.03]]
answered 7 hours ago
Michael E2Michael E2
153k12208493
answered 7 hours ago
Michael E2Michael E2
153k12208493
answered 7 hours ago
Michael E2Michael E2
153k12208493
answered 7 hours ago
Michael E2Michael E2
153k12208493
153k12208493
$begingroup$
this was what I was looking for
$endgroup$
– Javier
7 hours ago
add a comment |
$begingroup$
this was what I was looking for
$endgroup$
– Javier
7 hours ago
$begingroup$
this was what I was looking for
$endgroup$
– Javier
7 hours ago
$begingroup$
this was what I was looking for
$endgroup$
– Javier
7 hours ago
add a comment |
$begingroup$
eqns = x'[t] ==
x[t] - y[t] -
x[t] (x[t]^2 + y[t]^2) + (x[t] y[t])/Sqrt[x[t]^2 + y[t]^2],
y'[t] ==
x[t] + y[t] - y[t] (x[t]^2 + y[t]^2) -
x[t]^2/Sqrt[x[t]^2 + y[t]^2];
sol = NDSolve[Join[x[0]==1.5, y[0]==1.5, eqns], x, y, t, 0, 50];
ParametricPlot[x[t], y[t]/.sol//Evaluate, t, 0, 50, PlotRange->All]
answered 8 hours ago
b3m2a1b3m2a1
29.9k360176
$endgroup$
add a comment |
$begingroup$
eqns = x'[t] ==
x[t] - y[t] -
x[t] (x[t]^2 + y[t]^2) + (x[t] y[t])/Sqrt[x[t]^2 + y[t]^2],
y'[t] ==
x[t] + y[t] - y[t] (x[t]^2 + y[t]^2) -
x[t]^2/Sqrt[x[t]^2 + y[t]^2];
sol = NDSolve[Join[x[0]==1.5, y[0]==1.5, eqns], x, y, t, 0, 50];
ParametricPlot[x[t], y[t]/.sol//Evaluate, t, 0, 50, PlotRange->All]
answered 8 hours ago
b3m2a1b3m2a1
29.9k360176
$endgroup$
add a comment |
$begingroup$
eqns = x'[t] ==
x[t] - y[t] -
x[t] (x[t]^2 + y[t]^2) + (x[t] y[t])/Sqrt[x[t]^2 + y[t]^2],
y'[t] ==
x[t] + y[t] - y[t] (x[t]^2 + y[t]^2) -
x[t]^2/Sqrt[x[t]^2 + y[t]^2];
sol = NDSolve[Join[x[0]==1.5, y[0]==1.5, eqns], x, y, t, 0, 50];
ParametricPlot[x[t], y[t]/.sol//Evaluate, t, 0, 50, PlotRange->All]
answered 8 hours ago
b3m2a1b3m2a1
29.9k360176
$endgroup$
eqns = x'[t] ==
x[t] - y[t] -
x[t] (x[t]^2 + y[t]^2) + (x[t] y[t])/Sqrt[x[t]^2 + y[t]^2],
y'[t] ==
x[t] + y[t] - y[t] (x[t]^2 + y[t]^2) -
x[t]^2/Sqrt[x[t]^2 + y[t]^2];
sol = NDSolve[Join[x[0]==1.5, y[0]==1.5, eqns], x, y, t, 0, 50];
ParametricPlot[x[t], y[t]/.sol//Evaluate, t, 0, 50, PlotRange->All]
answered 8 hours ago
b3m2a1b3m2a1
29.9k360176
answered 8 hours ago
b3m2a1b3m2a1
29.9k360176
answered 8 hours ago
b3m2a1b3m2a1
29.9k360176
answered 8 hours ago
b3m2a1b3m2a1
29.9k360176
29.9k360176
add a comment |
add a comment |
Thanks for contributing an answer to Mathematica Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
Use MathJax to format equations. MathJax reference.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f199145%2fhow-to-plot-an-unstable-attractor%23new-answer', 'question_page');
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
$begingroup$
Try using
NDSolveinstead$endgroup$
– b3m2a1
8 hours ago