Advanced Search

Question: Why isn't this plotting?

Question:Why isn't this plotting?

GambiaMan 15 Maple
ode dsolve
November 29 2015
0

Hello,

I have an issue with this code. The issue lies in the warning which is outputed. However, I don't know how to resolve it. Any help would be greatly appreciated! Thank you in advance!

Kind regards,

Gambia Man

> 

with(plots):
> 

a := 1.501*10^9:
> 

Eqns := diff(xH(t), t) = vxH(t), diff(vxH(t), t) = -G*(Mh+Msat)*xH(t)/(xH(t)^2+yH(t)^2)^(3/2), diff(yH(t), t) = vyH(t), diff(vyH, t) = -G*(Mh+Msat)*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*(Mh+Msat)*(1-e)/(a*(1+e))):
> 

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

soln := dsolve({Eqns, ICs}, numeric):
> 

plots:-odeplot(soln, [xH(t)/a, yH(t)/a], 0 .. 10*T, scaling = constrained, labels = ["x/a", "y/a"], numpoints = 2000);

Warning, cannot evaluate the solution further right of 453574.15, probably a singularity

  

 
> 

``

 

Download Hyperion_Orbit.mw

 
Please Wait...