strange attractor plot

up vote
2
down vote

favorite

s = NDSolve[{

    Derivative[1][x][t] == -y[t] - z[t],

    Derivative[1][y][t] == x[t] + 0.1 y[t],

    Derivative[1][z][t] == .01 + z[t] (x[t] - 14),

    x[0] == z[0] == 0, y[0] == 0}, {x, y, z}, {t, 0, 400},     MaxSteps -> [Infinity]]; Show[ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. s], {t, 0, 400},    PlotPoints -> 2000,    PlotStyle -> Directive[Thick, RGBColor[0, 0, 0]],    ColorFunction -> (ColorData["SolarColors", #1] &)],   Graphics3D[{ColorData["SolarColors"][0],     Sphere[First[({x[t], y[t], z[t]} /. s) /. t -> 0], 1]}],   RotationAction -> "Clip", Boxed

-> False, SphericalRegion -> False,   Axes -> False, ImageSize -> 500]

When I do this code it only shows a picture of part of the Rossler attractor https://en.wikipedia.org/wiki/R%C3%B6ssler_attractor, even though the derivatives are correct. What's going wrong?

asked Nov 17 at 17:24

Forever Mozart

1153

1

Add the PlotRange -> All option to Graphics3D.
– Rohit Namjoshi
Nov 17 at 17:29

@RohitNamjoshi still not working for me
– Forever Mozart
Nov 17 at 17:43

1

That is odd. What version of Mathematica are you running? I am on 11.3.0 for Mac OS X x86 (64-bit) (March 7, 2018). When I add that option I see this.
– Rohit Namjoshi
Nov 17 at 17:49

@RohitNamjoshi It is the latest version. Can you send me the complete code you entered? I may be doing something wrong (I'm new to this).
– Forever Mozart
Nov 17 at 17:56

I copied exactly what you posted and added PlotRange -> All as the last argument to Plot. Perhaps you have some previously bound symbol that is interfering. Try evaluating ClearAll["Global*"]` first.
– Rohit Namjoshi
Nov 17 at 18:12

|
show 2 more comments

up vote
2
down vote

favorite

s = NDSolve[{

    Derivative[1][x][t] == -y[t] - z[t],

    Derivative[1][y][t] == x[t] + 0.1 y[t],

    Derivative[1][z][t] == .01 + z[t] (x[t] - 14),

    x[0] == z[0] == 0, y[0] == 0}, {x, y, z}, {t, 0, 400},     MaxSteps -> [Infinity]]; Show[ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. s], {t, 0, 400},    PlotPoints -> 2000,    PlotStyle -> Directive[Thick, RGBColor[0, 0, 0]],    ColorFunction -> (ColorData["SolarColors", #1] &)],   Graphics3D[{ColorData["SolarColors"][0],     Sphere[First[({x[t], y[t], z[t]} /. s) /. t -> 0], 1]}],   RotationAction -> "Clip", Boxed

-> False, SphericalRegion -> False,   Axes -> False, ImageSize -> 500]

When I do this code it only shows a picture of part of the Rossler attractor https://en.wikipedia.org/wiki/R%C3%B6ssler_attractor, even though the derivatives are correct. What's going wrong?

asked Nov 17 at 17:24

Forever Mozart

1153

1

Add the PlotRange -> All option to Graphics3D.
– Rohit Namjoshi
Nov 17 at 17:29

@RohitNamjoshi still not working for me
– Forever Mozart
Nov 17 at 17:43

1

That is odd. What version of Mathematica are you running? I am on 11.3.0 for Mac OS X x86 (64-bit) (March 7, 2018). When I add that option I see this.
– Rohit Namjoshi
Nov 17 at 17:49

@RohitNamjoshi It is the latest version. Can you send me the complete code you entered? I may be doing something wrong (I'm new to this).
– Forever Mozart
Nov 17 at 17:56

I copied exactly what you posted and added PlotRange -> All as the last argument to Plot. Perhaps you have some previously bound symbol that is interfering. Try evaluating ClearAll["Global*"]` first.
– Rohit Namjoshi
Nov 17 at 18:12

|
show 2 more comments

up vote
2
down vote

favorite

s = NDSolve[{

    Derivative[1][x][t] == -y[t] - z[t],

    Derivative[1][y][t] == x[t] + 0.1 y[t],

    Derivative[1][z][t] == .01 + z[t] (x[t] - 14),

    x[0] == z[0] == 0, y[0] == 0}, {x, y, z}, {t, 0, 400},     MaxSteps -> [Infinity]]; Show[ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. s], {t, 0, 400},    PlotPoints -> 2000,    PlotStyle -> Directive[Thick, RGBColor[0, 0, 0]],    ColorFunction -> (ColorData["SolarColors", #1] &)],   Graphics3D[{ColorData["SolarColors"][0],     Sphere[First[({x[t], y[t], z[t]} /. s) /. t -> 0], 1]}],   RotationAction -> "Clip", Boxed

-> False, SphericalRegion -> False,   Axes -> False, ImageSize -> 500]

When I do this code it only shows a picture of part of the Rossler attractor https://en.wikipedia.org/wiki/R%C3%B6ssler_attractor, even though the derivatives are correct. What's going wrong?

asked Nov 17 at 17:24

Forever Mozart

1153

s = NDSolve[{

    Derivative[1][x][t] == -y[t] - z[t],

    Derivative[1][y][t] == x[t] + 0.1 y[t],

    Derivative[1][z][t] == .01 + z[t] (x[t] - 14),

    x[0] == z[0] == 0, y[0] == 0}, {x, y, z}, {t, 0, 400},     MaxSteps -> [Infinity]]; Show[ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. s], {t, 0, 400},    PlotPoints -> 2000,    PlotStyle -> Directive[Thick, RGBColor[0, 0, 0]],    ColorFunction -> (ColorData["SolarColors", #1] &)],   Graphics3D[{ColorData["SolarColors"][0],     Sphere[First[({x[t], y[t], z[t]} /. s) /. t -> 0], 1]}],   RotationAction -> "Clip", Boxed

-> False, SphericalRegion -> False,   Axes -> False, ImageSize -> 500]

When I do this code it only shows a picture of part of the Rossler attractor https://en.wikipedia.org/wiki/R%C3%B6ssler_attractor, even though the derivatives are correct. What's going wrong?

plotting

asked Nov 17 at 17:24

Forever Mozart

1153

asked Nov 17 at 17:24

Forever Mozart

1153

asked Nov 17 at 17:24

Forever Mozart

1153

asked Nov 17 at 17:24

Forever Mozart

1153

asked Nov 17 at 17:24

Forever Mozart

1153

1

Add the PlotRange -> All option to Graphics3D.
– Rohit Namjoshi
Nov 17 at 17:29

@RohitNamjoshi still not working for me
– Forever Mozart
Nov 17 at 17:43

1

That is odd. What version of Mathematica are you running? I am on 11.3.0 for Mac OS X x86 (64-bit) (March 7, 2018). When I add that option I see this.
– Rohit Namjoshi
Nov 17 at 17:49

@RohitNamjoshi It is the latest version. Can you send me the complete code you entered? I may be doing something wrong (I'm new to this).
– Forever Mozart
Nov 17 at 17:56

I copied exactly what you posted and added PlotRange -> All as the last argument to Plot. Perhaps you have some previously bound symbol that is interfering. Try evaluating ClearAll["Global*"]` first.
– Rohit Namjoshi
Nov 17 at 18:12

|
show 2 more comments

1

Add the PlotRange -> All option to Graphics3D.
– Rohit Namjoshi
Nov 17 at 17:29

@RohitNamjoshi still not working for me
– Forever Mozart
Nov 17 at 17:43

1

That is odd. What version of Mathematica are you running? I am on 11.3.0 for Mac OS X x86 (64-bit) (March 7, 2018). When I add that option I see this.
– Rohit Namjoshi
Nov 17 at 17:49

@RohitNamjoshi It is the latest version. Can you send me the complete code you entered? I may be doing something wrong (I'm new to this).
– Forever Mozart
Nov 17 at 17:56

I copied exactly what you posted and added PlotRange -> All as the last argument to Plot. Perhaps you have some previously bound symbol that is interfering. Try evaluating ClearAll["Global*"]` first.
– Rohit Namjoshi
Nov 17 at 18:12

Add the PlotRange -> All option to Graphics3D.
– Rohit Namjoshi
Nov 17 at 17:29

@RohitNamjoshi still not working for me
– Forever Mozart
Nov 17 at 17:43

That is odd. What version of Mathematica are you running? I am on 11.3.0 for Mac OS X x86 (64-bit) (March 7, 2018). When I add that option I see this.
– Rohit Namjoshi
Nov 17 at 17:49

@RohitNamjoshi It is the latest version. Can you send me the complete code you entered? I may be doing something wrong (I'm new to this).
– Forever Mozart
Nov 17 at 17:56

I copied exactly what you posted and added PlotRange -> All as the last argument to Plot. Perhaps you have some previously bound symbol that is interfering. Try evaluating ClearAll["Global*"]` first.
– Rohit Namjoshi
Nov 17 at 18:12

|
show 2 more comments

1 Answer
1

active

oldest

votes

up vote
5
down vote

accepted

Add PlotRange -> All

Show[ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. s], {t, 0, 400}, 

  PlotPoints -> 2000, 

  PlotStyle -> Directive[Thick, RGBColor[0, 0, 0]], 

  ColorFunction -> (ColorData["SolarColors", #1] &)], 

 Graphics3D[{ColorData["SolarColors"][0], 

   Sphere[First[({x[t], y[t], z[t]} /. s) /. t -> 0], 1]}], 

 RotationAction -> "Clip", Boxed -> False, SphericalRegion -> False, 

 Axes -> False, ImageSize -> 500, PlotRange -> All]

enter image description here

answered Nov 17 at 18:26

Rohit Namjoshi

63319

add a comment |

Your Answer

StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\$","\$"]]);
});
});
}, "mathjax-editing");

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',
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
});

}
});

draft saved

draft discarded

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Post as a guest

Name

Required, but never shown

StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f186198%2fstrange-attractor-plot%23new-answer', 'question_page');
}
);

Post as a guest

Name

Required, but never shown

1 Answer
1

active

oldest

votes

1 Answer
1

active

oldest

votes

up vote
5
down vote

accepted

Add PlotRange -> All

Show[ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. s], {t, 0, 400}, 

  PlotPoints -> 2000, 

  PlotStyle -> Directive[Thick, RGBColor[0, 0, 0]], 

  ColorFunction -> (ColorData["SolarColors", #1] &)], 

 Graphics3D[{ColorData["SolarColors"][0], 

   Sphere[First[({x[t], y[t], z[t]} /. s) /. t -> 0], 1]}], 

 RotationAction -> "Clip", Boxed -> False, SphericalRegion -> False, 

 Axes -> False, ImageSize -> 500, PlotRange -> All]

enter image description here

answered Nov 17 at 18:26

Rohit Namjoshi

63319

add a comment |

up vote
5
down vote

accepted

Add PlotRange -> All

Show[ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. s], {t, 0, 400}, 

  PlotPoints -> 2000, 

  PlotStyle -> Directive[Thick, RGBColor[0, 0, 0]], 

  ColorFunction -> (ColorData["SolarColors", #1] &)], 

 Graphics3D[{ColorData["SolarColors"][0], 

   Sphere[First[({x[t], y[t], z[t]} /. s) /. t -> 0], 1]}], 

 RotationAction -> "Clip", Boxed -> False, SphericalRegion -> False, 

 Axes -> False, ImageSize -> 500, PlotRange -> All]

enter image description here

answered Nov 17 at 18:26

Rohit Namjoshi

63319

add a comment |

up vote
5
down vote

accepted

Add PlotRange -> All

Show[ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. s], {t, 0, 400}, 

  PlotPoints -> 2000, 

  PlotStyle -> Directive[Thick, RGBColor[0, 0, 0]], 

  ColorFunction -> (ColorData["SolarColors", #1] &)], 

 Graphics3D[{ColorData["SolarColors"][0], 

   Sphere[First[({x[t], y[t], z[t]} /. s) /. t -> 0], 1]}], 

 RotationAction -> "Clip", Boxed -> False, SphericalRegion -> False, 

 Axes -> False, ImageSize -> 500, PlotRange -> All]

enter image description here

answered Nov 17 at 18:26

Rohit Namjoshi

63319

Add PlotRange -> All

Show[ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. s], {t, 0, 400}, 

  PlotPoints -> 2000, 

  PlotStyle -> Directive[Thick, RGBColor[0, 0, 0]], 

  ColorFunction -> (ColorData["SolarColors", #1] &)], 

 Graphics3D[{ColorData["SolarColors"][0], 

   Sphere[First[({x[t], y[t], z[t]} /. s) /. t -> 0], 1]}], 

 RotationAction -> "Clip", Boxed -> False, SphericalRegion -> False, 

 Axes -> False, ImageSize -> 500, PlotRange -> All]

enter image description here

answered Nov 17 at 18:26

Rohit Namjoshi

63319

answered Nov 17 at 18:26

Rohit Namjoshi

63319

answered Nov 17 at 18:26

Rohit Namjoshi

63319

answered Nov 17 at 18:26

Rohit Namjoshi

63319

add a comment |

draft saved

draft discarded

draft saved

draft discarded

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Post as a guest

Name

Required, but never shown

Post as a guest

Name

Required, but never shown

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Post as a guest

Name

Required, but never shown

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Post as a guest

Name

Required, but never shown

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Post as a guest

Name

Required, but never shown

Name

Required, but never shown

Name

Required, but never shown

This page is only for reference, If you need detailed information, please check here

搜尋此網誌

Gfrktyl