Question: The sample option in odeplot

Hello,

I've been trying to collect points at certain intervals in my plot, however, the sample option doesn't seem to be working. Can anyone lend me a hand? Any help is greatly appreciated! Thank you in advance!

Kind regards,

Gam

with(plots):

a := 1.501*10^9:

Th := sqrt(4*Pi^2*a^3/(G*(Mh+Msat))):

HyperionOrbit := proc (`θIC`, `ωIC`, n) local a, Mh, Msat, G, e, beta, M, Eqns, ICs, soln; global `ωH`, Th; a := 1.501*10^9; Mh := 5.5855*10^18; Msat := 5.6832*10^26; G := 6.67259/10^11; e := .232; beta := .89; M := Mh+Msat; Eqns := diff(theta(t), t) = omega(t), diff(omega(t), t) = -G*Msat*beta^2*(xH(t)*sin(theta(t))-yH(t)*cos(theta(t)))*(xH(t)*cos(theta(t))+yH(t)*sin(theta(t)))/(xH(t)^2+yH(t)^2)^2.5, diff(xH(t), t) = vxH(t), diff(vxH(t), t) = -G*M*xH(t)/(xH(t)^2+yH(t)^2)^(3/2), diff(yH(t), t) = vyH(t), diff(vyH(t), t) = -G*M*yH(t)/(xH(t)^2+yH(t)^2)^(3/2); ICs := xH(0) = a*(1+e), yH(0) = 0, vxH(0) = 0, vyH(0) = sqrt(G*M*(1-e)/(a*(1+e))), theta(0) = `θIC`, omega(0) = `ωIC`; soln := dsolve({Eqns, ICs}, numeric, maxfun = 0); plots:-odeplot(soln, [theta(t), omega(t)/`ωH`], 0 .. n*Th, numpoints = 2000, labels = ["θ(t)","ω(t)/ωH"], axes = boxed, style = point, sample = [seq(i, i = 0 .. n*Th, Th)]) end proc:

HyperionOrbit(.5, 1.8*`ωH`, 10)

Error, (in plot/options2d) unexpected option: sample = [0, 1876321.326, 3752642.652, 5628963.978, 7505285.304, 9381606.630, 11257927.96, 13134249.28, 15010570.61, 16886891.93, 18763213.26]

 

``

``

 

Download Poincare_section.mw

bia Man

Please Wait...