On Wed, Jan 15, 2014 at 11:33:58AM +, Oscar Benjamin wrote:
> On Wed, Jan 15, 2014 at 11:08:51AM +1100, Steven D'Aprano wrote:
> >
> > Unfortunately, I don't think there is anyone here who is very familiar
> > with scipy.
>
> There's Eryksun and me at least. I think that simple scientific Pyth
On Wed, Jan 15, 2014 at 11:08:51AM +1100, Steven D'Aprano wrote:
>
> Unfortunately, I don't think there is anyone here who is very familiar
> with scipy.
There's Eryksun and me at least. I think that simple scientific Python
questions get answered round here.
If the science itself or the librarie
On Tue, Jan 14, 2014 at 06:40:03PM +, Grace Roberts wrote:
> Hi,
> I'm a python beginner, currently using scipy's 'odeint' to compute to a set
> of ODEs (obtained by splitting Newton's law of gravity ma=-Gm1m2r/r^3 in two
> ordinary differentials). The idea is that I'm modelling the orbit of
On Wed, Jan 15, 2014 at 1:42 AM, eryksun wrote:
>
> The trajectory is decaying into the planet. In real life it hits the
> surface.
Not quite. A radius of 1.4 km is inside the planet, so that's
unrealistic from the start. If it starts at the surface of Mars, at
around 3,400 km, then the satellite
On Tue, Jan 14, 2014 at 1:40 PM, Grace Roberts
wrote:
>
> -For certain initial conditions the programme displays impossible orbits,
> showing the satellite making immediate sharp turns or jumping up and down in
> velocity. The exact same code can also end up displaying different graphs
> when run
On Tue, Jan 14, 2014 at 06:40:03PM +, Grace Roberts wrote:
[...]
> I've attached a quick print-screen of the code.
I've managed to have a look at the screenshot, there are no obvious
problems, but I'm not a Scipy expert. If you do get an answer from a
Scipy mailing list, I would really appre
Hi Grace,
I see that you're asking about the behavior of scipy.integrate.odeint():
http://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.odeint.html#scipy.integrate.odeint
Unfortunately this might be too domain specific of a question for
folks here at Tutor. I don't see anyth
Hi Grace, and welcome,
My reply is interleaved with your comments.
On Tue, Jan 14, 2014 at 06:40:03PM +, Grace Roberts wrote:
> Hi,
> I'm a python beginner, currently using scipy's 'odeint' to compute to
> a set of ODEs (obtained by splitting Newton's law of gravity
> ma=-Gm1m2r/r^3 in two
On Tue, Jan 14, 2014 at 1:40 PM, Grace Roberts
wrote:
> Hi,
>
> I'm a python beginner, currently using scipy's 'odeint' to compute to a set
> of ODEs (obtained by splitting Newton's law of gravity ma=-Gm1m2r/r^3 in two
> ordinary differentials). The idea is that I'm modelling the orbit of a
> sate
